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Abstract 

Mechanical systems are ideal candidates for studying quantum behavior of macroscopic objects. To this end, a mechanical 
■ resonator has to be cooled to its ground state and its position has to be measured with great accuracy. Currently, various 
routes to reach these goals are being explored. In this review, we discuss different techniques for sensitive position 
O ' detection and we give an overview of the cooling techniques that are being employed. The latter include sideband cooling 
I and active feedback cooling. The basic concepts that are important when measuring on mechanical systems with high 
accuracy and/or at very low temperatures, such as thermal and quantum noise, linear response theory, and backaction, 
are explained. From this, the quantum limit on linear position detection is obtained and the sensitivities that have been 
achieved in recent opto and nanoelectromechanical experiments are compared to this limit. The mechanical resonators 
that are used in the experiments range from meter-sized gravitational wave detectors to nanomechanical systems that can 
' ^ only be read out using mesoscopic devices such as single-electron transistors or superconducting quantum interference 
[ devices. A special class of nanomechanical systems are bottom-up fabricated carbon-based devices, which have very high 
^ • frequencies and yet a large zero-point motion, making them ideal for reaching the quantum regime. The mechanics of 
^ ' some of the different mechanical systems at the nanoscale is studied. We conclude this review with an outlook of how 
C state-of-the-art mechanical resonators can be improved to study quantum mechanics 
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1. Introduction 



Mechanics is probably the most well-known branch of physics as everyone encounters it in every-day life. It describes 
a wide range of effects: from the motion of galaxies and planets on a large scale, the vibrations of a bridge induced 
by traffic or wind, the stability of a riding bicycle, to the trajectories of electrons in an old-fashioned television on a 
microscopic scale. In the early days of physics, mainly objects that could be seen or touched were studied. Until the 
beginning of the twentieth century it was thought that the three laws of motion obtained by Newton described the dynamics 
of mechanical systems completely. However, the development of better telescopes and microscopes enabled the study of 
mechanical systems on both much larger and smaller length scales. In the early 1900s, the rapid developments that led to 
the theory of special and general relativity and quantum mechanics showed that the laws of classical mechanics were not 
the whole truth. 

Relativistic corrections turn out to be important for objects with large masses or with velocities approaching the speed 
of light. It is therefore an important factor in astrophysics, where one studies the dynamics of heavy objects like galaxies 
and black holes or the bending of light by the curvature of space-time. When the masses and velocities of the objects 
involved are made smaller and smaller, the relativistic corrections eventually vanish and one obtains the classical laws of 
motion 

Quantum mechanics, on the other hand, is particularly well suited to describe the mechanics of objects at the other end 
of the length-scale range, i.e., (sub)atomic objects. In the beginning of the twentieth century, quantum theory successfully 
explained the photo-electric effect, black-body radiation, and the atomic emission spectra. Quantum mechanics is different 
from classical and relativistic mechanics in the sense that objects are no longer described by a definite position, but by a 
wavefunction. This wavefunction evolves deterministic ally according to the Schrodinger equation and its absolute value 
squared should be interpreted as the position probability-density function, the so-called Born rule |Q. To find the object 
at a particular location one has to measure its position. This process, however, inevitably disturbs the evolution of the 
wavefunction [|^, ^] . 

Quantum mechanics does not only describe processes at the (sub)atomic scale successfully, but it also explains the 
microscopic origin of many macroscopic effects such as the electronic properties of solids, superfluidity and so on. Unlike 
in relativity where one can simply take the limit m, v — > 0, in quantum mechanics it is still not entirely clear how the 
transition from quantum mechanics to classical mechanics exactly happens Although Ehrenfest's theorem implies 

that the expectation values of the momentum and potential energy obey Newton's second law [^, the classical laws of 
motion cannot in all cases be recovered by simply taking /i — > as has recently been shown in calculating the quantum 
dynamics of nonlinear resonators [0]. 

Another issue that is still debated is how quantum mechanics should be interpreted [|[ ||]: as the truth, as a tool to 
calculate outcomes of an experiment, or as an incomplete theory? This question has been asked since the early days 
of quantum theory, resulting in the famous Einstein-Podolsky-Rosen paper [^, but has still not been answered. The 
interpretation of quantum mechanics and its transition to the macroscopic world are related and can be reformulated into 
the question "Can a macroscopic object be put in a quantum superposition?" This was first illustrated by Schrodinger 
in 1935 with the famous dead-or-alive cat gedanken experiment. Superpositions of small objects are readily observed, a 
good example are the singlet and triplet spin states in a molecule, but this becomes increasingly difficult for larger and 
larger systems mainly due to decoherence [ |l0| , pi] ]. Superpositions have for example been created using the circulating 
current in superconducting quantum interference devices (SQUIDs) JT^ and with fullerenes [ [l3| , 

Nanomechanical devices [|l5[ [l^] are interesting candidates to further increase the size of systems that can be put 
in a superposition i jlTI , [l^ ]. These systems are the logical continuation of micromechanical devices which are made 
using integrated-circuit technology, but then on a much smaller, nanometer scale. Micro-electromechanical systems 
(MEMS) are currently widely used as, for example, accelerometers in airbags, pressure sensors, and in projectors. When 
scaling these devices down to the nanometer scale, their resonance frequencies increase and at the same time their mass 
decreases. From an application point of view nano-electromechanical systems (NEMS) may enable single-atom mass- 
sensing, mechanical computing and efficient signal processing in the radio-frequency and microwave bands. From a 
scientific point of view these devices are interesting as they can be cooled to temperatures so low that the resonator is 
nearly always in its quantum-mechanical ground state. Figure[l| shows some examples of miniature mechanical resonators. 

In the last decade, many groups have pursued demonstration the quantum limit of motion. The rapid progress in 
the development of sensitive optical techniques and mesoscopic electronics have led to detectors that have sensitivities 
that are approaching the quantum limit on position detection [^. Moreover, very recently the long-sought quantum 
mechanical behavior of mechanical resonators has become reality in two different experiments indicating that the quantum 



regime of mechanical motion has now been entered: Selected by Science magazine as research breakthrough | |21[ | of 
the year 2010, the groups of Cleland and Martinis demonstrated quantum mechanical behavior of a 6 GHz mechanical 
resonator by coupling it to a superconducting flux qubit [|2| . At dilution refrigerator temperatures the first ten energy 
states of the harmonic mechanical oscillator, including the ground state, could be probed. The measurement itself is 
performed on the superconducting qubit which acts as a two-level quantum system whose response changes when the 
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Figure 1 : Different types of MEMS/NEMS used at Delft University of Technology, (a) Piezoresistive cantilever. A deflection of the cantilever changes 
the electrical resistance between the two sides. This resonator geometry can be used for massdetection. (b) Bistable buckled beam. Blue and red 
indicates the two stabl e, po sitions of the beam, which can be used to encode digital information [hsp. (c) A beam resonator that is embedded in the loop 
of a dc SQUID (red) | |l9| |. A magnetic field couples the position of the beam to the magnetic flux through the loop, (d) Suspended carbon nanotube 
(white) as a flexural resonator 02(1]. 



occupation of the mechanical resonator state changes. A large coupling between the two quantum systems is achieved 
by using the piezoelectric properties of the mechanical resonator material. In a different approach Teufel et al. use a 
superconducting microwave cavity to actively cool a mechanical drum resonator to such low temperatures that it is in the 
quantum mechanical ground state for most of the time. The drum resonator is integrated in the superconducting resonant 
circuit to provide strong phonon-photon coupling and the measurements use concepts developed for optical cavities to 
achieve efficient cooling. 

These two breakthrough experiments not only show that non-classical behavior can be encoded in the motion of a 
mechanical resonator but also open a new exciting research field involving quantum states of motion. In this review we 
will summarize the main theoretical and experimental discoveries that have led to the demonstration of quantum motion. 
Our main focus will be on the mechanics, i.e., on the mechanical properties of resonators and on the different optical and 
electronic detection schemes that have been developed to measure their displacements. Since measuring always means 
that the detector has to be coupled to the resonator, we will also discuss the coupling and the consequences it may have 
on the measurement itself. Furthermore, the advantages and disadvantages of the various approaches will be discussed. 



7.7. Ground-state cooling and quantum-limited position detection 

The concept of the quantum limit on mechanical motion detection and its implications became relevant in the 1970s 



when more and more sensitive gravitational wave detectors [ [24( | were designed (see for example Refs. [gSp and |26|), 
raising questions on the (im)possibility to violate Heisenberg uncertainty principle ^ |^, ^]. Nowadays, these issues are 
important when measuring micro- and nanomechanical devices with very sensitive detectors or at very low temperatures. 
There are two important considerations when approaching the quantum limit. First, the thermal occupation is impor- 

as n = (kBTR/hfR) - 1/2 where /« is the resonance frequency and Tr is the resonator 



tant, which is defined in Sec. 3.3 



temperature. A value for n that is below 1 indicates that the resonator is in its ground state most of the time (see Sec. 3.3 
for a more detailed description of the thermal occupation). Second, the zero-point motion sets the ultimate limit on the 
resonator displacement. At high temperatures and in the absence of actuation, Brownian motion determines the resonator 
position. As temperature decreases, the displacement decreases to the point that the quantum regime is reached (« < 1). 
Zero-point motion of the undriven beam remains with an amplitude of uq = h/Amnfn with m the resonator mass. To 
observe the zero-point motion, one has to detect it. This sounds trivial, but it turns out that a measurement on a quantum 
system inevitably disturbs it. Quantum mechanics sets a limit to the precision of continuously measuring the position of 
the resonator, the standard quantum limit. A detector at that limit is therefore called quantum-limited. Detection of the 
zero-point motion thus requires resonant frequencies higher than about 1 GHz (where n = 1 corresponds to a temperature 



of 50 mK, which can be reached in a dilution refrigerator) in combination with a position detection scheme that meets the 
quantum Hmit. Alternatively, a low-frequency resonator can be actively cooled to its quantum mechanical ground state. 
In experiments, both approaches are pursued and we discuss the status of both of them. 

1.2. Bottom-up and top-down nanomechanical devices 

Table |l] provides an overview of the different groups that have performed experiments with mechanical resonators that 
approach the quantum limit in position-detection sensitivity or that have been cooled to a low resonator temperature. The 
micro- and nanomechanical devices listed in this table are made using so-called top-down fabrication techniques, which 
are also employed in the semiconductor industry. Different groups use different types of resonators: doubly clamped 
beams, singly clamped cantilevers, radial breathing modes of silica microtoroids, membranes, micromirrors and macro- 
scopic bars. We have indicated the resonator mass, which ranges from 20 ag to 1000 kg, and the resonator frequency, 
which ranges from about 10 Hz to a few GHz. From these numbers the zero-point motion mq has been calculated; for 
top-down devices it is generally on the order of femtometers. 

These top-down nanoscale structures listed in Table |l| (e.g. beams, cantilevers and microtoroids) are made by etching 
parts of a larger structure, for example a thin film on a substrate, or by depositing material (evaporating, sputtering) on a 
resist mask that is subsequently removed in a lift-off process. In both cases, patterning of resist is needed, which is done 
using optical or electron-beam lithography. State-of-the-art top-down fabricated devices have thicknesses and widths of 
less than 100 nanometer. 

A major drawback of making smaller resonators to increase their frequency, is that the quality factor decreases [p^. 
The quality factor is a measure for the dissipation in the system. A low quality factor or "Q-factor" means a large 
dissipation, and this is an unwanted property for resonators in the quantum regime. The associated decoherence of 
quantum states then limits the time for performing operations with these states. For example, in the experiments of 
Ref. where Q is on the order of a few hundred, the time for the manipulation of quantum states is limited to only 6 ns. 
The decrease in Q-factor with device dimensions is often attributed to the increase in surface-to-volume ratio Jl^ |7l| , 
|72| . An explanation comes from the fabrication which may introduce defects at the surface during the micro-machining 
processes that are involved. These defects provide channels for dissipation, resulting in the low Q-factor. 

Having this in mind, a different approach is to use the small structures that nature gives us, to build or assemble 
mechanical resonators. Bottom-up devices are expected not to suffer from excessive damping, as their surface can be 
defect-free at the atomic scale. Examples are inorganic nanowires[| carbon nanotubes and few-layer graphene. The last 
two are examples of carbon-based materials. Using these bottom-up materials, mechanical devices with true nanometer 
dimensions can be made with the hope that surface defects can be eliminated. High Q-values are therefore expected for 
these devices, which, combined with their low mass, make them ideal building blocks for nano-electromechanical systems 

Table || shows the properties of mechanical resonators that have been made so far using bottom-up fabricated devices. 
Their frequencies are high: by choosing the right device geometry resonances in the UHF band (300 MHz - 3 GHz) are 
readily made, as Table || shows. When comparing the quaUty factors and zero-point motion of these devices, it is clear 
that the nanowire performance is more or less comparable to top-down fabricated devices since their thickness is of the 
order of 10-100 nm, about the size of the smallest top-down fabricated devices. 

Due to their low mass m and high strength (see the mechanical properties listed in Table in Sec. ^) the frequencies of 
carbon-based resonators are high and their zero-point motion uq large. Note, that in Table ^uq was given in femtometer, 
whereas in Table ^ it is given in picometer Figure ^ illustrates this point more clearly. The quantum regime with a large 
Mo and small n is positioned in the upper right corner When cooled to dilution refrigerator temperatures, carbon-based 
resonators would therefore be in the ground-state while exhibiting relatively large amplitude zero-point fluctuations. 

Position detectors for bottom-up NEMS, however, are not yet as sophisticated as those for the larger top-down coun- 
terparts. Consequently, non-driven motion at cryogenic temperatures (i.e., Brownian or zero-point motion) nor active 
cooling has been reported for carbon-based NEMS. The devices always need to be actuated to yield a measurable re- 
sponse. This is at least in part due to their small size, which makes the coupling to the detector small as well. Nevertheless 
impressive progress in understanding the electromechanical properties of bottom-up resonators has been made in recent 
years using so-called self-detecting schemes. In these schemes, the nanotube both acts as the actuator and detector of its 
own motion. Ultra-high quality factors have now been demonstrated for carbon nanotube resonators at low temperatures 
[|T]] as well as a strong coupling between electron transport and mechanical motion [ 106 ^ . 



1.3. Carbon-based materials 

In this Report we pay special attention to the mechanics of bottom-up materials as this has not been reviewed in such 
detail as the mechanics of the silicon-based top-down devices. In particular we will focus on the carbon-based materials 
as they have extraordinary mechanical and electrical properties. 



'Sometimes "nanowire" is also used for top-down fabricated devices. Here, tliis term is used exclusively for grown wires. 
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Table 1 : Overview of recent key experiments with micro- and nanomechanical resonators in chronological order. Several types of resonators and 
detection methods are used by different groups in the field and are measured at dift'erent temperatures T . The table shows the resonance frequency f^, 
quality factor Q and the mass m of the resonator. From this, the zero-point motion uq is calculated. 





Group 


Resonator 


fR(MHz) 


Q 


T(K) 


m(kg) 




uo (fm) 


Ref. 


1 


Roukes 


SiC beam 


1029 




500 


4.2 


3.4- 


10" 


17 


15 


1 


28 




2 


UCSB 


GaAs beam 


117 




1700 


0.030 


2.8 ■ 


10" 


15 


5.1 


1 


29 




3 


Schwab 


SiN/Au beam 


19.7 




35000 


0.035 


9.7 ■ 


10" 


16 


21 


1 


30 




4 


LMU 


Si/Au cantilever 


7.30- 10" 


_^ 


2000 


295 


8.6 ■ 


10" 


12 


12 


1 


31 




5 


Schwab 


SiN/Al beam 


21.9 




120000 


0.030 


6.8 ■ 


10" 


16 


24 


1 


32 




6 


LKB Paris 


Si micromirror 


0.815 




10000 


295 


1.9 


10" 


7 


7.3 ■ 10"^ 


1 


33 




7 


UCSB 


AFM cantilever 


1.25 ■ 10 


_2 


137000 


295 


2.4 ■ 


10" 


11 


5.3 


1 


34 




8 


Vienna 


Si02/Ti02 beam 


0.278 




9000 


295 


9.0 ■ 


10" 


12 


1.8 


I 


35 




9 


LKB Paris 


Si micromirror 


0.814 




10000 


295 


1.9 


10" 


7 


7.4- 10"^ 


1 


36 




10 


MPI-QO 


Silica toroid 


57.8 




2890 


300 


1.5 ■ 


10" 


11 


0.10 


1 


37 




11 


Roukes 


SiC/Au cantilever 


127 




900 


295 


5.0 ■ 


10" 


17 


36 


1 


38 




12 


LIGO 


Micromirror 


1.72 ■ 10" 


_4 


3200 


295 


1.0 


10" 


3 


7.0- 10"^ 


1 


39 




13 


LIGO 


Micromirror 


1.27 ■ 10" 


_5 


19950 


295 


1.0 


10" 


3 


0.026 


1 


40 




14 


LMU 


Si micromuTor 


0.547 




1059 


300 


1.1 ■ 


10" 


14 


37 


1 


41 




15 


JILA 


Gold beam 


43.1 




5000 


0.25 


2.3 ■ 


10" 


15 


9.2 


1 


42 




16 


IBM 


Si cantilever 


2.60- 10" 


_3 


55600 


2.2 


3.2 ■ 


10" 


13 


100 


1 


43 




17 


NIST 


Si cantilever 


7.00 ■ 10" 


_3 


20000 


295 


1.0 ■ 


10" 


10 


3.5 


1 


44 




18 


LKB Paris 


Micromirror 


0.711 




16000 


295 


7.4 


10" 


4 


1.3 ■ 10""* 


1 


45 




19 


Harris 


SiN membrane 


0.134 




1100000 


294 


3.9 ■ 


10" 


11 


1.3 


1 


46 




20 


Vienna 


Si cantilever 


0.557 




2000 


35 


4.0 ■ 


10" 


11 


0.61 


1 


47 




21 


ANU 


Mirror on beam 


8.48 ■ 10" 


_5 


44500 


300 


6.9 


10" 


4 


0.10 


1 


48 




22 


MPI-QO 


Silica toroid 


74.0 




57000 


295 


1.0 ■ 


10" 


11 


0.11 


1 


49 




23 


JILA 


Al beam 


0.237 




2300 


0.040 


2.0 ■ 


10" 


15 


133 


1 


50 




24 


Roukes 


SiC/Au beam 


428 




2500 


22 


5.1 ■ 


10" 


17 


20 


1 


51 




25 


IBM 


Si cantilever 


4.95 ■ 10 


_3 


22500 


4.2 


2.0 ■ 


10" 


12 


29 


1 


52 




26 


Delft 


AlGaSb beam 


2.00 




18000 


0.020 


6.1 ■ 


10" 


13 


2.6 


1 


19 




27 


AURIGA 


Al bar 


8.65 ■ 10" 


_4 


1200000 


4.2 


1.1 


■ 10 




3.0- 10"« 


1 


53 




28 


AURIGA 


Al bar 


9.14 ■ 10" 


_4 


880000 


4.2 


1.1 


■ 10 




2.9 ■ 10"^ 


1 


53 




29 


Alberta 


Si cantilever 


1040 




18 


295 


2.0- 


10" 


17 


20 


1 


54 




30 


Tang 


Si beam 


8.87 




1850 


295 


1.3 ■ 


10" 


15 


27 


1 


55 




31 


JILA 


Al beam 


1.53 




300000 


0.050 


6.2 ■ 


10" 


15 


30 


1 


56 




32 


JILA 


Al beam 


1.53 




10000 


0.050 


6.2 ■ 


10" 


15 


30 


1 


56 




33 


LMU 


SiN beam 


8.90 




150000 


295 


1.8 ■ 


10" 


15 


23 


1 


57 




34 


Queensland 


Silica toroid 






545 


300 






8 


6.7 ■ 10"^ 


1 
1 


JO 




35 


Tang 


Si cantilever 


13.86 




4500 


295 


4.5 ■ 


10" 


16 


37 


1 


59 




36 


Vienna 


Si cantilever 


0.945 




30000 


5.3 


4.3 ■ 


10" 


11 


0.45 


1 


60 




37 


MPI-QO 


Silica toroid 


65.0 




2000 


1.65 


7.0- 


10" 


11 


0.043 


1 


61 




38 


Painter 


Si ph. crystal 


8.2 




150 


360 


4.3 ■ 


10" 


14 


4.9 


1 


62 




39 


Oregon 


Silica sphere 


118.6 




3400 


1.4 


2.8 ■ 


10" 


11 


0.050 


1 


63 




40 


LIGO 


Susp. mirror 


1.23- 10" 


-4 




300 


2.7 


■ 10 




1.6- 10"4 


1 


64 




41 


Painter 


Silica double tor. 


8.53 




4070 


300 


1.5 ■ 


10" 


13 


2.6 


1 


65 




42 


JILA 


Al beam 


1.04 




160000 


0.015 


1.1 ■ 


10" 


14 


27 


1 


66 




43 


MPI-QO / LMU 


SiN beam 


8.07 




10000 


300 


4.9- 


10" 


15 


15 


1 


67 




44 


Schwab 


SiN/Al beam 


6.30 




1000000 


0.020 


2.1 ■ 


10" 


15 


25 


1 


68 




45 


MPI-QO / LMU 


SlN beam 


8.30 




30000 


300 


3.7- 


10" 


15 


16 


1 


69 




46 


UCSB 


AIN FBAR 


6170 




260 


0.025 


2.8 ■ 


10" 


12 


0.022 


1 


22 




47 


Delft 


AlGaSb beam 


2.14 




24000 


0.015 


6.1 ■ 


10" 


13 


2.5 


1 


70 




48 


NIST 


Al drum 


10.69 




360000 


0.020 


4.8- 


10" 


14 


4.0 


1 


23 
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Table 2: Overview of recent experiments with bottom-up resonators. Several types of resonators are used: carbon nanotubes (CNT), nanowires (NW), 
single-layer graphene (SLG), and few-layer graphene (FLG) or graphene oxide (FLGO) sheets. The table shows the resonance frequency fa and quality 
factors at room and cryogenic temperature (Qrt and Qcryo resp.). T^\^ is the lowest temperature at which the resonator is measured, m is the mass of 
the resonator and C is its length. From these data, the zero-point motion uq is calculated. 







fR 






Train 


I 


m 




Un 




Group 


Type 


(MHz) 


Qrt 


Qcryo 


(K) 




(kg) 




(pm) 


I 


lef 




Cornell 


I^IN 1 


j3 






^nn 


1 7^^ 


7.4 


10 


21 


4.52 [ 


73 


1 


Ucilt 


L,iN 1 




1 nn 

lUU 




^nn 

jUU 


1 l'^ 


1.0 


10 


20 


3.66 [ 


20 


1 


Berkeley 


PMT 
1 




Aztn 






U.J 


5.3 


10 


20 


0.67 [ 


74 


1 




PMT 




o 
o 








5.8 


10- 


20 


0.22 [ 


75 


1 






1 


9n 
zu 




jUU 


U.ZOD 


1.1 


10- 


19 


0.69 [ 


75 


1 






J 1 J 










4.0 


10- 


22 


6.03 [ 


75 


1 




1 


ID / 


ZUU 


zuuu 


J 


n Q 


1.4 


10 


21 


6.03 [ 


76 


1 


Berkeley 


1 




1 000 






U.ZUJ 


1.0 


lU 


21 


/I n 1 r 
4.UI 1 


u. 


1 


rnlTpnh 
^^tli iCLll 


CNT 








\j 




4.8 


IQ- 
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Figure 2: Comparisnii nf thgresonance frequency and y prn-n nint motion hq of top-down (Table |l|) and two kinds of bottom-up devices: Suspended 
carbon nanotubes fcol Rll, |l06| | and graphene nanodrums | |l07| | . The resonance frequencies and zero-point motion are much larger for the bottom-up 
devices. 




diamond 



graphite 



carbon nanotube 



Figure 3: The structure of the different allotropes of carbon, (a) Diamond has two intertwined face-centered cubic lattices, (b) Graphite consists of 
stacked planes of hexagonally ordered carbon atoms. A single plane is called a "graphene sheef . (c) A Cgo buckyball molecule, (d) A single-walled 
carbon nanotube, which can be viewed as a graphene sheet that has been rolled up and sewn together, (c) and (d) Reprinted from Nanomedicine: Nan- 
otechnology. Biology and Medicine, 4 3, M. Foldvari, M. Bagonluri, Carbon nanotubes as functional excipients for nanomedicines: I. pharmaceutical 
properties, 173-182, Copyright (2008), with permission from Elsevier. 



Carbon exists in many different forms, ranging from amorphous coal to crystalline graphite and diamond. Diamond 
has a face-centered cubic structure as shown in Fig. ^a and is one of the hardest materials known. Its Young's modulus 
(Table H) is extremely high: about 1 TPa. Graphite has a very different crystal structure: it consists of stacked planes of 
carbon atoms in a hexagonal arrangement (Fig. ||b). Its Young's modulus for in-plane stress is nearly as high as that of 
diamond, but it is much lower for out-of -plane stress, as Table ^ will indicate. This difference is caused by the nature of 
the bonds holding the carbon atoms together. Atoms in one of the planes are covalently bonded to each other, whereas 
different planes are held together by the much weaker van der Waals force. 

Graphite and diamond were already known for millennia, but in the last decades novel allotropes o f car bon were 
discovered. First, in 1985 C6o molecules, called B uckm insterfullerenes or "buckyballs", were synthesized [108]. Then in 
the early 1990s carbon nanotubes were discovered [109|. These consist of cylinders of hexagonally ordered carbon atoms; 
similar to what one would get if one were to take a single layer of graphite and role it up into a cylinder In 2004 another 
allotrope called "graphene" was identified [110, 1 11 1. This is a single layer of graphite, which is, unlike a nanotube, flat. 
Graphene is usually deposited onto a substrate using mechanical exfoliation, and high-quality sheets of mm-size have 
been made using this technique [112]. Another way of making graphene devices is to grow it directly on a substrate, 
and in a semi-industrial process meter-sized sheets have been reported [113]. Although truly two-dimensional structures 
are not energetically stable ]114], graphene can exist due to fact that it contains ripples that stabilize its atomically thin 
structure J 1 15 116,117,118]. Figure ^shows the structure of the four different carbon allotropes. 

Although nanomechanical devices have been made out of diamond using top-down fabrication techniques [119, 12C , 
[tT , 121], in this Review we focus on bottom-up fabricated carbon-based NEMS made from few-layer graphene or sus- 
pended carbon nanotubes. 
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1.4. Outline of this review 

This Review consists of this introduction, three main Sections, one on nanomechanics, one on backaction and cooling, 
and the last one focusses on different types of detectors. These are supplemented with a final Section summarizing some 
prospects and future directions in the field of quantum electromechanical systems or QEMS in short. In Section || we will 
discuss mechanics at the nanoscale. From continuum mechanics, the general equations of motion are obtained. These 
are illustrated by studying the dynamics of a number of nanomechanical devices, such as beam and string-like resonators, 
buckled beams and suspended carbon nanotubes. We will demonstrate that the dynamics of a particular vibrational mode 
of the resonator is that of a harmonic oscillator. Another important point is that in nanoscale devices, tension is a crucial 
property that must be included in the analysis. As we will show it also provides a unique tool for the tuning of resonators 
when their thickness is on the atomic to nanometer scale. 

In Section ^ we discuss the properties of the (quantum) harmonic oscillator and we study the effects of backaction. A 
measurement always influences the measured object itself and this is called backaction. It has important consequences for 
linear detection schemes, and we will show that backaction ultimately limits their position resolution to what is known 
as the standard quantum limit. This limit will be derived in several ways. Backaction, however, can also be used to 
ones advantage as it provides a way to cool resonator modes. This can be done by backaction alone (self-cooling). 
Alternatively, the resonator temperature can be lowered by adopting active cooling protocols. An overview of the two 
most popular protocols will be given, including a summary of the main achievements. 

In general, two distinct approaches for position detection of mechanical oscillators are used: optical and electrical. In 
Sec. Q we give an overview of the different detection schemes. We will discuss the use of optical cavities in various forms, 
optical waveguides and their analogues in solid state devices in which electrons play the role of photons. In addition, we 
will summarize the concepts behind capacitive and inductive actuation and detection as well as the self-detection schemes 
that are used in bottom-up NEMS devices. Special attention will be paid to the achieved position resolution and to the 
limitations that prevent the standard quantum limit to be reached. Furthermore, we explain the mechanisms behind the 
backaction and, when possible, quantify the coupling between resonator and detector. 
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2. Mechanics at the nanoscale 



Nearly, all objects, including macroscopic and nanometre-sized systems, have particular frequencies at which they 
can resonate when actuated. These frequencies are called the eigenfrequencies or normal frequencies; when actuated at 
an eigenfrequency, i.e. on resonance, the amplitude of vibration can become very large and this can have catastrophic 
consequences. Examples include the breaking of a glass by sound waves and the collapse of the Tacoma bridge which 
was set in motion by the wind flowing around it. 

Continuum mechanics provides the tools to calculate these frequencies. For simple geometries such as cantilevers, 
doubly-clamped beams or thin plates, the frequencies of flexural or torsional modes can be calculated analytically. How- 
ever, in many NEMS experiments more complicated geometries are used, such as microtoroids [122|, suspended photonic 



crystal structures [52 123 1, or film bulk acoustic resonators (FBARs) ||22, 124] and the vibrational mode may have more 



a complicated shape than the simple flexural or torsional motion. To calculate the eigenfrequencies and mode shapes in 
these cases one has to rely on numerical calculations. Popular implementations include the software packages ANSYS 
and COMSOL which are based on the finite-element method. These packages have additional advantages as they provide 
a means to model other properti es of the nanomechanical system: This includes electrostati c int eractions (see e.g. Refs. 
[ 125 , 126, 5T^), thermal effects [127, 122|, and electromagnetic (optical) properties [^ 128 |. 

In this Section, we will summarize some main results of continuum mechanics needed to describe the experiments 
discussed in this Report. We will review analytical expressions for the eigenfrequencies of simple structures such as 
beams, buckled beams and strings. This material has been described in several textbooks (see e.g. the book by A. Cleland 
[[l5|]). Less attention has been devoted to thin beams or plates with nanometer-sized cross sections made from, for example, 
carbon nano tubes or graphene. These resonators are in a different regime than top-down devices that generally have larger 
sizes. In particular, the deflection of these nano-resonators can exceed their thickness or radius so that tension-induced 
nonlinear eff'ects start to play a role. We derive the equations for describing these nonlinear eff'ects in nanobeams and 
show that the induced tension can be used to tune the frequencies over a large range. We end this section with a discussion 
on the mechanics of (layered) graphene resonators that can be viewed as miniature drums. 



2.1. Continuum mechanics 

To describe the motion of mechanical objects, the dynamics of all particles (i.e., atoms and electrons) which make up 
the oscillator should, in principle, be taken into account. It is, however, known that for large, macroscopic objects this is 
unnecessary and that materials can be accurately described as a continuum with the mechanical behavior captured by a few 
parameters such as the elasticity tensor. Molecular dynamics simulations [129, 130, 131, 132, 133, 134| and experiments 
[ pO , 135, 107] demonstrate that even for nanometer-sized objects continuum mechanics is, with some modifications, still 
applicable. This means that the dynamics of the individual particles is irrelevant when one talks about deflections and 
deformations; the microscopic details do, however, determine the material properties and therefore also the values of 
macroscopic quantities like the Young's modulus or the Poisson ratio. 

The basis of continuum mechanics lies in the relations between strain and stress in a material. The strain tells how the 
material is deformed with respect to its relaxed state. After the deformation of the material, the part that was originally 
at position x is displaced by u to its new location x + u. The strain describes how much an infinitesimal line segment is 
elongated by the deformation u(jt:, y, z) and is given byfl [|136[ : 



1 / dui du j du„ du„ 



2 \ dxj dxi dxi dxj 



(1) 



This definition shows that strain is symmetric under a reversal of the indices, i.e., - jjj. The diagonal elements (i.e., 
/ = i) of the first two terms are the normal strains, whereas the off-diagonal elements (/ + j) are the shear strains. Eq. ^is 
exact, but the last, non-linear term is only relevant when the deformations are large [[ij] and will not be considered in this 
work. 

To deform a material external forces have to be applied, which in turn give rise to forces inside the material. When the 
material is thought of as composed of small elements, each element feels the force exerted on its faces by the neighboring 
elements. The magnitude and direction of the force depend not only on the location of the element in the material but also 
on the orientation of its faces (see Fig. Bp). The force 6F on a small area 6A of the element is given by: 



6Fi 



CTijnjdA, 



(2) 



■^In this Report, the so-called Einstein notation 1 136] for the elements of vectors and tensors is employed. When indices appear on one side of an 
equal sign only, one sums over them, without explicitly writing the summation sign. For example, .v, = RjjXj reads as x, = Rij^j- The index runs 
over the three cartesian coordinates (x,y,z), where xi = x, X2 = y and X3 = z- Finally, the symbols x, (x, y and z) denote the unit vectors in the fixed 
rectangular coordinate system (which form a basis) so that a vector r can be expressed as r = r/Xj. 
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Figure 4: (a) Visualization of the stress tensor on a cubic element AV. The force per unit area is the inner product of the stress tensor and the normal 
vector of the surface n. (b) Deformation of a plate under plane stress. The original plate (dotted) is deformed by the stress cTyy. 



where n is the vector perpendicular to the surface and cr is the stress tensor Now consider an element of the material with 
mass Am and volume AV. When it is moving with a speed v = u, its momentum Ap is: 



Ap 



J Am 



V dm •■ 



f p\dV, 

Jav 



(3) 



where p is the mass density. The rate of change of momentum equals the sum of the forces working on the element. These 
forces include the stress cr at the surface and body forces Fb that act on the volume element AV: 



— = r FbdV + r o-dA. 

df Jav Js(av) 



(4) 



Examples of body forces are gravity, with Fb = pgz and the electrostatic force q&, where g is the gravitational acceleration, 
q is the charge density, and £ is the local electric field. Using the Green-Gauss theorem, the integral over the boundary of 
the element can be converted into an integral over the volume: J^^j^y^ crdA - J^^ dcrij/dxi ■ x, dV. Eq. ^ should hold for 
any element because so far nothing has been specified about the shape or size of the element. This then yields Cauchy's 



first law of motion [136]: 



puj 



dxi 



(5) 



A similar analysis for the angular momentum yields Cauchy's second law of motion: cr,j = ctj,-. With these equations (and 
boundary conditions) the stress distribution can be calculated for a given applied force profile ^b(x,y,z). 



2.2. Elasticity 

The stress tensor gives the forces acting inside the material, whereas the strain tensor describes the local deformation 
of the material. These two quantities are, of course, related to each other When the deformations are not too large, the 
stress and strain tensor are related linearly via the the fourth-rank elasticity tensor E: 



Eijkijki 



(6) 



The properties of the stress and strain tensor imply that Eiju - Ejiu - Eijik - Ekuj, so E has at most 21 independent 
elements out of a total of 3x3x3x3 = 81 elements. This makes it possible to express Eq. ^ in a convenient matrix 
representation: 



a- XX 




CTyV 








O-xz 




O-yz 




_ O-xx 





ExX.XX Ex_xyy E XXZZ ^ XXXZ. 

ExXVV Eyyyy E yy^y E yyx^ 



F 



^ XXZZ 

Exxxz ^yyxz ^zzzx ^xzzx 

E xxyz ^yyyz ^ zzzy ^ xzzy 

ExXXV ^ \'V\'X ^ZZXY E yxxZ 



^xxyz 
^ yyyz 



'-xzzy 
p 

^ \xxz 



p 



E'xyyz ^xyyx 



7xx 

yyy 

Jzz 

2rxv 



(7) 



or in short hand notation^: [cr] = [£][y]. The inverse of the elasticity tensor is called the compliance tensor C, which 
expresses the strain in terms of the stress: 



Jij = Cijkio-u,oi [y] = [C][cr]. 



(8) 



The number of independent elements of E is further reduced when the crystal structure of the material has symmetries 
[O, 137, 138]. The most drastic example is an isotropic material, whose properties are the same in all directions. In this 



'Note that there are three different notations for the elasticity tensor: E is the actual tensor with elements Eij^. Finally, the elements can also be 
written in a matrix [E\ . 
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Table 3: Mechanical properties of materials that are used in nanomechanical devices. Most materials have a density p around 3 ■ 10^ kg/m^^ and a 
Young's modu lus of the order of 10^ GPa. The carbon-based materials graphite and diamond are slightly lighter, but much stiffer. Compiled from Refs. 
^ and Q. 



Material 


p (lO^ke/m^) 


E (GPa) 


V 


Silicon 


2.33 


130.2 


0.28 


Si3N4 


3.10 


357 


0.25 


SiC 


3.17 


166.4 


0.40 




2.65 


85.0 


0.09 


SiOa (amorphous) 


2.20 


~ 80 


0.17 


Diamond 


3.51 


992.2 


0.14 


Graphite (in-plane) 


2.20 


920 


0.052 


Graphite (out-of -plane) 


2.20 


33 


0.076 


Aluminum 


2.70 


63.1 


0.36 


Gold 


19.30 


43.0 


0.46 


Platinum 


21.50 


136.3 


0.42 


Niobium 


8.57 


151.5 


0.35 


GaAs 


5.32 


85.3 


0.31 


InAs 


5.68 


51.4 


0.35 



case, only two independent parameters remain: the Young's modulus E and Poisson's ratio v. The compUance matrix is 
in this case given by: 



l/E 


-v/E 


-v/E 











-v/E 


l/E 


-v/E 











-v/E 


-v/E 


l/E 




















1/G 




















1/G 




















1/G _ 



[C] = 



where G is the shear modulus. By inverting [C], the elasticity matrix is obtained: 



G = 



2h-2v 



(9) 
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Ev 











Ev 


Ev 


E(l - V) 
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G 



(10) 



When a plane stress cr„ is applied, a material will be stretched in the x-direction, as illustrated in Fig. The resulting 
strain - cr.Kx/E induces stress in the y- and z-directions, which are nulled by a negative strain (i.e., contraction) in 
these directions that is v times smaller than the strain in the x-direction. This follows directly from the structure of 
the compliance matrix. In the opposite situation where a plane strain is applied, the stress can directly be calculated 
using the elasticity matrix, Eq. 10 [136|. The Young's modulus and Poisson's ratio of materials that are frequently used 
for nanomechanical devices are indicated in Table ||. The Young's modulus of most semiconductors (Si, GaAs, InAs), 
insulators (Si3N4, SiC, SiOa) and metals is of the order of 100 GPa. This is much larger than the values for soft materials 
such as polymers (typically between 0.1 and 1 GPa) which have also been used for nanomechanical devices [139], but 
still smaller than that of diamond and graphite. Their Young's modulus of slightly less than 1 TPa combined with a low 
mass density makes these carbon-based materials ideal to build high-frequency resonators. Also, the large spread in the 
mass density of the metals should be noted. 

For non-isotropic materials, the Young's modulus and Poisson's ratio depend on the direction of the applied stress 
and are defined as: = l/Cuu and v/j = -Cajj/Caa (/ j) [141 1. Graphite is highly anisotropic, and its mechanical 



properties are important for studying many carbon-based nanomechanical devices. It consists of la yers of carbon atoms 
that are stacked on top of each other (see Fig. ||b and |]b) with an inter-layer spacing c = 0.335 nm [142|. The individual 
layers are called graphene sheets. The unit cell of graphene consists of a hexagon with a carbon atom on each corner, as 
illustrated in Fig. Each of the six carbon atoms lies in three different unit cells; a single unit cell thus contains two 
carbon atoms. The sides of the hexagon have a length c/„ = 0.14nm, so that the unit cell has an area of 5.22 ■ 10"^*^ m^, 
and the two-dimensional mass density is p2d = 6.8 ■ 10"^ kg/m^ . The graphene planes are not located exactly above each 
other but every other layer is shifted by half the unit cell, or equivalently, it is rotated by 60° around an axis through one 
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Figure 5: (a) The unit cell of graphene with the dimensions indicated, a is the length of the two translation vectors ai 2 



ia[±l, V3,0]. and rf„ is the 



distance between two carbon atoms. The area of the u nit ce ll is \ V3a^ = 5.22 ■ 10 m^. (b) Bending of a few-layer graphene sheet. The equilibrium 



distance between the graphene layers is c = 0.335 nm [142| 



of the carbon atoms. Three of the six atoms are on top of the atoms in the other layer and the other three are located at 
the center of the hexagon below them. The six-fold rotational symmetry ensures that the elast i c pro perties are the same 
when looking in any direction along the planes, i.e., they are isotropic in those directions [137, 138]. On the other hand, 
the mechanical properties for deformations perpendicular to the planes are quite different. It is therefore convenient to 
introduce the in- and out-of-plane Young's modulus, and respectively, and the corresponding Poisson's ratios Vq 
and Vj.. They are defined such that the compliance matrix is given by^ 



[C] 



!/£□ 


-yulEu 


-vJE^ 











-VulEn 


!/£□ 


-vJE^ 











-vJE^ 


-yJE^ 


IIE^ 




























































1/G^ 



(11) 



where|j£a = 0.92 TPa, E^ = 33GPa, Gq = 1.8 GPa, = 0.44 TPa, Vq = 0.052 and = 0.076 [|137|]. The large 
in-plane stiffness is the reason that carbon nanotubes and graphene have very high Young's moduli of about 1 TPa, which 
makes them one of the strongest materials known. Note, that these six elastic constants are not independent as the in-plane 
shear modulus is given by = (£1111 - £'ii22)/2 for a material with hexagonal symmetry. 



2.3. Energy, bending rigidity and tension 

In the previous Section the relation between the stress and strain in a material was given. Here, we focus on the energy 
needed to deform the material. From this, the equations of motion are derived. For small deformations, the potential 
energy U depends quadratically on the strain. It should be invariant under coordinate transformations [136|, leading to 



U 



= [/' dV with U' = ^-Eijuyijyu, 



(12) 



where U' is the potential energy den sity. From this definition it follows that the stress is given by cr,^ = dU' /dyij. For an 
isotropic material Eq. O reduces to [138|: 



U' ^- 



1 E 



2 1 + V 



(13) 



Although Eq. [T^ is valid for any mechanical system in the Unear regime, it is not straightforward to analyze a system this 
way. Therefore, we focus on two simple and frequently used geometries where the equation of motion can be obtained 
without too much effort, namely plates and beams. 



Note that this definition is slightly different from the conventional definition of the Poisson's ratio in an anisotropic material that was given on page 

The values of the elastic constants ri pxienri on the giiality of the graphite samples. Therefore, slightly diff'erent values can be found in the literature. 
Compare, for example, the data in Refs. ||l43|| and [ 144 1 with the values in Ref. ||l37|] 
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(a) 




(c) 




Figure 6: (a) Bending of a plate witli tliickness h. The top part of the plate is extended, whereas the bottom is compressed. The black dashed line 
indicates the neutral plane, (b) The normal vector of the top surface of a slightly deflected plate makes a small angle with the unit vector z (dotted), (c) 
If tension T is present in a plate, a net vertical force F results when the displacement profile has a finite curvature. 



2.3.1. Plates 

A plate is a thin object that is long and wide, i.e., h ^ € < w. When a torque is applied to it, it bends, as illustrated 
in Fig. The top part, which was initially at z = hjl, is extended whereas the bottom part of the beam, originally at 
z - -hjl, is compressed. There is a plane through the plate where the longitudinal strain is zero: the so-called neutral 
plane. T he ve rtical displacement of this plane is indicated by u{x, y) and for small deflections it lies midway through 
the plate |138], which we take at z = 0. Because of the small deflection and small thickness this is called the thin-plate 
approximation. 

Consider the top (or bottom) face of the plate that is shown in Fig. |6|b: because there is no material above (below) that 
surface, there cannot be a normal force F„ at this surface (except at the clamping points). In other words, the perpendicular 
stress components vanishes: F„ , - crijiij - 0. For a thin plate, the normal vector n at the top and bottom face points in the 
z-direction i.e., n = +z to first order in the displacement u or, equivalently, in the radius of curvature 7?^'. The condition 
for vanishing stress thus becomes: cr„ - a-y, - cr^ = at the faces. This not only holds at the faces but also inside the 
material because the plate is thin and the stress cannot build up. For an isotropic material the displacement and strain 
fields that satisfy these requirements are [ |138[ ]: 

-z ■ du/dx, y„ = -z ■ d^u/dx^ 
-z ■ du/dy, yyy = -z ■ d^u/dy^ 

Yzz - zy/(l - v) ■ V^M 
0, y.vv = ■ d^ujdxdy. 



Ux = 

Uy - 

u, — 



7.XZ 



Jyz 



(14) 



The vertical displacement of the material u, is thus equal to the deflection u for every value of z. Besides this vertical 
displacement there is also a horizontal displacement (mj and m,.) induced when u changes. In that case, the material 
displaces in different directions above and below the neutral plane (Fig. ^a) as expressed by the proportionality with z. 
Moreover, the strain components averaged over the thickness h, Jijix, y), are all zero, as the contributions above and below 
the neutral plane cancel each other. 

To proceed, Eq. |l^ is inserted into Eq. |l3| and the integration over z in Eq. |2]is carried out. This yields the energy 
needed to bend the plate: 



Eh' 



24(1 



v2)// 



dx^ dy^ 



\ dxdy J 



dx?- dy^ ^ 



Ax Ay 



(15) 



It is possible that the plate is not only bent, but that it is also under a longitudinal tension T - /cdz (positive for 
tensile tension, negative for compressive tension). The tension is tangential to the surface and from Fig. ^ it is clear that 
the longitudinal tension results in a restoring force in the z-direction when the plate is bent, i.e., when d'ujdx^ + 0. The 
tension deforms the plate, as indicated in Fig. The displacement results in a strain field y^^ - \{du^ldxa + dualdxp + 
du/dxa ■ duldxp), where Greek indices run over the x and y coordinate only. Using Eq. |lj the work done by applying the 
tension is calculated. The resulting stretching energy is: 



f/r = ^ XT ^"'^'^"'^ Axdy, where Tap = J" cr„/j dz = hcTa, 



(16) 



The equation of motion for the vertical deflection of the plate is obtained when the variation of the total potential energy 
U - Ub + Ut + Uf is considered {Uf - - Jj Fudxdy includes the eflrect of an external force per unit area F in the 
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Figure 7: Cross-sections of a rectangular beam with tliickness h and widtli w (a) and a cylindrical beam with radius r (b). 



z-direction) for an arbitrary variation in the displacement u ^ u + 6u. This yields the equation of motion for the plate 
PI: 

ph-^ + \dV' - —T,p—juix,y) = Fix,y), (17) 

where p is the mass density of the material|. The first term in Eq. |l^ is the inertial term and the last term is the 
external force acting on the plate. The term in brackets is the force resulting from the deformation of the plate. Here, 
D - Eh^ /12(1 -v^)is the so-called bending rigidity of the plate, the prefactor in Eq. |l5| that quantifies how much energy 
it costs to bend a unit area of the plate. The tension makes Eq. |l^ nonlinear as a displacement of the plate elongates it 
and thereby induces tension, i.e. Tap - Tap{u(x' ,y')'\(x,y). However, for small deformations the displacement-induced 
tension is small and will be overwhelmed by the bending rigidity or by tension induced by the clamping. Then the tension 
is independent of u and Eq. |l7| is linear. 

2.3.2. Beams 

Table |l] shows that many nanomechanical devices are doubly-clamped beams or cantilevers instead of plates i.e., they 
have a width that is much smaller then their length. For a beam, the normal components of the stress on the sides should 
also vanish, i.e., cr,,,, - cr„ = and the integration of Eq. [l^overy can be done directly. This implies that the tension 
Tap is only in the x-direction. Moreover, the material displaces in the x-direction so that T = T^.x is independent of x. 
Combining all of this yields the equation of motion in the "thin-beam" approximation: the Euler-BernouUi equation with 
tension: 

df- dx'^ dx?- 

where the crosssection area A equals wh for a rectangular beam and for a cylindrical beam A - np- (see Fig. The 
structure of Eq. |l8| is similar to that of the the equation of motion for a plate, Eq. |l7|. The first term is the acceleration 
and the term on the right-hand side is the external force. The restoring force due to the bending stiffness depends on 
the fourth order derivative of the displacement w.r.t. x. For the tension, this is a second-order derivative. The bending 
rigidity^ D - Eh?wl\2(\ - v^) can be written as the product of the Young's modulus and the second moment of inertia 
D — EI /{ I - v^), where the small correction (1 - v^) is often omitted JT^. For a rectangular beam the second moment of 
inertia is / = /i^w/12; for a (solid) cylinder with radius r it is / = [jlj]. The bendin g rigi dity depends strongly on 



the dimensions of the device. It can be as small as 10 Nm for a carbon nanotube (Sec. |2.4.4| ) or as large as 10 Nm 
for a millimetre-sized miiTor fbSn. 



2.4. Examples 

2.4.1. Cantilevers and doubly-clamped beams 

A cantilever is a beam that is clamped on one side (x = 0) and free on the other side {x - t). Because a cantilever 
is not fixed on the, say, right side, the tension is zero apart from more exotic cases where electrostatic forces [ (74| , |80| ] or 



surface tension [ 145, 146] acts on the free end. With T - Q, the eigenmodes u„ and (angular) eigenfrequencies a>„ satisfy: 



ujlpAun^D^. (19) 

The solution to this equation is a linear combination of the regular and hyperbolic sine and cosine functions (sin(A;jc), 
cos{kx), sinh(A;jic), and cosh(A;jic) resp.), where - uP'DlpA. Their coefficients are determined by the boundary conditions. 



^In principle, the first term should read phifu/di^ + ph^ d I dP' {d^ u I d>r + d'^u/dy^)/\2 as the material is also moving in the x and y direction. These 
corrections are, however, negligible when h/£ <K 1. 

'Note that the units of the bending rigidity, tension and external force are different from the case of a plate due to the integration over the y-coordinate. 
D is given in Nm^ instead of Nm(= J) and T is now in N instead of N/m. The external force F is given per unit length instead of per unit area. From 
the context it should be clear what the meaning of the different symbols is. 
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Table 4: Normalized eigenfrequencies and average mode deflection p„ for the first 5 flexural eigenmodes for cantilevers and doubly-clamped beams. 
fi„ = { ^„(x) dx indicates the displacement of the mode averaged along the length of the resonator per unit deflection. This number is important 
to calculate the detection efficiency for detectors that couple over the whole length of the beam (pCq). a„ and /?„ are the solutions of Eq. Ell and |l8| 
respectively which determine the eigenfrequencies. For large n the solutions approach a„ — > (n + 1 /2)n and j8„ — ► (n + 3/2)ff. 



Mode 




Cantilever 






Beam 




n 


a„ 


./«/./() 


Mil 


A, 


./«//() 







1.875 


1.000 


0.783 


4.730 


1.000 


0.831 


1 


4.694 


6.267 


0.434 


7.853 


2.757 





2 


7.855 


17.547 


0.254 


10.996 


5.404 


0.364 


3 


10.996 


34.386 


0.182 


14.137 


8.933 





4 


14.137 


56.843 


0.141 


17.279 


13.344 


0.231 



At the fixed end, the displacement is zero and the beam is horizontal: m(0) = and m'(0) = 0, where ' denotes differen- 
tiation with respect to x. At the free end, the force in the z-direction and the torque vanish, so u"({) - and u"'{t) - 0. 
There are thus four boundary conditions and four unknown coefficients. This system always has a trivial solution u„ - 
where all four coefficients are zero. There are, however, certain values k - k„ where one of the four boundary conditions 
is automatically satisfied. These values correspond to the eigenmodes of the flexural resonator. Using the other three 
boundary conditions, three of the coefficients are expressed in the fourth one (which we call ca). For a cantilever we 
define a„ = kn( and the n-th eigenmode is: 



X X 

Unix) - C4 1 sin(a„--) - sinh(ff„--) ■ 



sin(a„) + sinh(Q'„) 
cos(a„) + cosh(Q'„) 



X X 

cos(a„-) - cosh(a„-) 



(20) 



Note, that if u„(x) is an eigenmode of the cantilever, then c ■ u„{x) is one as well, for every value of c. In the following part 
we will therefore use the normalized eigenfunctions which satisfy ^ ' ^^dx - 1. The eigenfrequencies are given 
by: 

cos{a„)cosh{a„) +1=0, w„ = 2nf„a\r^ yjOlpA. (21) 

This equation can only be solved numerically. The first few solutions a„ are indicated in Table ^ and Fig. ^ shows the 
corresponding mode shapes. Unlike in a string under tension, a type of resonator that will be discussed in the next section, 
which has /„//o = (n+l), the eigenmodes do not have a harmonic spectrum. Also, the average displacement of the mode, 
defined as ju„ = ^ ^„ dx, decreases with increasing mode number n. This is because a part of the cantilever is moving 
upwards and an other part is moving downwards. 

The analysis for the flexural eigenmodes of a doubly-clamped beam, which does not have a free end but which is 
clamped at both x = and x = t, closely follows that for a cantilever. The difference is that the boundary conditions at 
X = { are now u{l) = and u'if) = 0. With /?„ = k„{ this yields: 



X X 

Unix) = C4 1 sin(y8„-) - sinh(;0„-) ■ 



sin(je„) + sinh(y3n) 
cos(j6„) + cosh(yS„) 



X X 

cos(j6„-) - cosh(y8„-) 



(22) 



and 



cos(j8„)cosh08„)- 1 =0, 



2nf„ 



(23) 

Table ^ shows that yS„ > a„, so that the eigenfrequencies of a clamped-clamped beam are higher than that of a cantilever 
with the same dimensions. This is because the additional clamping makes it stiffer. Moreover, due to symmetry /i„ 
vanishes for the odd modes of a beam. Figure |^ shows the first three mode shapes. 

Equations |1] and ^ show that when the resonator is made shorter while its transverse dimensions are kept the same, 
the eigenfrequencies increase due to the C^'^ term in the characteristic frequency scale Uch = '<jD/pA [147 148 1. On 
the other hand, the crosssectional area and the bending rigidity depend on the size of the resonator The exact scaling is 
a>ch h/fi and ojch rjfi for a rectangular and cylindrical resonator respectively . When making all dimensions of a 
resonator smaller [|8]|, the resonance frequencies thus increase inversely proportionally with the size. This scaling is the 
reason that nanomechanical resonators can have very high resonance frequencies of more than 1 GHz 



2.4.2. String resonators 

So far we have not considered the effect of tension in the examples of flexural resonators. In certain types of materials 
there can be so much stress that the bendingrigidity hardly contributes to the restoring force. The resonator is then a 
string under tension instead of a beam. Eq. |l8| shows that this is the case when T » Djl^, or when this is rewritten in 
the (longitudinal) strain y, y » (h/{)^/l2. Tension is thus more important in resonators with a large aspect ratio {/h. 
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Figure 8: Schematic of a singly-clamped cantilever (a) and a doublv-clamped beam (c). (b) and (d) show the shape of the first three (n = 0, 1, 2) flexural 
modes of these flexural resonators. This shape is calculated using Eqs. and p2| respectively. 



The eigenfrequencies of a string under tension are /„ = -\jT/pA x (n + l)/2£ and the corresponding modeshapes are 
^„(x) - y/2sm(7Tnx/(). 

Tension in the resonator arises when materials with different thermal expansion coefficients or diff'erent lattice con- 
stants are used. An example of the former is given by Regal et al. [|50|], where a 50;um long aluminum beam resonator 
is annealed at 150 to 350 °C, thereby increasing the resonance frequency from 237 kHz to 2.3 MHz. Stress can also be 
induced by growing heterostructures with different lattice constants. If the layer with the resonator has a smaller lattice 
constant that the layer underneath it, then the resonator is strained. By engineering the heterostucture, different amounts 
of tension can be induced in the resonator [149, 150|]. Another way of inducing tension is by placing the resonator on a 



flexible substrate that can be bent. As the top part of the substrate is elongated, the resonator becomes strained. With this 
technique, the resonance frequency of a resonator has been tuned by more than a factor 5 [ 1 5 1 1 . The most commonly used 
high-stress mater ial is silicon nitride. Under the appropriate growing conditions thin films wi th st r esse s of ~ 1 GPa can 
be obtained [ 152 1 and resonators made using this material have extremely high quality factors jl53[ 1 54 ] . Recently, high- 
stress SiN has therefore become a very popular material for nanomechanical resonators [^, 57, 69t 151, 7^ . Another 
important type of string resonators are suspended carbon nanotubes. Since in this system the displacement can r eadily be 
of the same order as the radius of the tube, nonlinear effects are important. These systems will be treated in Sec. 2.4.4 . 

Although high Q-values seem to be a general observation for strained beams, the mechanism behind the increase 
is not completely clear. Recall that the Q-factor is proportional to the ratio between the energy stored in the oscillator 
and the energy dissipated per cycle. Supported by experiments, the authors in Ref demonstrate that stress does 
not substantially change the dissipation rate in strained beam mechanical oscillators but rather significantly increases the 
elastic energy stored in the resonator They argue that the microscopic origin of the damping lies in localized defect states 
in the material. In the abovementioned experiments by Regal et al. the quality factor of an aluminum beam resonator was 
increased by a factor 50 by annealing it [|o|. Although the induced tension also increased the resonance frequency by a 
factor 10, the quality factor increased more. This means that in this case the annealing did not just increase the energy in 
the resonator, but it actually reduced the damping rate (y^ - cor/Q), which is the ratio of the resonance frequency ojr and 
the quality factor Also in this system the physical mechanism behind this remarkable increase is still unclear 



2.4.3. Buckled beams 

In the case of a string resonator the tension is tensile, i.e., the resonator is elongated. It is also possible that the 
resonator is compressed. In that case, the tension is negative. When the tension exceeds a critical value, it is energetically 
favorable for the beam to have a non-zero flexural displacement; the beam buckles. Buckled beams form an important 
class of (nano)mechanical resonators, and can be formed when the resonator consists of layers of different materials 



[155, 156, 19p, or using thermal expansion | 

The starting point for the analysis of buckled beams is the Euler-Bernoulli beam equation (Eq. |18|) with no external 
force, i.e., F - 0. The boundary conditions are the same as for a doubly-clamped beam: m(0) = u({) - and m'(0) = 
u'({) - [|l9, 157]. The negative tension has two contributions: first the beam is compressed by the fact that it is clamped. 
Figure ^ shows that the length of the free beam (A)) and that of the clamped beam ({) are different when the clamping 
point exert a longitudinal tension on the beam. This is the so-called residual tension, Tq - EAjo, with yo - ({ - A))/^o- 
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Figure 9: (a) The distance between the clamping points ( of a doubly clamped beam differs from the length of the free beam fg. If the induced tension 
is large enough, it is energetically favorable for the beam to displace, releasing strain energy at the cost of bending energy. This is called buckling. Due 
to the displacement u(x) the length of the beam is extended to L(t). (b) The total potential energy Ur + Ub for a beam below (dotted) and above (solid) 
the buckling threshold. In the former case the potential energy only has a single minimum at Umax = 0, whereai^the latter has two minima at non-zero 
deflection, (c) An example of the calculated eigenfrequencies of a 50,um-long beam with D = 3.20 ■ 10"'^ J ||l9||. The mode shapes at the position of 
the dots are shown. 



The second, positive, contribution comes from the stretching of the beam when it flexures. The resuhing length of the 
flexed beam is denoted by L = ^ {1 + {du/dxf'}^^^ dx. Combining both effects gives for smaU (du/dx <c 1) deflections: 



EA r^duf 

In the absence of driving, both the displacement and the tension are time-independent and the static deflection Mdc satisfies: 

d^UAr d'^UAc EA I dudc\^ 

- = = ^ 27 Jo (17) ^''^ 

This only has a non-trivial solution Mdc ^ when Tdc - n^Tc, where Tc = -An^D/L^ is the critical tension at which the 
beam buckles and n is an integer The solution is then Mdc(j^) = Mmax[l - cos(27rnx/^)]/2. When an initially unstrained 
beam is compressed slightly, work is done and the energy stored in Ut increases. The potential energy has a single 
minimum around M^ax = in this case; see the dotted line in Fig. ||b. When the compressive residual tension Tq is made 
more negative than it becomes energetically favorable for the beam to convert a part of Ut into the bending energy 
Ub- The potential energy now has two minima at non-zero static displacements, as illustrated by the solid line in Fig. 
The beam buckles to a displacement that keeps the tension exactly at Tc (for n-1). The value of the displacement depends 
on the residual tension that caused it: Mmax = +2{/7t ■ ([Tc - T()]/EAy^^ for To < Tc < 0. Note, that in the absence of a 
static force, the beam does not have a preferential direction of buckling. 

To find the eigenmodes of the buckled beam, we do not only focus on the static deflection, but we also include the 
dynamics of the displacement. When driving the modes of the beam, the total deflection u is the sum of the static (Mdc) 
and an oscillating part (Mac)- The time-dependent displacement satisfies: 

dul^ d'^Uac 5^Mac O^UAr 



with: 

dUdc 5Ma, 



_EA r' 

" ^ Jo ■ 



dx dx 



dx. (27) 



Note, that both sides of Eq. ^ are linear in Mac and that the static displacement Mdc thus acts as an effective spring constant 
for ac motio n, as indicated by the r.h.s. of Eq. |6[ The eigenfrequencies of buckled beams were calculated by Nayfeh 
et al. [157|. As an example the eigenfrequencies and modes for the beam used in Ref. [|l9[ are shown in Fig. 



At zero buckling the frequency of the fundamental mode is zero as the potential energy is quartic in the displacement: 
the quadratic terms in Ub and Uj cancel each other exactly. The frequency of the fundamental mode increases with 
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increasing buckling due to the contribution of Tac- The first higher mode has an eigenfrequency ui/ln - 1 .44 MHz and is 
independent of Umax, as the mode is anti-symmetric around the node, giving T-^c = 0. When u^nx is increased to 0.92 /urn, 
the two lowest modes cross and the fundamental mode is higher in frequency than the first odd mode.Q Thus, when the 
length of the beam, bending rigidity and the buckling displacement are known, the eigenfrequencies and modeshapes of a 
buckled beam can be calculated. 



2.4.4. Nanobeams 
In Sec 



2.4.2 we showed that tension can overwhelm the influence of bending rigidity in nanomechanical devices 
with high aspect ratios. For these thin wires the displacement can be of the order of the resonator thickness and in this 
case nonlinear effects connected with the deflection-induced tension become important. Suspended carbon nanotubes 
resonators are prototypical examples where the induced tension can be so large that the tubes are tuned from bending 
(beam-like) to tension dominated (string-like) [p^ 



158] 



The induced tension is a key feature of thin resonators as it can be used to electrically tune their frequency over a large 
range. In this Subsection we derive the equations for the frequency tuning and give some typical numbers for single-wall 
CNTs. We will start with considering the device geometry of a suspended CNT resonator and derive some of the basic 
equations describing the electrostatics of the problem. The analysis is, however, also applicable to other thin string-like 
resonators such as multi-wall CNTs, suspended graphene nanoribbons, or long suspended nanowires made from inorganic 
materials (see Table ||). 

Fig. sh ows an atomic force microscope image of a suspended single-wall carbon nanotube in a three terminal 
geometry [159|. The nanotube is connected to source and drain electrodes, enabling electrical transport measurements. 
The tube is suspended above a gate electrode at a distance hg, which cannot only be used to change the electrostatic 
potential on the tube, but also to drive the resonator [160| and to induce tension. An electrostatic force can be applied to 
the tube by applying a voltage Vg between a gate electrode (Fig. |l0|b) and the nanotube. The potential energy||^ depends 
on the capacitance between the tube and the gate Cg, and is Uf = -CgVg /2. The potential energy depends on the distance 
between the gate and the tube hg - u{x) via the gate capacitance, which means that there is a force acting acting on the 

nanotube. The electrostatic potential energy is written as: Uf--J^^ Cg(x)Vg/2dx, where Cg(x) is the capacitance per 

„{ 

unit length, and the potential energy equals by definition Uf = - Jq Fu dx. The electrostatic force per unit length is thus 

Fit) = {dcglduVj. 

To calculate the displacement dependence of the capacitance, we first consider the spatial profile of the electrostatic 
potential. Under the assumption that the screening effect of the source and drain electrodes is negligible, the tube is 
viewed as an infinitely long grounded solid cylinder, suspended above a conducting plate at an electrostatic potential 
(p{z = 0) = Vg. The potential profile for m = is given by [ 162 1: 



C^(y,z)^Vg-Vg 



1 



arccosh(/i„/r) 



In 



z + 



+ r 



z- ^hl - r2 



(28) 



The field lines associated with this potential are shown in Fig. |10|b. The deflection of the nanotube is included by replacing 
/ig with hg - u. After dividing the induced charge by the gate voltage, the capacitance per unit length Cg{x) is obtained]^ 



Cgix) = 



2n£o 



Inen 



arccosh([/!g-M(x)]/r) arccosh(/!g/r) - arccosh^(/!g/r) 



m(x). 



(29) 



The approximation in Eq. g9| is allowed because the displacement u is typically much smaller than hg. This in contrast to 
top-down fabricated devices, where higher order terms can be important and electrostatic softening of the spring constant 
might occur [163, ^7|]. This softening is thus due to nonlinear capacitance terms, which should be contrasted to the 
nonlinearities due to the displacement-induced tension that we will discuss further on in this Section. For thin resonators 
under high built-in tension, however, the nonlinear capacitance terms can become relevant again. This has been nicely 
illustrated in graphene resonators, which show a small downward shift in the frequency around zero-gate voltage [104]. 
We will not consider these effects in the following, but we finally note that the electrical softening is closely related to the 



optical spring that will be discussed in Sec. 3.5.3 



The modes are classified by their shape and not by the ordering of eigenfrequencies. The fundamental mode is the mode without a node. 
"This is the potential energy for the tube, in contra st to the energy stored in the capacitor: +Cg Vg/2. The difference in the sign is because the voltage 
source performs work when the capacitance changes 1 161 1, whi ch sh ould also be taken into account. 

'^This expression might appear different from those in Refs [ 135 1 and 1 158|, but note that arccosh(A') = ln(jc + Vjc^ - 1) » ln(2A:) for A" 2> 1. 
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Figure 10: (a) an AFM image of a suspended nanotube device connected to the source and drain electrode. The tube is suspended above the trench only. 
The suspended part of this device is C = 1 .25 /im long and the radius of the tube is r = 1 .4 nm. (b) Field lines of the electrostatic potential induced by 
the gate electrode. When the distance between the tube and the gate hg changes, the gate capacitance Cg changes. 



The electrostatic force per unit length is now: 

1 dcg , neoVlit) 
Fit) ^ --^Vl(t) ^ ——= ^ . (30) 

^ -J/ig - arccosh^^/ig/r) 

Often, the gate voltage consists of two parts: a static part V^^ and a time-dependent part cos(wf) to drive the nano- 
tube at frequency / = cjlln. The experimental condition V^*^ <K Vt'^ ensures that terms proportional to (yi'^)^ are 
negligible. The force is then the sum of a static and driving contribution: F - F^^^ + F^c cos(a>t), with F^^^ - 

neo/ihl - r2)'/2arccosh2(/!g/r) ■ (Vf )2 and f ac(f) = ^eo/(/?g - r^y^^aTccosh^[hg/r) ■ 2yf y|'=(f). 

The bending mode vibrations are described by the Euler-Bernoulli beam equation with tension included, i.e., Eq. [T^ . 
When the amplitude of the oscillation Mac is small compared to the larger of the tube's radius and the static displacement, 
terms proportional to ul^ are negligible and the tube is in the linear regime]^ Similar to the analysis presented for the 
buckled beams, the equation of motion (Cf. Eq. 18) is separated: 

^^-"--^ =^dc, (31) 

-CO pAu^c + icoyui^c + D-^-^ - Td^-^-^ -T^^-^ = /^ac- (32) 

The first equation describes the static displacement of the tube that is induced by the dc gate voltage. This equation is 
independent of Mac- On the other hand, the ac displacement, given by Eq. ^ depends on the static displacement and also 
on the static tension Tdc- Similar to the case of the buckled beam, the static tension has two contributions, as indicated 
by Eq. |5| The first one is the residual tension due to the clamping as the length of the suspended part is not necessarily 
equal to the length when it would not be clamped. For example, a nanotube resonator could be strained during the growth 
process or lay slightly curved on the substrate before suspending it. The second contribution is the displacement-induced 
tension: The gate electrode pulls the resonator towards it, thereby elongating it. Moreover, the oscillator experiences a 
time-dependent variation in its length and Tac is the part of the tension that is linear in Mac- Both effects ar e included in 
Eq. As the tension contains the static displacement, it has to be solved self-consistently [158 135 159 1 with Eq. 31 
to find the static displacement. The resulting static tension, ac tensi on a nd dc displacement are then inserted into Eq. 32 

to find the eigenfrequencies <y„ and the response function Mac(j^, <y) [164|. 

To analyze the system of Eqs. ^ |T and ^ it is useful to take a closer look at their scaling behavior [135|. In 
this section, we use the convention that primed variables indicate scaled (dimensionless) variables. An obvious way to 
normalize the coordinate x is to divide it by the tube length: x' - x/{, so that the equation for the static displacement Eq. 
pT] becomes: 

g'^Mdc £^Tdc d^Udc _ {"^F^^. _ ^ 

dx'^ D 5x'2 " D - ^ ^ 

On the right hand side, a natural length scale for the static displacement, l^c, appears. However, scaling the displacement 
with lie is not handy because Zjc equals zero at zero gate voltage. Therefore Mdc (and Mac) are scaled by the radius of the 
tube: M^^ = Mdc/'"- Moreover, the tension has become dimensionless, resulting in an equation for the static displacement 
where the number of parameters has been reduced from 5 to 2: 

" (34) 



5x'4 "'^ dx- 



'^This refers to the dynamical behavior. The static displacement is actually nonlinear when the tube is in the strong bending limit 1 158]. 
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Figure 11: (a) static displacement profiles for T'^^ = (top), and T'^^ = 10^ (bottom). (b,c) The calculated static displacement at the center of the 
nanotube (b) and the corresponding tension (c) for various value of the residual tension T^. The limits for the tension for small and large static forces 
are indicated. 




Figur&^12: Gate-voltage tuning of the eigenfrequencies of the first (a) and second (b) flexural mode of a 1 /im-long suspended carbon nanotube (see 
Table H). The residual tensions are To = -91, 0, + 230 pN for the solid, dashed and dotted lines respectively, which corresponds to the dimensionless 
values T' = -39.4, 0, + 100. 



where 



T-dc - — = To + — (^) , r„ = — , and = l,Jr. (35) 



The definition of T'^^ with the dependence shows that tension becomes more and more important when the length of the 
device increases. Figure pT^ shows the dc displacement profiles for different values of the static tension. In the case where 
the bending rigidity dominates (top panel) the profile is rounded at the edge, whereas for high tension (lower panel) the 
profile is much sharper. The tension and center deflection are calculated by self-consistently solving Eq. ^ with Eg . ^ 
and are plotted in Fig. [lip and c. Two different slopes can be distinguished in the double-logarithmic plot of Fig. |l l| b. 
These correspond to the weak and strong bending regime [|158p, where Wdc is proportional to F^j^ and F^^ repectively. 

The two regimes cross at T^* = 6y/T0 ~ 50.2, Z^* = 36 ■ 70-'''* ~ 871. The gate voltage at which 1'^^ - Z^* is called 
the cross-over voltage, V*. AFM measurements of th e gate-induced displacement of multi- walled carbon nanotubes have 
confirmed this scaling behavior experimentally [135|. 

To calculate the gate-tuning of the resonance frequency fuiVg'^), the scaling analysis is also applied to the equation for 
the ac displacement, Eq. ^% One immediately finds the length scale /^c - i'^F^.cID for the ac force and 7"^^ = T^ci^/EI. 



As in Sec. 2.4.1 , Wc/, = {D/pAy^^/i^ is again the characteristic frequency scale for the bending mode vibrations. Next, 
one has to solve Eqs. ^ and |3j to obtain the static tension and dc displacement. Then, the boundary conditions are 
imposed to find the resonance frequencies. Figure |lja shows the calculated eigenfrequencies, plotted against the static 
pulling force for dififerent residual tensions. The higher the residual tension is, the higher the resonance frequencies are at 
a given V^^. The value - -39.4 ^ indicates that the thin resonator is close to buckling; this is visible by the nearly 
vanishing resonance frequency of the first mode at low V^'^. When the static force is increased, all resonance frequencies 
increase and the diff'erences between the curves due to the diff'erent residual tensions become smaller. 

To relate the dimensionless quantities in this Subsection to physical ones, the system dimensions are needed. Table 
H shows the estimated sizes and the calculated values of several parameters for two typical nanotube devices, one with a 
length ^ = 1 jjm and one with a length { - 200 nm. Note, that with our definition of u„, the mass and spring constant 
appearing in the zero-point motion mq (Eq. |Tj) and the equipartition theorem (Eq . ^4| ) are equal to the total mass m and 



kg - mw„ respectively. There is no need to introduce an effective mass; see Sec. 3.1. The table indicates that for a long 



device tension is important and the diff'erence between /« and /o is large. For a short resonator, the mechanical properties 
are mainly determined by the bending rigidity and /« x /g. The higher value of V* indicates that a larger gate voltage has 
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Table 5: Data for two suspended cai'bon nanotubes with different lengths C. The values of the parameters are calculated with a nanotube radius 
r = 1.5 nm, and distance to the gate electrode ftg = 500nm. The resonance frequency fit, tension Tjc and static displacement u^^.{f/2) are evaluated 
for a gate voltage Vg*^ = 4 V and zero residual tension, m, fo and D are the mass, resonance freque ncy of the nan otube without tension and the bending 
rigidity respectively. Furthermore, ijc and kn = mw^ are the static and dynamicspring constants |165, 
uq the zero-point motion (see Sec. |l|for its definition and the discussion in Sec. hi). 
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to be applied to tune this resonator to the strong bending regime. 

As a final note, we consider the difference between a nanobeam resonator wit h a h igh frequency and one with a low 
frequency that is tuned using a gate vol tage t o the same frequency as the former [159]. The current associated with the 
motion in mixing experiments (see Sec. 4.3.3 ) ^ Q |2[ |8j] is proportional to the length-averaged displacement 

amplitude yu„M„(w); it is thus proportional to the length scale 4c that determines the amplitude. As a consequence, the 
mechanical signal drops rapidly with decreasing length, making the measurement of single-walled carbon nanotubes with 
/o > 1 GHz (corresponding to a device with { < 0.2 /im and r = 1 nm) challenging as the signal is about lOOx smaller 
compared to a fo = 100 MHz device with £ ^ 0.6 fj.m. The latter tube can also operate at 1 GHz frequency by tuning it 
with a gate-induced tension of T'^^ - 5 ■ 10^. In this case the signal decreases too, but only by a factor of 10. A tension of 
T'^^ = 5 ■ lO^ corresponds to a strain of about 0.2%, which is larger than the val ues in the abov ementioned experiments, 
but is still smaller than the strain at which single-walled carbon nanotubes break [168, 169| , 170 1. 



2.4.5. Nanodrums 

So far, we have only considered one-dimensional resonators such as beams and strings. An example of a two dimen- 
sional resonator is a graphene nanodrum. This device consists of a hole that is etched in a substrate and that is covered by 
a (few-layer) graphene flake. Surprisingly, an one-atom layer can be suspended and holes with a diameter of 100 /im have 
been reported [171]. A much smaller version of these devices (a nanodrum) has been used to study the (im)permeability 
of graphene to gases [100| and to measure the bending rigidity of and tension in the flake using an atomic force micro- 
scope [107, 172]. In the latter experiments, an atomic force microscope tip is used to apply a force F^ip to the flake as 
illustrated in Fig. pj|a. This versatile technique has also been applied to other geometries [173, 174 1 and nanomaterials 
[175, 176, 177]. The point (ro, ^o) where the force is applied can be varied and the resulting deflection of the nanodrum 
is measured. The restoring force that opposes the applied force has several contributions. First of all, there is the bending 
rigidity of the flake D, and secondly, tension may be present in the flake. 

Since graphite is highly anisotropic, the analysis of the bending rigidity of an isotropic material in Sec. has to be 
generalized]^ Using the compliance tensor in Eq. [lljthe rigidity for bending along the sheets (see Fig. ||b) is calculated: 



D = £q/!712(1 



(36) 



which only contains the in-plane elastic constants. However, when the number of graphene layers becomes small, correc- 
tions to Eq. ^ have to be made. Consider the situation in Fig. ||b where a few-layer graphene sheet is bent with a radius 



'"iThe analysis is in principle also valid for other layers two-dimensional membranes; for a thin isotropic membrane Eq. | 
following it are still valid if the appropriate bending rigidity D is taken. 



and the derivation 
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of curvature Re- In the continuum case, the bending energy is given by Ub/{W = jDjR^. - j E^{z/Rc)^ dz, when 
taking v = for simpUcity. In the case of a small number of layers A^, the continuum approximation in the z-direction is 
no longer valid and the stress is located only at the position of the sheets z, - c(i - [N + l]/2), where c = 0.335 nm is the 



inter-layer spacing (see p. 12). The integral over z is replaced by a sum and the bending rigidity becomes: 



Dm = 



Enh^ 



N-1 



12(1 -v2^) 



(37) 



where the thickness is set by the number of layers h = Nc, which reduced to Eq. ^ in the limit N ^ oo. According to Eq. 

the bending rigidity vanishes for a single layer. However, molecular dynamics simulatio ns have shown that a single 
layer of graphene still has a finite bending rigidity: of the order of one eV [117, 178[ , 179| (compare this to the value 
for a double layer (A^ = 2) calculated with Eq. ^ D2 - 54 eV). The rigidity of a single layer comes from the fact that 
electrons in the delocalized 7r-orbitals, located below and above the sheet, repel each other when the sheet is bent [179]. 

The deflection of the nanodrum satisfies the equation for the deflection of a plate, Eq. 17, and when the force applied 
by the AFM tip is assumed to be located at a sing le point (xo,yo), one gets [|l38[ |l8^ [l8l] |l82[|: 



DV^ - - — Tap— I u{x,y; Xfi,yo) = F^^dix - xo)6{y - yo). 

OXa OXj} J 



(38) 



Here, the V-operator and the partial derivatives djdxi are working in the xy-plane only, as the z-dependence is absorbed 



r- 1 W? 

in the bending rigidity (see Sec. [2.3| ) and where the tension Tap - Jq crap dz. This equation is difficult to solve in its 
most general form, but fortunately some simplifications can be made: The tension tensor can have both normal and 
shear components. It is, however, always possible to find two orthogonal directions where the shear components are zero 
[ |136( ]. When we assume that the tension is uniform then these directions are independent of position, so without loss of 
generality the x and y-axis are taken along the principle directions of the tension. When the difference in tension in the x 
and y direction, AT - (r^ - Tyy)/2, is smaU, first the solution for isotropic tension (Tap ~ T6ap) is obtained and then the 
correction due to the finite AT" can be calculated 1 165 1. Here, we will focus on the situation where AT - 0. For a circular 



membrane, it is convenient to use polar coordinates and the equation for the displacement reads|^ 

F ■ 

(dV^ - TV^)u(r, 6- ro, 0o) = -^6(r - ro)6(e - 0q). 



(39) 



The solution is written as: 



u(r, 6; ro, 9o) = ^ Rm(r, ro) cos(ot6' - m6o). 



(40) 



m=() 



Inserting this into Eq. M yields for the radial coefficients: 



Roin ro) = Aok(Ar/R) + BoKo(Ar/R) + 
Coln(r//?) + A)+<'(r;ro), 
Rm(r, ro) = A„I„(Ar/R) + B,nK,„(Ar/R) + 

CUr/R)-'" + D„,{rlRr + R^^\n ro) (m > 0), 



(41) 
(42) 



where /,„ and K„ are the Bessel functions of the first and second kind respectively, and A - ^jTR■^|D is a dimensionless 
parameter that indicates the importance of the tension in comparison with the bending rigidity of the flake. 

The flake with radius R is clamped at the edge of the circular hole so the boundary conditions are u{R) - and 
dM/dr|r=R = 0. Furthermore, the deflection at the center is finite and smooth (i.e., dM/dr|r=o = 0). The set of coefficients 
{Am,Bm,Cm,D,„) Can be calculated analytically. Figure [l3|b shows the deflection profiles calculated where the force is 
applied at different distances ro from the center. The deflection of the flake is clearly reduced when the AFM tip is 
moved away from center of the nanodrum. This indicates that its local compliance fcy'(ro,6'o) - du{ro,6o; rQ,6o)/dFtip 
decreases. As the tension is assumed to be isotropic, k^^ is independent of 9q and its radial profile, shown in Fig. [T3| c 
for different values of the tension, contains all the information. In analogy with the displacement profile of a bending 
and tension-dominated carbon nanotube (Fig. [ll|a), the profile is rounded at the edge of the hole for vanishing tension 
(A - 0), whereas for large tension (A — > 00) the compliance profile becomes much shapei]^ at the edge and diverges at 



'*In the limit /i — > oo, tlie term in Eq. [39 vanishes and only a second order differential equation remains. Therefore, the boundary conditions 
du/dr|r=o,i; = are discarded. 
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Figure 13: (a) Schematic overview of the nanodrum. A few-layer graphene flake is suspended over a circular hole with radius R. A force, is applied at the 
point (ro, So) using an AFM tip. This results in a deflection of the nanodrum (b) Colormaps of the calculated deflection profile (Eq. B% of a nanodrum 
with vanishing tension. The force is applied at the location of the cross and the color scale is identical in all four panels: white corresponds to a large 
deflection and dark gray to no deflection, (c) The calculated radial compliance profile for different values of the tension, with A- = TR- ID. 



the center for a point force. In pract ice the tip has a finite radius of curvature which prevents that the spring constant 
of the flake A;/(ro = 0) vanishes 1 172 |. By co mparing the experimentafly measured profile with the calculated ones, the 
values for D and T can be determined | 107 |. Knowing the tension (and in principle also the bending rigidity, but this 
is negligible for thin flakes) an estimate of the resonance frequency can be made. The fundamental eigenfrequency of a 
circular drumhead is /o = ^[Tjp2dl2nR ■ Vg, where vq ~ 2.4048 is the first zero of the zero-order Bessel function of th e 
first kind. Measurements have indicated a typical value ofT~ 0.3 N/m for exfoliated single layer graphene flakes [172|. 
With this value and the value for the two-dimensional mass density of single-layer graphene, p2d - 6.8 ■ 10"'' kg/m^ (see 
page we obtain a resonance frequency of 100 MHz for a nanodrum with a diameter of 2/? = 5 /vm. By decreasing the 
size of the hole, the frequency can be increased. For a hole with a diameter of 500 nm, the frequency akeady exceeds 1 
GHz. The latter resonator has a mass of only m - 5 ■ 10 '^kg and a very large zero-point motion of mq = 0.4 pm. This 
makes graphene nanodrums ideal devices for QEMS. 
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3. Backaction and cooling 



In this Section we consider the influence of the detector, the so-called backaction and we discuss ways to cool the 
resonator to the ground state. To describe these eff'ects in detail concepts such as thermal noise, Brownian motion, and the 
eff'ective resonator temperature are introduced. Next we will discuss detector backaction and show that it plays an impor- 
tant role, especially when continuously probing the resonator properties. Backaction is therefore often an unwanted, but 
unavoidable element in measurements on mechanical systems. It determines how precise the position can be monitored; 
ultimately, quantum mechanics poses a limit on continuous linear detection, the standard quantum limit. We will derive 
this limit in different ways. The standard quantum limit can be circumvented by employing different detection schemes, 
such as square-law detection and backaction evading measurements. These will however not be discussed in depth in this 
Review. 

Backaction does not just impose limits, it can also work to one's advantage: It can squeeze the resonator motion 



[ |183| , |184[ |185| , |186[ |187| |, couple and synchronize multiple resonators [|188|, |189p, and backaction can cool the resonator as 
we will show. Besides this self-cooling by backaction, other active cooling schemes have been developed. In this Section 
we will discuss two different cooling schemes, namely active feedback cooling and sideband cooling. It is interesting to 
note that, at present, cooling has only been performed on top-down devices because the experiments to observe thermal 
motion of bottom-up devices are more challenging. 

3.1. From modes to harmonic oscillators 

Before starting the discussion of backaction, first an important point has to be addressed. In the previous Section 
we have given the equations of motion for several nanomechanical systems. Solving these equations gives the frequency 
and the displacement profile of a particular mode. An important conclusion of describing small displacements around 
the equilibrium position of NEMS is that each mode can be viewed as a harmonic oscillator To show this, we start with 
expanding the displacement in the basis[^ formed by the eigenfunctions ^„ [191, 192, 159, 164|: 



M(x,f) = 2M*'"(f)^„W. (43) 

n 

Inserting this into the equation of motion and taking the inner product with ^„ yields the displacement of mode n. For 



example, for the nanobeams discussed in Sec. 2.4.4, Eq. B3 yields: 



(44) 



The left hand side shows that the frequency response of each mode is equal to the response function of a damped driven 
harmonic oscillator. The same conclusion is reached for the other examples: all their modes can be describes as harmonic 
oscillators. The mathematical reason behind this is that the equation of motion for small deformations of a mechanical 
system can be written in the form mefrM(r, f) = -yM(r, f) + Jl[u(r, t)] for some Hermitian operator X.. Its eigenfunctions 
are the mechanical modes and these form a complete orthogonal basis. After expanding the displacement in this basis, a 
set of uncoupled harmonic oscillators results. An important question that one should ask after the transformation from the 
spatial modes to the harmonic oscillators is: What is the effective mass njetf of the oscillator? This question might seem 
trivial at first, but the concept of the effective mass has given rise to much confusion in the past years. From the equation 
of motion it follows that the effective mass of mode n equals m^g^,, - Jy p^^ dV. The effective mass thus depends on the 
normalization of the basis functions, or equivalently on the definition of the mode displacement u„. In this Review we 
have adopted the convention that the basis functions ^„ are orthonormal, i.e. ^ dx = 1. In this caseQ the effective 
mass equals the total mass of the system (i.e., meff.„ - m). On the other hand, one can also use different normalizations, 
for example using the average displacement {t ^„ dx - 1) or using the maximum displacement (max^„ = 1). In these 
cases, the effective mass differs from the total mass and it depends on the exact mode profile. To illustrate the confusion 
this may create, consider the flexural modes of a beam resonator: With the latter two definitions, every mode has a 
different effective mass. Moreover, tuning a flexural resonator from the bending to stretching-dominated regime (Sec. 



Z.4.4 ) changes the effective mass of its modes. These complications are avoided by using orthonormal basis functions, 
wherer' 1. 

Although some groups study the interaction between different coupled modes in nano and micromechanical systems 



[ 193 , 194| , |195| , 164 , 124[ 196 |, most experimental and theoretical work focuses on a single mode only and the harmonic 



oscillator describes the dynamics of the entire mechanical system. In this case, we make no distinction between the 



'It can be shown that these functions form a basis for functions that satisfy the homogeneous boundary conditions [ 19C , 
'^For simplicity it is assumed that the mass density is constant. 
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resonator (i.e., the entire mechanical system) and the mode that is studied. However, we note again that the actual 
displacement of the resonator at a certain position x is given by m*"'^„(x) and not by m*"' itself. 

Before reviewing the properties of the classical and quantum harmonic oscillator, we stress, however, that knowledge 
about the displacement profile remains important when analyzing NEMS experiments. In particular, different detectors 
or driving forces may couple differently to the displacement profile and could detect therefore different modes. In most 
cases both the driving and detection mechanisms couple to the average displacement of the resonator In this case, anti- 
symmetric modes are not visible in nanomechanical experiments as these have a vanishing value of the length-averaged 



displacement yu„ (seeEq. M4j and Table W). An example is a nanotube resonator with frequency mixing readout (Sees. 2.4.4 



and 4.3.3 ) that is either coupled to a back gate or to a local gate; the former couples uniformly to the whole nanotube, so 
that the detected signal is proportional to and consequently only symmetric modes can be detected. A local gate may, 
depending on its position, couple to all modes. 

3.2. The harmonic oscillator 

The harmonic oscillator is probably the most extensively studied system in physics. Nearly everything that returns 
to its equilibrium position after being displaced can be described by a harmonic oscillator. Examples range from the 
suspension of a car, traffic-induced vibrations of a bridge, and the voltage in an electrical LC network, to light in an 
optical cavity. For large amplitudes, the oscillator can become nonlinear. We will not study that situation in this Report, 



but instead we refer the reader to the large body of literature on this subject, see e.g. Refs. 1 197 198]. We will now 



proceed with describing the classical and quantum harmonic oscillator in more detail and reviewing their basic properties. 

3.2.1. The classical harmonic oscillator 

In a harmonic oscillator, the potential energy depends quadratically on the displacement u from the equilibrium posi- 
tion: 

V{u) = ]^kRU^, (45) 

where Ur is the spring constant. The parabolic shape of the potential results in a force that is proportional to the displace- 
ment. When damping and a driving force Fit) are included, the equation of motion reads: 

mu — —Uru - mjRii + F(t), (46) 

for a resonator with mass m and damping rate jr. When the oscillator is displaced and released, it will oscillate at 
frequency ojr with a slowly decreasing amplitude due to the damping. The quality factor Q - (jJrIjr indicates how many 
times the resonator moves back and forth before its energy has decreased by a factor e. 

The harmonic oscillator responds linearly to an applied force; in other words, it is a linear system. Any linear system 



is characterized by its impulse response or Green's function |199]. For the harmonic oscillator, the impulse response 
huoit), is the solution to Eq. ^with F{t) - kR6(a>Rt): 

huoit) = sin(wRf)e ©(wrO, (47) 

where ©(f) is the Heaviside stepfunction. The impulse response function^ describes how the resonator reacts to a kick at 
time f = and is plotted in Fig. [l4| a. The initial displacement is zero, but the kick gives it the resonator a finite velocity 
at f = 0. The resonator then oscillates back and forth with a period 2nlu)R and these oscillations slowly die out due 
to the dissipation. With the impulse-response function the time evolution of the displacement for a force with arbitrary 
time-dependence F(t) can be obtained directly: 

u(t)^hHo{t)®F{t)lkR, (48) 

where the symbol ® denotes convolution. 

In many experiments the oscillator is driven with a periodic force Fit) - Fo cos(wf). After a short (~ y^') transient, 
the resonator oscillates with the same frequency as the driving force. This motion is not necessary in-phase with the 
driving signal. Both the amplitude and phase of the motion are quantified by the transfer function Hnoi'^) that is obtained 
by taking the Fourier transformation^ of the equation of motion, Eq. 46 



Hho(<^) -kR-^^ . (49) 



This is tlie Green's function for a liigli-Q resonator. For lower Q- values, the resonator oscillates at a slightly lower frequency ai'j^ = ior 

Vl- (1/22)2 

and the impulse response of an underdamped oscillator (i.e., one that has Q > 1/2) is: hnoit) = sm(ai'^t)exp(-ojRt/2Q) ■ [1 - (2Q)-^T^'^&(t). An 
overdamped resonator {Q < 1/2) returns to h = without any oscillatioQi and has a different impulse response. Throughout this Review it is assumed 
that Q 3> 1 so that n i^ « tj/; and hfjo is given by Eq. Wj. Note, that Eq. M is valid for all (positive) values of Q. 

-"By convention | |l9S| ], the Fourier transformation is defined as: X(a)) = ?^[x(t)] = jc(;) exp(-icjf) d; so that the inverse transformation is given 
by: xif) = r-^[Xico)] = ^ f^^ X{o)cxp(+iwt) do. 
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Figure 14: (a) The Green's function h[jo(t) and (b) frequency response Hho(M>) of a harmonic oscillator with Q ■ 
and the corresponding wave functions \l/„(u) (solid) of the harmonic oscillator for n = 0..7. 



10. (c) Eigenenergies E„ (dashed) 



The magnitude |////o(w)| and phase ZHho((^) are plotted in Fig. 143. When the driving frequency is far below the 
resonance frequency, the oscillator adiabatically follows the applied force: u(t) - F{t)/kR o Hho - 1. and so both 
\Hho\ and ZHho are small. The small motion is then almost in phase with the driving signal. The amplitude grows when 
sweeping the frequency towards the natural frequency. Exactly on resonance, the amplitude has its maximum \Hho\ - Q- 
The phase response shows that at the resonance frequency, the displacement lags the driving force by -n/l. When further 
increasing the driving frequency, the oscillator can no longer follow the driving force: the amplitude drops and the lag 
approaches -tt. The motion is then 180° out of phase with the applied force. The width of the resonance peak is related 
to the damping: the full width at half maximum of the resonance equals jr = cor/Q- 

3.2.2. The harmonic oscillator in quantum mechanics 

In quantum mechanics the harmonic oscillator is described by the Hamiltonian H - 12m + {moj^u^ [^, as the 
classical displacement coordinate u and momentum p = mil have to be replaced by the operators u and p - -ih ■ d/du. 
The displacement is described by a wave function (^(m) that satisfies the time-independent Schrodinger equation: 



2m ou'^ 2 



Eip. 



(50) 



{mcjRU — ip) I y/2mha>R and a = 
^), where h = a' a is the number 
h(L)R{n + 4), with eigenstates 



This equation is solved by introducing the creation and annihilation operators: a' 
{mcjRU + ip) I y]2mh<jjR respectively. The Hamiltonian then becomes H - ha)R{h 
operator that counts the number of phonons in the oscillator. The eigenenergies are £„ - ,^.^kv" ' 5 
|n). The corresponding wave functions i//n(u) are plotted in Fig. The lowest (n = 0) eigenstate has a non-zero 

energy Eq = jhcuiR, the so-called zero-point energy. Even when the oscillator relaxes completely, it still moves around the 
potential minimum at m = 0. The probability density of finding the resonator at position u, is given by \if/o(u)\^ when the 
resonator is in the ground state. The zero-point motion mq is the standard deviation of this probability density: 



1/2 



h 

^ 2mcoR 



(51) 



The zero-point motion is an important length scale that determines the quantum limit on continuous linear position mea- 
surement and is also related to the effective resonator temperature as the following Sections will show. 

3.3. Thermal and quantum noise 



In the previous Section, Sec. 3.2.2 it was shown that a resonator always moves because it contains at least the zero- 
point energy. In practise, except at the lowest temperatures, the zero-point motion is overwhelmed by thermal noise. 
Thermal noise is generated by the environment of the resonator. As an example, consider a resonator in air At room 
temperature, the air molecules have an average velocity of about 500 m/s. The molecules randomly hit the resonator and 
every collision gives the resonator a kick. These kicks occur independently of each other, so the resonator experiences 
a stochastic force F„(t). Other thermal noise sources are phonons in the substrate that couple to the resonator via the 
clamping points, fluctuating amounts of charge on nearby impurities and so on. The environment of the resonator is 
thus a source of random fluctuations on the oscillator. The force noise can be described by an autocorrelation function 
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Figure 15: (a) Simulated time-trace of the displacement (blue), amplitude (black) and phase (red) of a resonator that is driven by Gaussian white noise 
for Q = 50. The phase if is in radians and both the displacement u and ampHtude A are normalized by the root-mean-square displacement (h^)''^. The 
corresponding auto-con'elation functions and the displacement-noise power spectral density are shown in (b) and (c) respectively. The area under the 
curve in (c) equals the variance of the displacement. 



RF„F„i^^ -_^^fn{t')F„{t' + f)] (the symbol E denotes the expectation value) or by its power spectral density^ (PSD) 
S F„F„{'^) [200]. For white noise, the latter is independent of frequency and F„{t) has an infinite variance. For a given 
F„{t) the realized displacement is easily found using the Green's function, i.e., with Eq. ^ Figure shows a simulated 
time-trace of this so-called Brownian motion. The resonator oscillates back and forth with frequency lor while its phase 
and amplitude vary on a much longer timescale. The displacement can be written as uif) - A(f)cos[wsf -i- ip{t)\ (see 
Appendix ^ and the time-traces of the amplitude A and phase ip are plotted in Fig. p^ a as well. Figure |l5|b shows the 
calculated autocorrelation functions of the displacement, amplitude and phase. The displacement autocorrelation Ruu{t) 
displays oscillations with period ItiIior, whereas Raa and R^^ do not contain these rapid oscillations. All three functions 



fall off at timescales 
positive. 



QIlor. Note, that /?aa(0 does not go to zero for long times, because E[A] > as A(f) is always 



The displacement PSD is proportional to the force noise PSD and is given by [200|: 



5„„(w) = k,^\HHo{(^TS F„F.,- 



(52) 



When S f„f„ is white in the bandwidth of the resonator, which is typically assumed, S uu has the characteristic bell shape 
shown^in Fig. p3| :. Moreover, the PSD can be used to obtain the variance of the displacement: 



<«2) = 1 f 
2^ Jo 



4 kl 2 kl 



(53) 



The force noise PSD is related to the te mpera ture, and in equilibrium the resonator temperature equals the environmental 
temperature. The equipartition theorem [201 1 relates the variance of the displacement to the equilibrium temperature: 



1 2 1 2 1 



(54) 



The thermal energy k^T is distributed equally between the potential energy and the kinetic energy. By combining Eqs. ^ 
and B4I a relation between the force noise PSD and the properties of the resonator is found: 



S f^^f^X^) - ^k-gTmojul Q. 



(55) 



This so-called fluctuation-dissipation theorem [ [20 1| , [202| , [203| ] shows that on one hand the force noise PSD can directly be 
obtained from the resonator properties and temperature, without knowing its microscopic origin. On the other hand, the 
force noise determines the dissipation (i.e., quality factor) of the resonator. 

At equilibrium, the temperature of the resonator is proportional to the variance of its Brownian motion as indicated by 
Eq. However, out of equilibrium the force noise is no longer given by Eq. ^ and the resonator temperature can be 



^'The engineering convention for the single-sided power spectral density, S xx(t^) = Sxx{t^) + Sxxi—i^), is used. Here, Sxx(t^) = "FlRxx] is 
the double-sided PSD and Rxx is its autocorrelation function. The variance of X is given by (X^) = Rxx(P) = (2?!')"' ■ j_^S xx{'JJ)<i'JJ = (2?r)"' ■ 
Sxxit^) dw. 

^^This peak shape is often said to be Lorentzian, although formally that is not correct. The peak is proportional to [(w^ - oj^)^ + which 
can be approximated by ai'^-^ /[{ain - dSf + w^/Q^] for w » wr. The latter is indeed a Lorentzian, but the approximation is only valid for frequencies 
close to the resonance frequency of a high-Q resonator 
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Figure 16: The resonator temperature Tr plotted against the environmental temperature T . At low temperatures T the resonator temperature saturates 
at the zero-point energy: Tr = ^hcon/k^ and at high temperatures Tr = T (dashed line). The insets show the occupation probability P„ of the energy 
levels at /tBr/fiw« = 0.3, l/ln2 a 1.44 and 2.0. 



different from T. The effective resonator temperature Tr is defined as: 



_kR<y} kR r - doj kR r 2-f , .... 

Tr = — = — 5„„(w)— = — \Hr\ , (56) 

ks Kb Jo 2n ks Jo /tt 



which yields Tr = T in equilibrium. When the force noise is larger than that of Eq. £^ the effective resonator temperature 
is higher than the environmental temperature. When (u^) is smaller than its equilibrium value, Tr < T. Equation q6| shows 
that the resonator temperature can be obtained from the experimental displacement noise PSD: The resonator temperature 
is proportional to the area under the curve (Fig. p^). Note that the suggestive notation duj/ln - df is used in Eq. ^as 
in an measurement typically the real frequency / is on the horizontal axis and not the angular frequency oj. 

When the resonator is cooled to very low temperatures where k^T ~ Hcor, the classical description breaks down as 
the quantized energy-level structure (Fig. n4c) becomes important. Semi-classically, the thermal and quantum effects are 



combined by replacing the force noise of Eq. |5^ with the Callen and Welton equation [ ]204{ |: 

^mo) 1 / fwj 

■ - rux) coth 

Q 2 \2k^T 



— Amoj 1 / h(jj . 

S F„F„ (w) = — — ■ - hoj coth — — , (57) 



so that for g » 1 : 



(u ) - ui ■ coth « T = In 

° \2ksTj k^ 



' {U ) + Mg 

<m2> - ul 



(58) 



Then by inserting the first part of Eq. |8| in Eq . |5^ the resonator temperature is obtained^. The dependence of Tr on 
the environmental temperature is shown in Fig. |16|. At high temperatures {k-^T » huR) the resonator temperature is the 
temperature of the environment: Tr - T. At zero temperature the resonator temperature is determined by the quantum 
fluctuations: Tr = i^wr/^b- 



In thermal equilibrium, the energy levels of a harmonic oscillator have occupation probabilities that are given by [ 20 1 1 : 



e 'bM M - e 'bJ- j. (59) 



The average thermal occupation is « = 2,^0 ~ [exp(/iWR/^B7^) - 1] [201 1, which equals k^TRltuDR - The insets 



in Fig. g_6|show the occupation probabilities at three different temperatures T. At low temperature the resonator is in the 
ground state most of the time. At ^b 7^ = In 2 Hlor the probability finding the resonator in the ground state is exactly 50% 
and the average occupation is « = 1 . At any non-zero temperature, there is always a finite probability to find the resonator 
in an excited state. With the statement that "the resonator is cooled to its ground state" one actually means « < 1 . In Sec. 



3.5 we will explain in detail how this goal can be achieved using different cooling techniques, but now we will focus on 



the detection of the resonator position, in particular on the role of the detector, backaction and the standard quantum limit. 



-'Note, that some authors use Eq. £^ instead of Eq. ^ the definition of Tr, which implies that when (u^) = u^, Tr = 0. This in contrast to the 
definition used here where Tr = ishwR/k^ for (u^) = Ug. In the latter case Tr is not the actual temperature, but it is a measure for the (quantum or 
thermal) fluctuations. 
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3.4. Backaction and quantum limits on position detection 

Since the discovery of the Heisenberg uncertainty principle in 1927 it is known that quantum mechanics imposes 
limitations on the uncertainty with which quantities can be measured. This was first discovered for single measurements 
of conjugate variables, such as the position u and momentum p of a particle, or the components cr^, cr,. and <j, of the spin 
of a spin- 1/2 particle. 

When a system is initially in a superposition of the eigenstates of the operator corresponding to the quantity that is 
measured, a strong measurement gives one of the eigenvalues as the outcome |Q. The probability for measuring a value 
is |c,p when - 2, c, |yu,) was the expansion of the original state in the basis of eigenstates of the operator ji, with 
fi \ni) = Hi \ni). At the same time, the wave function collapses into the state corresponding to the measured value yu,: 
li/f) lyu,). To obtain the probability of a certain outcome of a measurement of a different quantity v, the state |yu,) has 
to be expanded in the basis |v,). The uncertainties in h and v satisfy: A/^ ■ Av > |[//, v]/2/|, where [//, v] - fiv - vfi is the 
commutator of fi and v. For fi - u and v - p, this yield the Heisenberg uncertainty principle for position and momentum: 
Au ■ Ap > h/2. Note that quantum mechanics does not forbid to determine the position with arbitrary accuracy in a single 
measurement. 



Most measurements are, however, not single, strong measurements, but weak continuous measurements instead [ 205 1 . 
As a measurement of the position disturbs the momentum of the resonator, a subsequent measurement of the position 
after a time Af inevitably is influenced by the previous measurement. This backaction is therefore important in continuous 
linear displacement detectors. In an experiment, backaction results in three different effects on the resonator: a frequency 
shift, a change in damping, and a change in the resonator temperature. If the resonator temperature is lower than the bath 
temperature, backaction has led to self-cooling, whereas a higher temperature indicates that there is a net energy ffow 
from the detector to the resonator. 

A way to circumvent backaction is to perform measurements on the position squared. As we show in Appendix ^ 
such a square-law detector probes the energy states of the resonator and these are not disturbed by the measurement itself 
as the energy operator equals, and hence commutes with, the Hamiltonian of the system. Therefore, this measurement 
scheme is called a quantum non-demolition (QND) or a backaction evading (BE or BAE) measurement [Q]. 

In the following section, continuous (linear) detectors and their backaction are discussed in detail and the quantum 
limits on position detection are explored. The analysis presented is largely based on the work by Clerk and co-workers 



[ 206| , 207]. We start with an analysis of generic linear detectors and what the effects of backaction are. Then we discuss 
three different routes to arrive at the quantum limit: the Haus-Caves derivation, the power-spectral density method, and 
the optimal estimator approach. These all have different ranges of applicability and rigor but lead to the same conclusion: 
the quantum limited resolution for continuously monitoring the resonator position is approximately the zero-point motion 
«{)■ 

The sensitivity and resolution of a detector will be important concepts in the following discussion. These notions, 
however, sometimes lead to confusion and, before defining them mathematically, we will first discuss the similarities 
differences between them. Both are a measure of the imprecision with which the position is measured. A noisy detector 
has a bad sensitivity and a bad position resolution, whereas a good detector has a good sensitivity and resolution. Coun- 
terintuitively, this means that the latter detector has lower sensitivity values, although one often says that it has a higher 
sensitivity. This is the first point that leads to confusion. The second is the difference between sensitivity and resolution. 
This can be understood as follows: Consider a resonator that is standing still and that is coupled to a detector that in- 
evitably introduces noise in its output signal. When measuring the position for a short period, the inferred position has 
a large uncertainty due to the detector noise. By measuring longer, the noise averages out and the uncertainty decreases. 
This uncertainty is the resolution of the detector and is measured in units of length. It thus depends on the duration of the 
measurement Af . For white detector noise, the resolution improves (i.e. its value decreases) as 1 / VAt. The proportionality 
constant between the resolution and the measuring time is the sensitivity, and it has the units of m/ VHz. In the discussion 
of Fig. [l^ it will be shown that the sensitivity is easily extracted from the noise spectrum of the detector output. 

3.4.1. Continuous linear detectors 

A continuous linear position detector gives an output that depends Unearly on the current and past position of the 
resonator Figure O shows the scheme of a generic linear detector. Note again, that the analysis is for a generic continuous 



linear detector and quantum effects only come into play in Sec. 3.4.2 We will discuss all its elements step by step. The 
output signal of the detector^ is related to the displacement by: 

v(t) = AAy(t) ® u(t) + v„(t). (60) 

Here, i,, is the responsivity of the detector, with /ly(f) - for f < due to causality. v„ is the detector noise and A is 
a dimensionless coupling strength. The coupling between the resonator and the detector is an important element when 



Continuous linear detectors are usually (implicitly) assumed to be time-invariant 1 208 1 . Unless stated otherwise, this is also assii rnpd ialhis Report. 
Examples where the linear detector is not time-invariant are frequency-converting, stroboscopic and quadrature measurements 1 208 , 4 N, |27| |. 
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Figure 17: The detector is coupled (via A) to the resonator. Its output v(;) depends linearly on u(t) with a response function A,, but also contains detector 
noise v„. The detector exerts a backaction force Fba = A^jg, on the resonator that contains a stochastic part Fbam and a linear response to ii(t). Both 
contributions add up with the thermal or quantum force noise F„. The resonator displacement u is obtained via the transfer function Hf; = H ho I hi- 




frequency 



Figure 18: The resonator displacement u is measured by the detector. The detector also adds imprecision noise v„ . The sum of this noise and the physical 
displacement v is recorded, (a) shows schematically time traces of the resonator displacement u, imprecision noise v„ and the apparent displacement v. 
Often, one is interested in the noise spectrum of the detector output, especially when the resonator displacement is its Brownian motion. This can be 
done using a spectrum analyzer and (b) shows that from a resulting spectrum the resonance frequency f^, the quality factor Q, the detector noise floor 
S„„„„ and the force noise S f„f„ are readily extracted. The signal-to-noise ratio is the height of the resonance peak divided by the height of the noise 
floor. It is therefore advantageous to have a large quality factor as this leads to a larger signal-to noise ratio. 



discussing different detectors in Sec. The choice of the coupUng A is sHghtly arbitrary as it could also be incorporated 
in Ay. However, in most detectors it is possible to make a distinction between the coupling to the resonator and the output. 
Different types of detector will be discussed extensively in Sec. ^ but in the case of an optical interferometer, v represents 
the number of photons arriving at the photon counter, in a single-electron transistor it is the current through the island, 
and in a dc superconducting-interference-device detector it represents the output voltage. The difference between u, v 
and Vn is illustrated in Fig. |l8|a. Note, that often the detector responds instantaneously to the displacement. In that case 
AAv{t) = dv/du (5(f) and AA{a)) = dv/du. The frequency response of the resonator is then flatQ 

To calculate the sensitivity of the detector, the noise at the output of the detector, v„, is refeiTed back to the input 
using the known response function Ay and gain A. This yields the displacement noise u„ at the detector input. The 
equivalent input noise PSD is S u„ii„ - S /|A/l,,p. This power spectral density is an important parameter that characterizes 
the detector. In nanoelectromechanical experiments, this noise floor is usually determined by the classical noise in the 
electronics of the measurements setup. In optical experiments, however, the noise floor can be shot-noise limited; quantum 
mechanics now sets the imprecision of the experiment. We will discuss the quantum limit in the next subsection. Note 
that for a flat frequency response u„ is simply given by v„/AAy. 

The detector does not only add noise to the measured signal, it also exerts a force FsAit) on the resonator This is the 
so-called backaction force. Backaction, in its most general definition, is the influence of a measurement or detector on an 
object. The detector backaction is a force on the resonator. This can be seen from the following argument: when there 
is no coupling between the resonator and the detector, i.e., A - Q, the Hamiltonian describing the total system is the sum 
of that of the Hamiltonian of the oscillator and of the detector Hamiltonian. When these are coupled (A 0), there is an 



interaction Hamiltonian of the form //,„, = -A<^detU- The backaction force is then Fba — -dHi„t/du = AfJ^/^,, [ [206[ ] . The 
backaction force is thus also proportional to the detector-resonator coupling A. The quantity O^c, is related to one of the 
internal variables of the detector; for example in an optical cavity it is proportional to the number of photons, in a SET to 
the electron occupation, and in a SQUID to the circulating current. 
The backaction force has three different contributions: 

• A deterministic force that is independent of the displacement. This changes the equilibrium position of the resonator. 



-^When the output of the detector has a delay time r, AAy(t) = dv/du S{t - r) and AAy(a)} = dv/du exp(-((t)r). 
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Figure 19: (a) The position resolution of the detector in units of the amplitude of the Brownian motion (uj) ' (classicaLiimit) for different coupling 
strengths A. The solid line is the total resolution; the dotted line is the contribution of the backaction force noise (Eq. p.5[ ) and the dashed line the 
contribution of the imprecision noise (Eq. |b.6| ). The total resolution is optimized at A = A^pi 0.0355. The resolution has been calculated for Qo = l(f, 

/iv = I, Af = 0, 5v„i.„ ■ ^0 = S<t,j„„iij„„ = S F„F„ and So^,,,,!., = 0. 



Without loss of generality it can thus be set to zero and it is not considered here. 

A force that responds linearly to the displacement: Fba.u = AAf{t) (g) Au{t). This changes the effective resonator 
response from to //^ , where //^"' - H^} + A^Afiai). An example of this force is the optical spring that we will 



encounter in Sec. 3.5.3. This part of the backaction can lead to cooling, see Sec. 3.5.1 [ 



• A stochastic force FBA.n = A(i>det,n that is caused by the fluctuations in the detector, ^cWt.n- Note, that this force 
noise and the imprecision noise v„ can be correlated, i.e., S'o^j^,,,,,, ^ 0. In any case, this contribution tries to heat the 
resonator since it adds up to the original thermal force noise.. 

In summary, this whole process of action and backaction is the one shown schematically in Fig. |l^: The resonator position 
is coupled to the input of the detector, which adds imprecision noise v„ and exerts force noise on the resonator. For small 



coupling A, the statistical properties of (i>ciet,n and v„ are independent of the resonator displacement [ 206 , 207 1 . 

One way to analyze the data is to measure the detector output using a spectrum analyzer. This way, information 
about the resonance frequency, quality factor, and imprecision and force noise PSD can easily be obtained, as illustrated 
in Fig. |l8|b. However, in a linear detection scheme one is usually interested in measuring the resonator position as 
accurate as possible, without perturbing the resonator considerably. Using optimal-control and estimation theory, the best 
estimate u for the resonator displacement in the absence of the detector, m,, is found. The resolution of the detector is 

Am - (K[uj - ip-fj , as explained in detail in Appendix]^ It quantifies the difference between the displacement that the 
resonator would have had when it was not measured and the one reconstructed from the detector output. By rewriting 

vl/2 



the resolution as Am - (E[(m2 - m2) + (m2 - m2)]) it becomes clear that there are two contributions |]209|]: the first one. 



which we name Am„, indicates how well the realized displacement is reconstructed from the detector output, whereas 
the second term, Am^a, quantifies the difference between m,- and the realized displacement, i.e., how much the motion is 
perturbed. Heuristically, we can understand that the first term is due to the imprecision of the detector, whereas the second 
is due to the backaction. 

The resolution is plotted in Fig. |l9|a as a function of the coupling strength A. In experiments this coupling strength is 
an important parameter and Sec. Q we will give typical numbers for the different detection schemes. We now discuss the 
general features of Fig. |l9|in more detail. With a low coupling Am is large (i.e., the detector has a low resolution) because 
of the large imprecision noise contribution Am„ (dashed line) of the detector. The backaction contribution is very small. 
An increase of the coupling reduces Am because by increasing the coupling, the mechanical signal becomes larger whereas 
5, remains the same, so the relative contribution of the imprecision noise decreases. The increase of the coupling also 
raises the backaction contribution, which still remains small. The resolution improves with increasing A up to the point 
where the optimal value A - Aopi is reached. A further increase of A makes the backaction force noise dominant, driving 
the resonator significantly, thus yielding a higher Am. 

3.4.2. The Haus-Caves derivation of the quantum limit 

The system analysis of the linear detector discussed above is valid for any - quantum limited or not - linear detector. 



An elegant way of deriving t he qu antum limit of a continuous linear position detector was given by Haus and Mullen piQ ] 



and was extended by Caves [ 208 1 . They consider the situation where the input and output signal of the detector are carried 

32 



by single bosonic modes, fl„ and a,, respectively. When the "photon number gain" of the detector is G = {a\,a^)l{a\,au), one 
might think that the modes are related to each other by a,, = ^^cii,. This is, however, not valid as this gives [ay, al] - G 
instead of the correct value, 1 [208, 207]. The actual relation is — y/Ga,, + v„. Here, v„ represents the noise added 
by the amplifier This operator has a vanishing expectation value ((v„) = 0) and is uncorrected with the input signal 
([ouyVn] - [ai,v„] - 0). Requiring [fly, ai] - 1 yields [v„,v,',] - I - G for the commutator and, more importantly, 
Aal > GAal + i|G - 1| for the noise in the number of quanta of the output mode [E07|. The first term is the amplified 



input signal (i.e., the resonator motion) and the second one is the noise added by the amplifier In the limit of large gain^ 
(G » 1), the equivalent input noise of the detector is A(flf*' )^ = Aa^/G - Aa^ > i. This means that a quantum-limited 
detector adds at least half a vibrational quantum of noise to the signal. 

As pointed out in Ref. [207| most practical detectors cannot easily be coupled to a single bosonic mode that carries 
the information of the resonator to the detector, because there is also a mode that travels from the detector towards the 
resonator. Therefore, the linear-system analysis at the beginning of this Section is used to further explore the quantum 
limits on continuous linear position detection. 



3.4.3. A quantum-limited detector 

In the previous discussion, no constraints were enforced on the detector noises Q>dei.n and v„ . If both noise contributions 
could be made small enough, the resolution would be arbitrarily good. This is unfortunately not possible. It can be shown 
that the power spectral densities must satisfyFI: 



or, equivalently, when this is referred to the input: 



(61) 



(62) 



These constraints enforce the quantum limit of the linear detector and should be considered as the continuous-detector 
equivalent of the Heisenberg uncertainty principle: Accurately measuring the position results severe force noise and vice 
versa. 



Clerk et al. [207] continue now by finding the gain where the total added noise at the input, i.e., 5„„„„(w) + 
\H'j^\^S Fba„FbaS'^)^ is minimized. As they already point out, this is not entirely correct because at every frequency a 
diff'erent optimal gain is required. Usually, only the optimal gain at w = is used and then the magnitude of the signal 
and detector noise are equal at that frequency. In that case, the imprecision noise and the backaction-induced displacement 
provide exactly half of the total added noise [207 3^ However, by optimizing the total resolution, the true optimal 
gain is found, see Fig. The resolution is optimized at A = 0.87 and reaches a value of 0.81 times the zero-point 
motion. 

All three methods (the Haus-Caves derivation, the total added noise at the input, and the optimal estimator) indicate 
that the detector adds about the same amount of noise as the zero-point fluctuations of the resonator itself. 

3.5. Cooling 

To prepare a nanomechanical system in the ground state, the thermal occupation of its normal modes should be 
minimal. The most direct approach is to mount an ultra-high frequency (/« > 1 GHz) resonator in a dilution refrigerator 
{T < 50 mK) so that « < 1. Such a resonator will, however, have a very small zero-point motion and the readout of tiny 
high-frequency signals at millikelvin temperatures is difficult. An alternative approach is to perform the experiments with 
lower-frequency resonators and/or at higher temperatures. The thermal occupation is then higher than one and cooling 
techniques have to be used to reduce the temperature of the resonator Tr well below the environmental temperature T . 
Figure |T] and Table || show how recent experiments are approaching the limit n < 1 over a range of frequencies that spans 
seven orders of magnitude. 

In most experiments either active feedback or sideband cooling is employed. In the former case, which is discussed in 
detail in Sec. 3.5.2, the position of the resonator is measured and the detector signal is fed back to the resonator to damp 
its motion. With sideband cooling (Sec. 3.5.3) the resonator is embedded in an optical or microwave cavity. Phonons 
can be removed from the resonator by up-converting a red-detuned photon to the resonance of the cavity, thus cooling the 
resonator Finally, also other cooling mechanisms like bolometric (photothermal) [^ ^ 211 1 and backaction cooling 
[B4I (see Sec 3.4) are used. We will now discuss these cooHng mechanisms in more detail. 



^*In the opposite limit wliere the detector does not have any net gain, i.e., G = 1, no additional noise is required by quantum mechanics. 
^^Here, it is assumed that the meas urement of v(t) does not result in an additional force noise on the resonator and that the detector has a large power 
gain. For more details, see Ref. [EOTI]. 
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Table 6: Overview of resonator temperature and cooling of recent experiments with micro- and nanomechanical resonators. The table shows the 
resonance frequency fa, the temperature of the environment T, the minimum resonator temperature T™'" and the corresponding number of quanta. The 
numbers of the first column correspond to the experiments hsted in Table hi 
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Figure 20: The quantum limit for a continuous linear position detector using the optimal estimation method. The resolution of the detector in units of the 
zero-point motion is plotted for different coupling strengths A. The solid line is the total resolution, Au; the dotted and dashed line are the contribution 
of backaction, &uba, and imprecision noise, Ah„, respectively. The inset shows that the total resolution depends on the cross-correlation coefficient 
^ FgAn'iBAn FB^„FBAn ' '^^ "n "n ) ' ' ^hc main paucl is calculatcd with an optimal cross-correlation of 0.36. For all values of A the total added noise is 
slightly less than uq. 



3.5.1. Backaction cooling 

To measure the displacement, the resonator is coupled to a detector As shown in Sec. 3A , this influences the resonator, 
and, in particular, this backaction adds force noise and it can damp the motion, which can lead to cooling. The damping 
rate of the resonator is then increased fr om it s intrinsic value yo - wo/Qo to_y« = yo +2bA' where jba = -cj^AHmlAf]/ a> 



is the damping induced by the detector [206|. The resonator temperature i; 



212]: 



7()T + JbaTba 
ro + ysA 



(63) 



The resonator temperature is thus the weighted average of the environmental temperature T and the so-called backaction 
temperature of the detector Tba- For strong resonator-detector coupling {jba ^ To) the effective temperature is Tr = Tba- 
When this is below the environmental temperature, the resonator is cooled by the backaction. Eq. 56 shows that the 
backaction temperature is determined by the force noise exerted on the resonator: Tba - S Fn^^FBAiJ^^BmyBA- Both 
S Fba„Fba„ ^nd Jba are proportional to so that Tba is independent of the resonator-detector coupling A. In other words, 
the backaction temperature is an intrinsic property of the detector 

Although it might not be immediately clear, this cooling mechanism corresponds to the usual notion of cooling: 
Cooling is done by coupling something to something else that is colder In the case of backaction cooling the cold object 
is the detector. Because the resonator is not actively cooled, but only brought into contact with the detector, backaction 
cooling is therefore also called "self cooling" or "passive feedback cooling". 



3.5.2. Active feedback cooling 



It was shown in Sec. 3.3 that the resonator temperature is proportional to its random motion. By reducing that motion 
the resonator gets cooled. One way to do this is using feedback. When the position of the resonator is measured and fed 
back to it via an external feedback loop, the motion can be amplified or suppressed. Feedback control was already used 
to regulate the motion of soft cantilevers [213, 214, 215| for magnetic resonance force microscopy [216|, when it was 
realized that this technique can also be used to cool a mechanical resonator towards its ground state [217|. Actually, the 
lowest resonator temperature to date {Tr = 1 .4 /vK, see Table ^ has been reached using this cooling method ||64|]. 

Feedback systems are usually analyzed within the linear system representation. Figure ^ shows a schematic of 
the process. The resonator, with frequency (Dq/Itt and Q-factor Qo, is driven by the thermal force noise F„{t) and its 
displacement u{t) is detected. The detector output contains not only the displacement but also imprecision noise M„(f)-0 
The apparent position v is the signal at the output of the detector and this is thus the sum of the physical displacement and 
the detector noise: v = u + u„. The information contained in v is used to apply a force FfB to the resonator that damps its 
thermal motion^. The relation between the feedback force and the apparent position is described by the linear system or 



-**This assumes that the effective spring constant kg - koA^Re[Af] does not change considerably. As shown in Sec. 3.5.2 it is in general much more 
difficult to alter the spring constant than to alter the damping. 

^'in this section we assume that the detector has unit gain, i.e., AAy(ii>) = 1. Then the imprecision noise at ffie output, v„, and the noise referred to the 
input, u„, are the same. 

Although this term was never important in active 



In principle, also the backaction force noise of the detector has to be included [213 



52, 21 



feedback cooling expeiiments so far, ffie effect can be included by using ffie damping and resonator temperature of the coupled resonator and detector 
instead of the intrinsic ones of ffie resonator alone. 
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Figure 21: Overview of meclianical resonators with low temperature or occupation, compiled from Table q The top panel shows the experiments with 
the lowest occupation numbers. These arejEached using conventional cooling (gray) and active feedback and sideband cooling (black). The dotted line 
is located at n = 1. As discussed on page M, once the thermal occupation is below this value, the resonator is cooled to the ground state. The bottom 
panel shows the starting temperature T (gray) and the final temperature 7™° (black) achieved by groups that h ave ar- tivply cooled their resonator below 
100 mK. Note the diagonal line in the LIGO experiment which is due to a strong optical spring eft'ect (see Sec. 3.5.3), so both the resonator temperature 
and frequency change when the laser power is increased |UCf| . 



filter with transfer function HfB{co)- The output of the filter is multiplied by a selectable gainQg. One can think of it as 



a knob to crank up the gain of an amplifier. This forms a closed- loop system 1 199 213] with the following equations of 
motion: 

mu(t) + ma)Qu{t) I Qo + mci^u{t) = F„(t) + Ffsit), (64) 

FpBit) = mcolg ■ hpBit) ® lu(t) + u„(t)]. (65) 

The presence of the feedback results in a different displacement for a given thermal noise realization F„(f). The feedback 
thus changes the resonator response from //« to the closed-loop transfer function H'^, given by: 

H'r' = -gkoHpB, or H'g = " . (66) 

Comparing this with the response of the resonator itself, cf. Eq. ^ shows that the real part of HfB modifies the 
resonance frequency from coo to ojr - uq ^JT^^gRe[Hp^(a>y\, whereas the imaginary part alters the damping rate from yo 
to jR = Jo - ^WQlm[///rB/w]. Using this closed-loop transfer function, the PSDs of the physical (i.e., the real resonator 
displacement) and observed displacement (i.e., the detector output) are obtained: 

Sf„fJimajlf+g^\HfB\^Su„,„ 
5„„(w) = 2' (67) 

|l -( — ) + -Jr—- gHFB(c^)\ 



SM = : ^ — ^^^^ . (68) 



The resonator displacement PSD S uu shows that the resonator indeed responds to the force noise with the modified transfer 
function //^ instead of Hr. The force noise that drives the resonator (i.e., the numerator of the right-hand side of Eq. |&^ ) 



The gain g and filter Hps are defined such that \HfB(pio)\ = 1- Both g and Hps are dimensionless. 
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Figure 22: Linear system representation of the active feedback cooling sclieme. The resonator displacement u is converted by detector to its output 
signal V. This adds imprecision noise h„ and the sum of this noise and the physical displacement is the signal that is measured, for example using a 
spectrum analyzer. This signal v is also fed back to the resonator to attenuate the Brownian motion. In the feedback loop a filter with response HfB and 
a variable gain g ■ ko are included. The resulting feedback force FfB adds up with the (thermal) force noise F„. The resonator's response to the applied 
forces is determined by its transfer function Hn = Huo/ko- 



Still contains the original contribution S f„f„ , but now it also has a contribution due to the imprecision noise that is fed back 
to the resonator Since the latter is always positive, the feedback loop adds additional force noise to the resonator The 
apparent position PSD is also modified: It is not simply the sum of S uu and S u„u„ because the feedback creates correlations 
between the imprecision noise M„(f) and the actual resonator position u{t). 

So far, the analysis was general for any linear feedback system and different implementations of the feedback filter 
HfB are possible, each with their advantages and drawbacks. Using optimal control theory, the best feedback filter can in 



principle be found [213|. In practise, often simpler, but therefore suboptimal, filters are used. The PSDs of the true and 



apparent resonator displacement are plotted in Fig. ^ for the two simplest feedback schemes; 

• Velocity-proportional feedback where the measured displacement is used to apply a velocity-dependent force on the 
resonator with /i/r^Cf) - -lJq^ -d/dt, HfB = -/w/wq- In this case the dampingrate increases from yo to yo'(l +^6o) 
as indicated by Eq. ^ Figure ^3| a shows that at low gains, 5,,,, is lowered and thus that the resonator is cooled. 
However, when the gain is increased further, the tails of S uu start to rise as the detector noise (the second term in 
the numerator of Eq. ^7|) is fed back into the resonator Above a certain value g = gmin, too much detector noise is 
fed back to the resonator and the resonator temperature increases again. Figure^ shows the resonator temperature 
as a function of the feedback gain. The minimum resonator temperature in the limit g » Q^^ is 



' i?,min 



mcolT- 2T I 

A 1-^^ ""«« = -^=' for g = gmin = VSNR/go. (69) 
\ ksQo VSNR 



The minimum resonator temperature is thus set by the signal-to-noise ratio (SNR = S mi((i>o)g=o/S u„u„ ) of the original 
thermal noise peak and the detector noise floor, as was illustrated in Fig. |l8|b. Finally, note that, unlike for backaction 
cooling, the resonator temperature does not saturate at a certain value when g — > oo. Eventually more and more 
noise is added and the resonator temperature keeps on increasing with increasing gain. 

• Displacement-proportional feedback where the displacement is directly fed back to the resonator, which is charac- 
terized by hfsit) - -6(t), HfB - -I- This changes the spring constant from maj^ to mco^^il + g) and the resonance 
frequency to too ■ (1 + g)^^^- The stiffening of the resonator reduces its thermal motion and hence its tempera- 
ture, but to achieve the same cooling factor as with the velocity-proportional feedback the gain should be Qo times 
larger This, however, also feeds back much more detector noise to the resonator in the usual situation where 
2o ^ 1- Cooling can therefore only be achieved when the SNR is large. Figure p3| b shows that only heating in- 
stead of cooling is achieved for the choice of the parameters used to perform the calculation. For high-Q resonators 
velocity-proportional feedback is superior to displacement-proportional feedback. 

It is important to note that the feedback creates correlations between the resonator displacement and the detector 
imprecision noise, which lead to a substantial change in the shape of the noise spectra. Figure ^ shows calculated noise 
power spectral densities for different gains. Without feedback (i.e., g = 0) the spectrum of the apparent motion is simply 
the sum of that of the harmonic oscillator, 5„„, and a constant background level due to the noise, 5„„„„. The Figure shows 
that, even though the peaks in S uu only become broader or shift in frequency, the peaks in S „, becomes distorted when 
the feedback is applied. The resonator temperature is no longer simply given by the area under the peak in S ,,,, due to 
the abovementioned correlations between u„ and u. For large gains it is even possible that the peak in the spectrum S w 
changes into a dip (Fig. Ha). A way to circumvent this problem of determining the resonator temperature, is to use a 
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Figure 23: Feedback cooling of a resonator with Qo = 100 using velocity-proportional (a) and displacement-proportional (b) feedback. The top panels 
show the PSD of the observed displacement (i.e., that of the detector output) S,.r and the bottom panels show the real (physical) displacement PSD, S uu- 
The gain is stepped from g = (black) to ^ = 0.2 (light gray). These plots are calculated with 5 = 10^ and all PSDs are scaled by 5 /r„f„ /(mWg)^. 



second detector to measure the resonator motion [pq]. The noise of the second detector is not correlated with that of the 
resonator and the measured PSD is again the sum of a constant background and the resonator motion. 

Active feedback cooHng experiments have mainly been done on resonators in the kHz range, where one can simply 
measure the position, differentiate and feed the resulti ng sig nal back to the resonator. Cohadon et al. demonstrated the 
first active cooling of a mirror using feedback control [219|. They use a high-finesse cavity with a coated plano-convex 
resonator as the end mirror A feedback force is applied using a 500 mW laser beam. An acousto-optical modulator is 
used to control the exerted radiation pressure on the resonator. Next, Kleckner and Bouwmeester actively cooled a 12.5 
kHz AFM cantilever from room temperature to 0.135 K [^]. A tiny plane mirror attached to the cantilever served as the 
movable end mirror of the cavity and again a second high-power laser was used to apply the velocity proportional force. 
Arcizet et al. cooled a millimetre-scale resonator in an optical cavity to 5 K by applying an electrostatic feedback force 
on the resonator [^^. By using a piezo element to apply feedback to an ultrasoft silicon cantilever cantilever cooling 
from 2.2 K to 5 mK was demonstrated by Poggio and coworkers [^. Feedback cooling has also been demonstrated with 
the gravitational-wave detectors AURIGA and LIGO [p4|]. In the former experiment the 2 ton heavy detector was 
cooled from 4 K to 0.17 mK. The motion of the bar is measured capacitively using a resonant electric circuit coupled to 
a SQUID amplifier The output of the amplifier is put through a low-pass filter to create the 7r/2 phase shift required for 




Figure 24: Resonator temperature for velocity-proportional feedback vs feedback gain g for a resonator with 2o = 100. The dotted line is for 
^qS u„ui,/S F„F„ = 10^ which gives a signal-to-noise ratio of 10- and the solid line is for k^S u^^u^JS f,^f,^ = 1, with SNR = 10*. For low gains the 
two curves overlap. When g > . 1. the detector noise already starts to heat the resonator for the low-SNR curve, whereas the high-SNR curve still 
decreases. As predicted by Eq. 6Sl the minimum temperature is lower in the latter case, and this occurs at a higher gain. 
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Figure 25: (a) schematic overview of an optical cavity. Tlie cavity is driven via an input laser with frequency a>j. The left miiTor is fixed, but the 
right mirror is a mechanical resonator that can move. The resonator displacement u determines the length of the cavity and thus the cavity resonance 
frequency 01^. A part of the circulating power is transmitted by the left mirror to a detector. The linewidth of the cavity is From T. J. Kippenberg, K. J. 
Vahala, Science 321 (2008) 1 172-1 176. Reprinted with permission from AAAS. (b) When the cavity is driven on its resonance, the intensity inside the 
cavity is largest, a detuning reduces the intensity. A displacement of the resonator shifts the cavity resonance (purple) and changes the intensity of the 
light inside the cavity (orange and purple dots), (c) The displacement dependence of the radiation pressure Fmj. When the resonator oscillates, the force 
reacts with a delay due to a finite value of as indicated by the ellipsoidal trajectories. 



velocity proportional feedback and then injected into the electronic circuit. In the latter experiment the center-of-mass 
motion of the four mirrors of a Fabry-Perot cavity with 4 km long arms is reduced from room temperature to only 1 .4 fik. 
This is done by adjusting the servo control of the mirrors, which creates a velocity-proportional feedback force in the right 
frequency range without affecting the very high signal-to-noise ratio. 

When the resonator frequencies are high, say above 1 MHz, delays in the feedback circuit start to play a role. The 
force is then applied when the resonator has already advanced and a purely velocity-proportional feedback will have a 
displacement-proportional component, degrading the cooling performance. The effect of a delay is even more dramatic 
when it equals half the resonator period, so that the Brownian motion is actually amplified instead of attenuated. Further- 
more, the bandwidth (or sampling speed in the case of a digital filter) of Hfs should be at least a few times (Dr, which is 
often not an issue for a>R < I MHz. However, active feedback cooling has been reported for MHz resonators by either 
adjusting the delay [ |5^ ] of the signal, or using a mixer circuit to down-convert the mechanical signal to a lower frequency 
[[70|]. Finally, note that feedback can also have other purposes than cooling. Examples are the regulation of the motion of 
soft cantilevers as mentioned at the beginning of this section, and feedback can be used to modify or null the nonUnearity 



of a resonator [ 22C ]. 



3.5.3. Sideband cooling 

Another commonly used cooling technique is sideband cooling [221, 222, 223 1. In this technique the resonator is 
embedded in an optical [|3[|l[|7[||,|3|,|3|,ra|7[|o[|3[^or microwave cavity @ ^ Js^ ||, These 
detection schemes will be discussed in detail in Sec. W but here a brief introduction is given. Figure [2^Thows a schematic 
drawing of an optical cavity where the right mirror is the mechanical resonator. Both mirrors have a low transmission 
so that a photon is reflected many times before it can go out through the left mirror, and then toward the detector. Such 
a cavity has many different optical eigenmodes, but here we focus on a single one and denote its resonance frequency 
by oJc- In analogy with the Q-factor and linewidth y« of a mechanical resonator, the cavity has an optical Q-factor Qopt 
and linewidth k. A laser sends light with frequency cOd into the cavity. This frequency can be different from the cavity 
resonance frequency Wc; the light is then detuned. Because the resonance frequency of the cavity is determined by the 
cavity length, a displacement of the resonator changes a>c. As illustrated in Fig. p5| b this leads to a change in the intensity 
(and phase) of the light in the cavity, which results in a change in the detector output. Optical cavities used this way are 
very sensitive position detectors for two reasons: first, it enhances the intensity of the light by a factor Qopi, and secondly 



it makes the intensity depend strongly on the displacement [222]. 



Each photon in the cavity carries a momentum pp^ - ha>d/c whose direction is reversed when it reflects off the mirror. 
Here c is the speed of light. The resonator thus experiences a kick of 2pph every time a photon reflects, the so-called 
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radiation pressure. A single round trip of a photon in a cavity of length L takes a time 2L/c, whereas the average time that 
a photon spends in the cavity is a:"'. The total transfer of momentum per photon is thus - HoJc/kL. The total force 
exerted on the resonator is also proportional to the number of photons present in the cavity, rtc, and equals Frad = hujcndL. 
Note, that tic is proportional to a: ' and to the input power. 

The radiation pressure depends on the displacement of the resonator This is illustrated in Fig. ^3^ : A small change in 
u changes the cavity frequency coc, which in turn leads to a proportional change in radiation pressure. This thus changes 
the effective spring constant of the resonator to - ko - dF^^^ldu; the so-called optical spring |224, 225, 226, 227, 228, 
^ , 22S , 230 1 . Similar to displacement-proportional feedback this can cool the resonator [p9| |64 , 23 1 1 and can even lead to 



bistability of the resonator position [ 225 ] . A much stronger cooling effect is, however, the fact that the number of photons 
does not respond immediately to a change in displacement, but that they can only slowly leak out of the cavity at a rate 
~ K. This was first realized [224| and demonstrated [232| by Braginskii and coworkers. When k » (Wq, a change in the 
displacement changes the number of photons instantaneously and Frad follows the orange curve in Fig. p5| :. However, a 
finite value of k causes a delay in the response of the radiation pressure as indicated by the ellipsoids in Fig. ^5| e. In the 
case of red-detuned driving (wj < w^. ; the ellipse is traversed counterclockwise) work is done by the resonator so that it 
loses energy, whereas for blue-detuned driving (w^ > cOc', clockwise trajectory) the resonator gains energy. The increased 
damping for re d detuning cools the resonator, as the backaction temperature associated with the detector is very low 
[ 233 , 234 , 223 ]. This process is called "dynamical backaction" and the cooling mechanism is called "sideband cooling" 
because the cavity is driven off-resonance, i. e., on a sideband^ | |236| . Note, that sideband cooling can also be described 
in the language of backaction cooling (Sec. 3.5.1 ) as the optical force responds with a delay to the displacement. This is 
thus a detector response Ay (Fig. [l^) with a real (the optical spring) and an imaginary part (the delay). 

The ultimate limit on the resonator temperature that can be reac hed with s ideband cooUng has been studied using the 
radiation-pressure Hamiltonian in the rotating- wave approximation [233 234 1: 



H - hAc^c + ha>oa a + hGouc Her + a). 



(70) 



Here, Gqm - uodtoddu is the optomechanical coupling rate and c ' (c) is the creation (annihilation) operator for a cavity 
photon (from (to) a photon with frequency w^). A - o),] - tOc is the detuning of the laser light with respect to the cavity 
resonance frequency. The quantum mechanical picture of sideband cooling is that a phonon together with a red-detuned 
driving photon can excite a photon in the cavity. This is a Ukely process because it up-converts the red-detuned driving 
photon to a frequency closer to the cavity resonance. This removal of phonons cools the resonator The opposite process 
is also possible: a cavity photon can emit a phonon and a lower-frequency photon, thereby heating the resonator. The 
rate of these two processes depends on the density-of-states of the cavity at u)d + wo and - wo respectively. If the 
detuning is at exactly at the mechanical frequency A = -wq and the cavity linewidth is small, the lowest temperatures are 
obtained. The process is analogues to the doppler cooling of cold atoms. In the good-cavity limit (w,,, > k) the lowest 
resonator occupation is Timm - ^ 1 [233, 234 , |222[] . Alth ough the good-cavity limit (also called the resolved- 

sideband regime) was reached a few years ago [|49[ p6|r|51|, |60|, |63 , 237 , ground-state cooling was not immediately 
demonstrated, due to the fact that the driving could not be increased to sufficiently high powers. The cooling power 
should be large enough to remove t he h eat coming from the environment to reach the ground state and the cooling power 
is proportional to the input power [223|. Table ^ and Fig. |1] show the final thermal occupation numbers that have been 
reached up to now. 

The effects of sideband cooling (i.e., frequency shift, change in damping, and cooling) were first considered in the 
context of gravitation al wave detecto rs, where the effect was more-or-less viewed as a technical point with limited ap- 
plications, see Refs. [224, 232, 238 [ and references therein. The first experiments with the aim of cooling mechanical 
resonators towards the ground state were reported in 2006 by Arcizet [ ^3| | and Gigan [ |3^ cooled resonators with frequen- 
cies of a few hundred kHz by modest factors of 30 and 10 respectively. These measurements were done using free-space 
optical cavities which were in the unresolved sideband regime (i. e., wh ere llir <sc k). Schliesser and co-workers used a 
mrrcotoroidal resonator vibrating at 58 MHz as cavity (see Sec. [4.1.3| ) and demonstrated cooling to 11 K [p7|], which 
corresponds to an occupation number of n = 4300. A much larger, gram-scale mirror was cooled to a much lower tem- 
perature Tr - 6.3 mK by the LIGO team [ ^ p^ . The resonator frequency is 1 kHz in this case (this is with the optical 
spring included, the bare resonance frequency is only 12.7 Hz). The thermal occupation n ~ 10^ is therefore higher than 
the abovementioned experiments by Schliesser. A major step forward was made in 2008 by the Kippenberg group, who 
were the first to reach the resolved sideband regime [p^. Their 73.5 MHz microtoroidal resonator has a cavity line width 
of only 3.2 MHz, placing this device deep in the good-cavity limit. The resonator was cooled to ~ 19 K. 

It was realized that a way to further cool resonators is to start at a lower temperature, so less phonons have to be 
refrigerated away to reach the groundstate. This approach was pursued by the group of Lehnert using superconducting 
stripline resonators (Sec. [4.1.5| ). Their first experiment was in the unresolved-sideband regime [^, but this was soon 



■A recent proposal uses a displacement-dependent cavity damping instead of the usual displacement-dependent cavity frequency [ 13^ ]. 
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superseded by measurements in the good-cavity limit [|5^, where a 0.237 MHz resonator was cooled from a 50 mK base 
temperature to 10 mK. Superconduting striplines are of course always operated at low temperature, but recently a lot of 
effort has been done to place optical cavities at cryogenic temperatures. Groblacher et al. precooled their microresonator 
to 5 K and subsequently cooled the 1 MHz resonator to 1.3 mK, where n - 32 [|^. Around the same time, Schliesser et 
al. cooled a resonator placed in a cryostat from 1 .7 K to 0.2 K where n = 63. Very recently, the first actively cooled device 
with an occupation number smaller than one was demonstrated by Teufel et al. [p3|]. They employed sideband cooling in 
a resonant superconducting circuit. 



3.5.4. Concluding remarks 

In the previous Subsections backaction cooling, active feedback cooling, and sideband cooling have been discussed. 
Although it might appear that these methods are unrelated, the converse is true: From the linear system representation of 
backaction (Fig. |l7|) and that of active feedback cooling (Fig. ^), it is clear that they are closely related. In the former 
case the delayed force, needed to achieve the largest cooling factors, is caused directly by the detector, whereas in the 
latter it is actively exerted by the experimenter The coupling A and the gain g play the same role in the respective pictures, 
and so do the detector force response Ap and the feedback filter Hfb- Also, the lowest temperature that can be reached is 
determined by the noise, <S>dei,n and u„ respectively. Note, however, that for backaction cooling the resonator approaches 
Tba when the coupling A — > oo, whereas in the case of active feedback cooling the resonator temperature diverges for 
§ — > oo, and the minimum temperature occurs at a finite value of the feedback gain. It was also shown that one can 
describe sideband cooling in a similar fashion: namely as the (delayed) response of the optical force to the displacement, 
and the same is true for bolometric forces where the delayed force is due to the finite heat capacity [ pi] , p4[ [21l| , 239 1. 
Finally, we no te th at there are cooli ng m echanisms in non-cavity systems, such as in superconducting single-electron 
transistors [32, 240 1 and double dots [241| that are formally identical to sideband cooling. 

Another important point is illustrated in the experiments done by Teufel and coworkers (lines 3 1 , 32 in Tables |l] and ||). 
When cooling the resonator, in this case using sideband cooling [|6|, the quality factor decreases. The resolution of the 
detector then degrades when the sensitivity stays the same because the resonator bandwidth jr increases. In other words, 
one has less time to average-out the imprecision noise. The opposite occurs when the damping of the resonator is reduced: 
the resolution is improved, but the resonator temperature increases significantly []6^. In particular, this happens in a cavity 
optomechanical system for a blue-detuned laser drive, but this effect has also been observed in single-electron transistors 
[^2[], superconducting interference devices [242] and in many other detectors. In this regime energy is transferred from 
the cavity to the mechanical system. The oscillation amplitude grows and the resonance becomes sharper as the resonator 
energy increases. When the coupling is strong enough, the damping rate can vanish (yo + Jba - 0) or even become 
negative. Beyond this instability the resonator exhibits self oscillations that are only bound by nonlinearities in either the 
resonator or the detector [ 224|, and this can lead to complex nonlinear dynamics. Theoretical work on the classical effects 
in this regime has been documented in Refs. [243, 244, 245, 246 1. Experimentally t his instab ihty has been observed in 
a va riety of micro and nanomechanical systems, including microtoroidal resonat ors [247 , 248 , 37| ], Fabry-Perot cavities 



3|, 228, 239 1, superconducting stripline resonators [ p6t and microspheres [249|. Devices that are self oscillating 



can be us eful f rom an application point of view as ultra-high frequency, low phase-noise oscillators [^1]], or as memory 



elements [250[. 
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Figure 26: Resolution (top panel) and sensitivity (bottom panel) of the experiments that ai'e Hsted in Table Q Experiments with optical (electrical) 
detection are shown in gray (black) respectively. The solid line indicates the standard quantum limit. As discussed in Sec. ^, detectors can have an 
imprecision resolution below the quantum limit at the price of a large backaction force noise. There is a trend that the resolution and sensitivity degrades 
(their value increases) with increasing resonator frequency. 



4. Detection methods 

In this section we discuss the main detection methods for mechanical motion in nanoscale systems or in mechanical 
devices which aim at reaching the quantum regime. Most methods involve a linear displacement detector, i.e., the output 
of the detector depends linearly on the displacement of the resonator. As shown in the previous sections, such schemes 
inevitably introduce backaction on the resonator position, which consequently leads to the fundamental limit on the 
sensitivity set by quantum mechanics. We will only briefly comment on quadratic (square-law) detectors, as they have 
not been studied to the same extent as linear detectors. In recent years, tremendous improvement in the sensitivity has 
been obtained using a variety of different detection methods. Using optical cavities, displacements sensitivities as small as 
10"^" m/ VHz (see Table have been reached; using mesoscopic electromechanical devices the sensitivity can be as good 
as 10"'^ m/ VHz. Fi gure Ma shows the resolution Au„ due to the imprecision noise S u^^u„ of the experiments listed in the 
Table. (For the exact definition of Am„ in this context we refer to Eq. B.7 in Appendix ^). Note, that this quantity does not 



include the effects of backaction force noise and can therefore be substantially smaller than the standard quantum limit as 
the Figure shows. This is most clearly demonstrated by the relatively low-frequency resonators that are read-out optically 
[ ^6[ , p3| , ^ p3| , |60[ |. Although the backaction force noise is present to ensure the SQL and thereby heating the resonator; 
at room temperature, this contribution is masked by the much larger thermal motion. Also in solid-state devices at low 



temperature, Am„ < mq has been achieved [|32| |66[ |68| and there an increase in resonator temperature due to backaction is 
seen 



The most straightforward method to analyze the data measured in an experiment is to record the output of the linear 



detector as a function of time; for a phot odetector in a cavity experiment (Sec. |4.1[ ) this would be an output current and for 



a dc SQUID position detector (Sec. |4.5| ) this is the voltage over the SQUID. This signal can be fed to a spectrum analyzer 
and from the measured spectrum, the resonance frequency and Q-factor can directly be obtained in the linear response. 
From the thermal noise spectra (cf. Fig. [T^) other parameters like 5„,„,^, S f,„iF,„i can be obtained as illustrated in the 
Figure. 

An important issue in measuring the resonator dynamics is the available bandwidth of the setup. In particular this 
holds for the solid-state devices that often have a high impedance (» 50 Q). The combination of high impedances and 
unavoidable stray capacitances can lead to RC times that are smaller then the resonator period 1 //«. There are several ways 
to circumvent this problem: frequency mixing, impedance matching using tank circuits, or using low-impedance devices 
such as superconducting quantum interference devices (SQUIDs) and microwave striplines. Yet another way relies on 
self-detection which yield dc information about the vibrational motion, e.g. rectification and spectroscopy measurements. 
These (non-linear) methods are mainly used in bottom-up devices. 

In the following sections, we will discuss different detection schemes focusing on the resonator-detector coupling 
and on the backaction mechanism. As we concentrate on detection schemes aimed at the quantum regime, we will, for 
example, not consider piezo-resistive detection schemes, which are very common in MEMS and for applications, but 
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Table 7: Overview of the sensitivity and resolution of recent experiments. The table shows the type of detector that was used, and the displacement- 
imprecision noise 5„„„„ from which the resolution Aei,, is calculated. The ratio of the resolution to the zero-point motion uq is also stated. The numbers 
of the first column coirespond to the experiments listed in Table N and Table M. 





Detector 


Su„u„ fm/ VMi 


Au„(fm) 


AUn/Uo 


Ref. 


1 


Magn. mot. 








1 


28 




2 


SET 


2.0 


658 


130 


1 


29 




3 


SET 


3.8 


113 


5.4 


1 


30 




4 


Opt. cav. 


2.0- 10^ 


4.8 ■ 10^ 


415 


1 


31 




5 


SET 


0.35 


5.9 


0.25 


1 


32 




6 


Opt. cav. 


4.0 ■ 10-* 


4.5 ■ 10-3 


0.62 


1 


33 




7 


Opt. cav. 


100 


38 


7.2 


1 


34 




8 Opt. cav. 1 


35 




9 


Opt. cav. 


4.0 ■ lO-'* 


4.5 ■ 10-3 


0.62 


1 


36 




10 


Opt. cav. 


3.0 ■ 10-^ 


0.53 


5.4 


1 


37 




11 


Piezoresist. 


39 


1.8- 10^ 


506 


1 


38 




12 


Opt. cav. 


0.2 


0.06 


8.3 


1 


39 




13 


Opt. cav. 


0.03 


9.5 ■ 10-4 


0.04 


1 


40 




14 


Opt. cav. 


1.0- 10^^ 


2.8 ■ 10"* 


778 


1 


41 




15 


APC 


2.3 


268 


29 


1 


42 




16 


Opt. cav. 


1.0- 10^ 


271 


2.7 


1 


43 




17 


MW res. 








1 


44 




18 


Opt. cav. 


2.7 ■ 10-^ 


2.3 ■ 10-'* 


1.8 


1 


45 




19 


Opt. cav. 


0.54 


49 


39 


1 


46 




20 


Opt. cav. 








1 


47 




21 


Opt. cav. 


45 


2.5 


25 


1 


48 




22 


Opt. cav. 


1.0- 10-^ 


0.05 


0.42 


1 


49 




23 


Stripline 


200 


2.5 ■ 103 


19 


1 


50 




24 


Magn. mot. 


12.8 


6.6 ■ 103 


339 


1 


51 




25 


QPC 


1.0- 10^ 


588 


20 


1 


52 




26 


SQUID 


10 


132 


51 




19 




27 


Capacitive 


3.0 ■ 10-5 


1.0- 10-'' 


0.34 


1 


53 




28 


Capacitive 


3.0 ■ 10-5 


1.2- 10-'' 


0.42 


1 


53 




29 


Opt. cav. 


1.0- 10^ 


9.5 ■ 10'' 


4.8 ■ 10^ 


1 


54 




30 


Opt. trans. 


72 


6.2 ■ 103 


232 


1 


55 




31 


Stripline 


600 


1.7- 103 


57 


1 


56 




32 


Stripline 


45 


696 


23 


1 


56 




33 


Capacitive 


2.0- 10"* 


1.9- 10^ 


8.5 ■ 103 


1 


57 




34 


Opt. cav. 


1.5 ■ 10-^ 


0.2 


30 


1 


JO 




35 


Opt. trans. 


40 


2.8 ■ 103 


76 


1 


59 




36 


Opt. cav. 


0.03 


0.18 


0.4 


[ 


60 




37 


Opt. cav. 


1.5 ■ 10-3 


0.34 


7.9 


1 


61 




38 


Opt. trans. 


0.04 


11.7 


2.4 


1 


62 




39 


Opt. cav. 


0.26 


61 


1.2- 103 


1 


63 




40 


Opt. cav. 


1.0- 10-3 






1 


64 




41 


Opt. cav. 


0.07 


4.0 


1.5 


1 


65 




42 


Stripline 


4.8 


15 


0.57 


1 


66 




43 


Opt. cav. 


0.64 


23 


1.6 


1 


67 




44 


Stripline 


1.2 


209 


8.3 


1 


68 




45 


Opt. cav. 


0.6 


12.5 


0.76 


[ 


69 




46 


Qubit 








1 


22 




47 


SQUID 


2.0 


24 


9.3 


1 


70 




48 


MW res. 


0.23 


92 


23 


1 


23 
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rarely used for nanoscale experiments [|8 251 1. The main disadvantage of this method is that it requires a large current 
that is dissipated in the resonator. 



4.1. Cavities 



Table 8: Overview of the coupling factors reached in optical and microwave cavity experiments. As explained in the text goM = dojc/du quantifies 
how much the cavity frequency changes with the displacement and Gom = gOM X "0 is the vacuum coupling rate. The cavity frequency aj^/ln lies 
in the GHz range for the microwave cavities and is ~ 10^ THz for the optical cavities. *: is the cavity linewidth. In the experiment on the last line, a 
superconducting qubit is used instead of a cavity. In this case, the value of vacuum Rabi rate is listed as an estimate of Gom [Ml- 



^ 8oM /MHz\ Gom ,,, . <^c . k , , , , „ * 

Group -— — — (Hz) — (GHz) — (MHz) Lqm (m) Ref. 

In \ nm / In In 2n 
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1 1 A 1 
i . io ■ lU 


U. i J 




4.91 


n AO 




■ iU 


-3 




50 


1 


TTT A 


4.0U ■ lU 


M 1 A 




5.22 


U./j 


L.Li 


■ iU 


-3 




56 


1 


TTT A 


-3 on 1 r\-2 


W.oO 




7.49 


Z.oo 




1 

L\) 


6 




56 


1 


Schwab 


/ .jU ■ lU 


U. io 




5.01 


n AO 


ooo 


1 

iU 


6 


t 
1 


25: 


^1 




8.40 ■ 10-2 


2.12 




7.48 


0.60 


89.0 


- 10 


-6 








Aalto 


1.0 


16.5 




7.64 




7.64 


10 


-6 


1 


253| 


NIST 


56 


226 




7.47 


0.17 


133 


10- 


9 




2541 


NIST 


49 


198 




7.54 


0.20 


154 


10- 


9 




23 


\ 


Vienna 


1.13 • 102 


206 


2 


82 


10^ 


120 


2.50 


- 10 


-3 




35 


1 


LKB Paris 


1.17-102 


75 


2 


82 


10^ 


2.08 


2.40 


10 


-3 




36 


1 


LKB Paris 


1.48- 102 


0.020 


3 


70 


10^ 


0.26 


2.50 


- 10 


-3 




45 


1 


Harris 


2.1 


2.68 


2 


82 


lO-'' 


0.32 


0. 


13 






46 


1 


Vienna 


15 


6.98 


3 


85 


105 


0.20 


25.0 


- 10 


-3 




50 


1 


MPl-QO 


1.40-104 


599 


3 


85 


10^ 


19 


27.5 


- 10 


-6 




51 


1 


Painter 


1.23 • 10^ 


5.99 - 10^ 


1 


94 


10^ 


646 


1.58 


- 10 


-6 




52 


1 


Painter 


3.38 ■ 10^ 


8.79 - lO'* 


1 


97 


10^ 


113 


5.84 


- 10 


-6 




55 


1 


Vienna 


11 


2.70 


2 


82 


10^ 


0.22 


25.0 


- 10 


-3 


1 


237 


MPl-QO / LMU 


10 


145 


1 


94 


10^ 


4.90 


19.4 


- 10 


-3 




^7 


1 


Cornell 


9.35 ■ 10^ 




2 


00 


10^ 


2943 








1 


127| 


MPl-QO / LMU 


40 


660 


3 


53 


10^ 




8.82 


10 


-3 




|691 


MPl-QO 


1.20- 10"* 


1.08- 10^ 


3 


87 


10^ 


15 


32.3 


- 10 


-6 


1 


2551 


MPl-QO 




2.68 - 10^ 


1 


93 


10-5 


2.0- lO'* 








1 


25a 


Cornell 


6.00 - 10"* 


7.32- 10^ 


1 


89 


105 


1.0- 10"* 


3.14 


10 


-6 


1 


2301 


UCSB 




6.20 - 10^ 




6.17 


9.36 








[|2 
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The most popular and sensitive technique in optomechanics involve optical cavities. There are several reviews about 
this topic, see e.g. Refs. [257, 221, 258| , 222| and we will only discuss the basic features of these devices. In a cavity, 
photons bounce back and forth many times before they are absorbed or escape the cavity. During this time, the light 
interacts with the resonator and the total interaction strength depends on the number of times that the photon reflects 
inside the cavity; the so-called finesse T ■ The light that eventually comes out of the cavity contains the displacement 
signal, enabling position detection. Also, when the light interacts w ith the resonator, the photons exert a backaction force 
on the resonator; the radiation pressure that was introduced in Sec. 3.5.3 Diff'erent types of optical cavities exist that we 
will treat separately in the next Subsections. Cavities cannot only be realized in the optical domain where the frequencies 
are hundreds of THz, but they can also be made with superconducting resonant circuits that operate in the microwave 
frequency range (GHz). At the end of this section we will introduce these microwave cavities and compare them to their 
optical equivalents. 

An important figure of merit for cavity-optomechanical systems is the coupling strength, which is usually defined as 
KoM - ducldu. It indicates how much the cavity frequency shifts per unit motion. Table |^ shows the coupUng strength, 
the cavity frequency, its damping rate, and some derived quantities for recent experiments with optical and microwave 
cavities: The optomechanical coupling length is defined as Lqm - <^clgoM, which equals the physical cavity length in 
the case of a Fabry-Perot cavity. Lqm is a convenient quantity as the displacement normalized by it, directly gives the 
relative change in cavity frequency: ujLoM - {<^c{u) - (jl>c{Q))I(jl>c{Q)- To obtain a high optomechanical coupling it is thus 
advantageous to have a small cavity. In the experiments reported in Refs. |35 , 3g 



[45[ the Fabry-Perot cavity length is 

only a few millimeters. The coupling rate that appears in the cavity-resonator Hamiltonian (Eq. ^0|) has been introduced 
as Gom - gomuo, which indicates how much the cavity moves due to the zero-point motion of the resonator and this 
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Figure 27: Different implementations of optical cavities for position detection, (a) Schematic overview of an optical cavity with a movable membrane 
in it. The cavity consists of two fixed miiTors at a distance L. The membrane at position x can be located at difi"erent intensities of the standing wave 
in the cavity (green). The membrane (inset) is made of 50 nm thick silicon nitride and is held by a silicon substrate. The membrane reflectivity is 
= 0.42. Reprinted by permission from Macmillan Publishers Ltd: Nature 452, 72-75, copyright 2008. (b) Top: scanning electron micrograph of 
a microtoroid. Bottom: schematic illustration of the optical whispering gallery mode (WGM) and the mechanical radial breathing mode (RBM). The 
WGM encircles the toroid a nd is mainl y Inr ateri at the rim. In RBM vibrations the rim of the toroid moves outward, thereby elongating the WGM cavity 
length. Adapted from Refs. 1 261 1 and [ 122|. Reprinted by permission from Macmillan Publishers Ltd: Nature 421, 925-928, copyright 2003. (c) Artist 



impression of a double-disk cavity. The cavity consists of two silica disks separated by a layer of amorphous silicon. Similar to the microtoroid in panel 
(c), WMGs exist in each of the two disks. They from bonded and antibonded modes (modes with even and odd parity respectively). The mechanical 
mode of interest is that where the gap between the two disks changes. Reprinted figure with permission from Q. Lin et al, Phys. Rev. Lett. 103 
(2009) 103601. Copyright 2009 by the American Physical Society, (d) The setup used to study the mechanics and the optical properties of the zipper 
cavity. Near-infrared (NIR) laser light goes through an erbium-doped fiber amphfier (EDFA) and a variable optical attenuator (VOA). The polarization 
is adjusted in the fiber-polarization controller (FPC) and the hght is spht in two paths: one through a fiber Mach-Zehnder interferometer (MZI) and 
the other through a tapered fiber that couples the hght to the zipper cavity. Finally, the photons in both paths are detected with photodetectors (PD). 
Reprinted by permission from Macmillan Publishers Ltd: Nature 459, 550-555, copyright 2009. 



coupling rate is sometimes called the vacuum coupling rate |259|. The role of this quantity becomes clear when one 
considers a quantum state in which both the mechanical resonator and the cavity are in an eigenstate of the individual 
subsystems. In the uncoupled system (i.e., for Gqm = 0) the system remains in this state since it is also an eigenstate of 
the total Hamiltonian. However, this is no longer the case in the coupled system. The state will then oscillate between 
states with more energy in the mechanical resonator and states with more photons in the cavity. The rate at which this 
happens is proportional to Gqm- Note that this is very similar to Rabi oscillations in a two-level system coupled to a cavity 
[E6Q]. Furthermore, the Hamiltonian in Eq. shows that the optomechanical coupling is related to the backaction: the 



force exerted by a single photon in the cavity is hgoM, which can be as large as 13 fN for the on-chip cavities that we will 



encounter in Sec. 4.1.4. Finally, we note that the coupling strengths goM and Gqm are solely determined by the device 



geometry and are independent of the optical and mechanical Q-factors. 
4.1.1. Fabry-Perot cavities 

A schematic picture of a typical Fabry-Perot setup has already been shown in Fig. and was discussed in Sec. 3.5.3. 
Briefly, the cavity consists of two mirrors, one of which can move, whereas the other one is fixed. When the moveable 
mirror is displaced, the cavity length changes and therefore also the cavity resonance wavelength. An important property 
of the cavity is the reflectivity of the mirrors. A large reflectivity means that the optical quality factor^ is high, i.e., the 
light bounces back and forth many times so that the interaction between mechanical motion and the photon is strongly 
enhanced. A small fraction of the light is allowed to come out of the cavity and is guided to a photo-detector. These can 
detect individual photons and therefore the position detection sensitivity can be shot noise limited. 

There are two main implementations to measure the light that comes out of the cavity. The first method, shown in Fig. 
Esk, uses a polarized beam splitter in combination with a /l/4-plate to separate the incoming laser light from the light that 



^^The quality factor indicate the number of oscillations of the optical field before it leaves the cavity or before it is dissipated. In optics, one usually 
speaks about the finesse 7^ of the cavity, which is the number of reflections of the field before it is lost. The finesse is related to the optical quality factor: 
'F = Acoc/k = Qx Aojc/ciJc- Here, Aixi^ is the free spectral range of th e cav ity i.e. the distance between subsequent cavity resonances {AoJc = nc/L). 
Typical numbers for range from 10^ to 10^ for high-quality cavities b23|. 
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comes out of the cavity. The advantage of this setup is its compactness. The second setup consists of two arms with cavities 
that are illuminated by a single laser via a beam splitter. Light that comes out of the two arms interferes destructively 
at the beam splitter when the length of the two arms is identical. However when one of the mirrors displaces, this is no 
longer the case and photons arrive at the photo-detector. This implementation is a so-called Michelson interferometer and 
is used to search for gravitational waves [262|. With two 4 km arms, displacements of kg-scale resonators (i.e., mirrors) 
are detectable to within 10 m/ VHz [|4]. 



For a Fabry-Perot cavity the coupling constant is goM = (^c/L, where L is the cavity length [122| as explained in 
Sec. 4. 1 For a cm-long cavity and light with a visible wavelength this gives goM ~ ^tt ■ 10 MHz/nm (see Table In the 
presence of a strong coherent pump laser, the effective coupling rate is enhanced by the square root of the number photons 
in the cavity -\/ri^. The optomechanical strong coupling regime where Gqm ■ n^J^ » lor, Kc was first reached in a 25 mm 
long Fabry-Perot cavity [237|, as evidenced by the mode splitting [263| of the mechanical resonator and the detuning of 
the cavity. Strong coupling has now also been demonstrated in other optomechanical systems [ 254 , 255 , 264 ] . 

Fabry-Perot cavities are very sensitive position detectors and have already been used for a while to approach the quan- 
tum limit on position detection, see for example Ref. [209|. Their properties (shot-noise limited sensitivity, backaction 
etc.) are therefore well known. The most sensitive implementation of the Fabry-Perot interferometer reaches the impres- 
sive sensitivity of 10"^" m/ VHz for the motion of the differential motion of the two mirrors forming a 0.25 mm long 
cavity with a finesse of 230 000 [E^]. Currently, a lot of effort is put in reaching true strong coupling where Gqm > k, ojr 



and in the observation of quantum backaction [265, 266| 



4.1.2. Movable membrane inside the cavity 

Another use of cavities in optomechanics, is the setup pioneered by the group of Harris, where a flexible membrane is 
positioned inside a rigid cavity [Q. A schematic drawing of this setup is shown in Fig. ^7| a. One of the advantages of this 
implementation is that the mechanical resonator and the mirrors of the cavity are separated, so that a high-quality cavity 
can be made without degrading the mechanical properties of the resonator The membrane, with a low (field) reflectivity 
Kc, can be positioned at different locations x inside that cavity, i.e., at nodes or at anti-nodes of the standing light-field 
waves. Depending on where the membrane is placed, the coupling is different: At an anti-node, the membrane strongly 
interacts with the cavity, whereas at a node, it does not. In a way this system can be viewed as two coupled cavities. 
For a particular optical mode, the light is predominantly on one side of the membrane exerting a radiation pressure on 
it from that side. From the above it is clear that the cavity frequency tOc should be a periodic function of the membrane 
position x„, and so is the optomechanical coupling. Using a one-dimensional model the coupling strength is obtained: 
gOM ~ 4;r|rc|c/(L/l) sin(4;7rjic,„//i) assuming Ird <K 1. This optomechanical system thus exhibits a tuneable coupling. 

At a node of the optical field (x,„ - nA/2), the linear coupling term between the cavity and the resonator vanishes and 
a small quadratic term remains. Interestingly, this quadratic term can be greatly enhanced by shghtly tilting (~ 1 mrad) the 
membrane [267|. In this case, avoided crossings of different optical cavity modes occur and near these avoided crossings 
the cavity frequency depe nds s trongly on the position of the membrane squared (i.e., square-law position detection). 
In the realization of Ref. [267| the tilting results in an increase in the coupling from 30kHz/nm^ to > 30 MHz/nm^, 
which might enable dir ect measurements of the quantization of the membrane's energy (see the discussion in Sec. |]and 
Appendix g) l|6||6|]. 



4.1.3. Optical cavities on the micro scale 

As we have discussed before, it is advantageous to have a cavity with a high finesse (or the related optic al qua lity 
factor, see footnote |3^) in combination with a high mechanical Q-factor. In the experiments described in Sec. 4.1.2 this 
was done by physically separating the mirrors from the mechanical resonator. A different approach is to use microtoroids 
(see Fig. ^7|b), which can have optical Q-factors in excess of 10^ or equivalently T > 10^ [261]. Light can be coupled 
into these devices via free-space evanescent coupling by positioning a tapered fiber close (~ 1 fj.m) to it | 270 |. The light 
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travels around the outer edge of the toroid in, what is called, a whispering-gallery mode (WGM) (see Fig. 
left panel). The light in this mode is strongly coupled to the mechanical vibrations of the toroid, in particular to its radial 
breathing mode (RBM). In this mode, the toroid expands and retracts in the radial direction, thereby changing its diameter 
slightly in time. The mechanical RBM frequencies are of the order of 10 to 100 MHz; the mechanical quality factor can 
be as high as 32000 and depends on the exact device geometry [ 27 1 ) . The coupling length is in this case the toroid radius, 
LoM - R [ 122], and the couphng constant goM - dJc/R can reach In ■ 10 GHz/nm (see Table |^). This high value can be 
understood from the much smaller dimensions of the toroid cavity: The cm-long Fabry-Perot cavities are now replaced a 
by toroid with a circumference of a few hundred jum. 



A very similar system are silica microspheres [270, 249, 63 1, where again the whispering-gallery mode couples 
strongly to the mechanical breathing modes with a frequency of the order of 100 MHz for a typical diameter of ~ 30 jim. 
The optical quality factor of these devices exceeds 10^ and typical values for the mechanical Q-factor are 10"*. Light is 
coupled into the WGM by focussing a laser beam close to the sphere. Similar to the fiber taper, the hght enters the sphere 
via the evanescent field and this causes a detectable phase shift in the transmitted light [^. 
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Cooling has been extensively studied in the microtoroid syste ms [|37[ , [1% p8| ] and an advantage is that they can be 
integrated and precooled to helium temperatures in cryostats |63|]. In the latter experiments, thermal occupation 
numbers as low as n = 70 and 37 have been achieved respectively. The small scale of these devices also bears another 
advantage in that it can be coupled to other mechanical resonator by placing them in the near vicinity of the toroid. 
An example is described in Refs. [ |67| , in which the flexural modes of a nanomechanical SiN string are probed via 
the toroid with a position sensitivity that is two times below the standard quantum limit. The coupling goM decreases 
exponentially with the distance between the string and the toroid and a maximum value of goM - InY. lOMHz/nm 
has been reported. Finally, it was demonstrated that microtoroid can be actuated electrostatically using gradient forces, 
making it an optoelectromechanical system [ 272 1 . A recent review summarizes the achievements in cavity optomechanics 
with whispering-gallery modes [ 122|. 
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4.1.4. On-chip optical cavities 

A clever way to further increase the optomechanical coupling has been realized by Painter and co-workers 
The basic idea is as follows: When placing two optical waveguides in close vicinity (submicron scale), symmetric 
and anti-symmetric optical modes form with a mode volume of the order of . The optical coupling length, Lqm, can be 
viewed as the length scale over which a photon's momentum is transferred, which, in this case, is reduced to a length scale 
of the order of the optical wavelength, A. Thus, the coupling can be estimated to be goM ~ Wc//l. Since the wavelength 
is about 1 .5 fim, the coupling can be at least an order of magnitude larg er than in the micro-cavities discussed above (see 
Table ||). Two implementations have been built. In one version [127, ^5[|, a pair of silica (or SiN) disks separated by 
nanometre-scale gaps was used as shown in Fig. |27| c. The coupling depends on the separation between the disks and 
for an air-gap of 138 nm, the coupling was found to be goM = 27r ■ 33GHz/nm. Efficient cooling of the mechanical 
mode was achieved |^], and static and dynamic mechanical wavelength routing wa s dem onstrated [ 231 1. Furthermore, 
Wiederhecker et al. demonstrated attrac tive and repulsive forces between the disks [ 127 | and optomechanical tuning of 
the cavity modes over more than 30 nm [23C[ in a similar device. 

An even larger coupling is obtained with two stoichiometric silicon nitride ladder structures with a photonic crystal 
structure ("zipper cavities") as illustrated in Fig. pTp. With a separation of 120 nm between the two waveguides, a 
coupling of goM — 27r ■ 123 GHz has been achieved [|62[|. The strong coupling yields a large optical spring effect (see 
Sec. 3.5.3), where the resonance frequency is mainly determined by the laser field instead of by the structural properties 
of the resonator This effect shows up as a change in resonance frequency if the input power and detuning are changed. 
Interestingly, the optical spring only acts on the differential motion of the beams. The common-mode vibrations are 
not affected by the light field. By varying the detuning of the driving light, Lin et al. could shift the differential mode 
through the resonance frequency of the common mode. The observed Fano-like lineshape indicates coherent mixing of 
the mechanical excitations [ |196[ | . 

Another advantage of optomechanical crystals is the control over the location of the optical and mechanical modes 
Since the spatial extent of modes differs, the coupling strength can be engineered. Thus, the simultaneous con- 
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finement of optic al (photonic crystals) and mechanical modes (phononic crystal) leads to strong, controllable light-matter 
interactions [273|. Finally, n ote th at in on-chip o ptomechanical devices the optical gradient force is typically much larger 
than the photon pressure [^, 274,^, 127, 275 1. 



4.1.5. Superconducting microwave cavities 

Instead of using an optical cavity, electromagnetic waves can also be confined in a superconducting microwave cavity 
[ [276( 1 . The electromagnetic cavity mode is sensitive to capacitive changes and this principle has been used to measure 
nanomechanical motion [ |50| , ^ ^]. The device used in the first realization of such a transmission-line position detector 
[ ^0[ [ is shown in Fig. 

The working principle is as follows: In the superconducting transmission line with an inductance and capacitance per 
unit length, radiation propagates. On one side microwaves are injected via a capacitively coupled feedline; the other end 
can be open or shorted to ground (so-called Ajl and AjA microwave cavities respectively). In the first case standing waves 
with a voltage anti-node at the end form, whereas in the latter case waves with a node at the end form. The mechanical 
resonator is connected to the stripline and positioned close to the ground plane. Its displacement couples to the cavity 
capacitance and this in turn changes the cavity resonance frequency which can accurately be probed via transmission or 
reflectivity measurements. In the realization of the group of Lehnert (i.e., the one shown in Fig. this is done via the 
transmission of the feedline [50 223, 5^ in the Schwab-group implementation the transmission through the cavity is 
measured [^, 252 1. 

Typical lengths of these on-chip resonators are in the millimetre range so that the microwave resonance frequency is 
a few GHz. The electric quality factor is typically 10"^; in the experiments, the dissipation is determined by the coupling 
to the feedline and not by internal losses (i.e., they are overcoupled) [|50|]. An overview of the key parameters of these 
devices is given in Table ^. Large optomechanical coupling is obtained for thin, long doubly clamped beams that are 
positioned as close as possible near the ground plane. A coupling of goM - 27r ■ 84kHz/nm is reported [^3| for a 170 
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Figure 28: Illustration of the stripline resonators used by Regal et al. (a) Chip with 6 meandering stripline resonators (pink) and the straight feedhne 
(green). The AjA resonators have a slightly different length, which allows for frequency multiplexing, (b) The resonators are coupled to the feedline by 
a capacitive elbow coupler, (c) and (d) zooms of the resonator. The 50 /jm long aluminum beam resonator is under tensile stress which is induced by 
annealing at an elevated temperature. A movement of the beam changes the capacitance of the stripline to the ground plane (blue) and thereby the cavity 
resonance frequency. Reprinted by permission from Macmillan Publishers Ltd: Nature Physics 4, 555-560, copyright 2008. 



nm wide and 140 nm thick beam, that is formed from 60 nm of stoichiometric, high-stress siHcon nitride and 80 nm of 
aluminum. The resonator is located 75 nm from the gate electrode. Further substantial improvements in the coupling 



strength are difficult with this particular geometry and fabrication technology. Sulkko et al. [253| used a focussed ion 
beam to create a very small gap of about 10 nm between a mechanical resonator and a resonant microwave circuit and 
obtained goM -2n ■ \ MHz/nm, thus increasing the coupling considerably (see Table [s]). 

An important difference with optical cavities is that there are no single photon detectors available for microwave fre- 
quencies. In optical systems these do exist, and the detection of the light in itself is quantum-limited. Present commercial 
microwave amplifiers, however, always add substantially more noise than required by quantum mechanics since they are 
not shot-noise limited. This problem can largely be overcome by using a Josephson parametric amplifier Teufel et al. 
( ]66| ] have made a nearly shot-noise-limited microwave interferometer and demonstrated nanomechanical motion detec- 
tion with an imprecision below the standard quantum limit. On the other hand, similar to their optical equivalents, the 



photon-pressure backaction of the stripline position-detecto rs on the mechanical resonator has been observed [q6[| and 



even backaction-evading measurements have been reported [|252||. In the latter experiment, a single-quadrature measure 



ment of motion with a sensitivity of four times the zero-point motion has been demonstrated. Sideband cooling has been 



performed in a series of experiments [q6[ ^ [252| ]. Very recently Teufel et al. were the first ones to demonstrate ground 
state cooling using a superconducting LC resonator [^], in which a movable membrane both acts as the capacitor of the 
LC circuit and as the mechanical oscillator ("drum resonator" geometry). They obtained a thermal occupation of n = 0.34, 
where the drum resonator is in the ground state for 57% of the time. In the same type of device they also observed strong 



coupling between the mechanical resonator and the LC resonator |254|. 

It should be noted that mechanical resonators have also been coupled to non-superconducting resonant circuits [ p4| . 



277 1 . Sideband cooling from room temperature to 45 K and 100 K respectively has been demonstrated in these systems. 
The first method differs from the abovementioned superconducting cavities in that it consists of a stripline resonator with 
lower frequency (100 MHz) and lower quality factor (234). In the second experiment, the resonator is an open-ended 
coaxial cable with a resonance frequency of 1 1 GHz and a 2 of 80. 

4.2. Transmission modulation based detection methods 

Mechanical motion can also be detected by measuring the transmission of electrons or photons through devices which 
embed a movable part. The motion of the mechanical resonator modulates the electron or photon transmission and the 
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Figure 29: (a) Readout using the light transmitted from one waveguide (left) to another one (right). When one of the cantilevers moves, the transmission 
of the light changes. This is a non-interferometric optical method to detect mechanical motion. Reprinted by permission from Macmillan Publishers 
Ltd: Nature Nanotechnology 4, 377-382, copyright 2009. (b) An atomic-point-contact displacement detector. A tunnel contact between the beam and a 
nearby electrode is made using electromigration. When the Au beam is far from the electrode, the tunneling probabihty for electrons (see inset) is low 
and the resistance is high. Closer to the electrode, the transmission increases and the resistance is lower. Reprinted figure with permission from N. E. 
Flowers-Jacobs, D. R. Schmidt, K. W. Lehnert, Intrinsic noise properties of atomic point contact displacement detectors, Phys. Rev. Lett. 98 (9) (2007) 
096804. Copyright 2007 by the American Physical Society, (c) Setup to measure multi-walled nanotube vibrations via field emission. Electrons are 
emitted by applying a voltage Va = -500 to -900 V between the nanotube and a nearby (3 cm) screen. The electrons hit the screen at a certain position 
and show up a bright spot. The insets at the bottom show the spots for a silicon-carbide resonator tha t is d riven off-resonance (left) and on resonance 
(right). In the latter case the spot is blurred due to the motion of the resonator. Adapted from Refs. |278| and [po|]. Reprinted with permission from 



S. T. Purcell et al. Tuning of nanotube mechanical resonances by electric field pulling, Phys. Rev. Lett. 89 (27) (2002) 276103. Copyright 2002 by the 
American Physical Society and S. Perisanu et al. Appl. Phys. Lett. 90 (4) (2007) 043113. Copyright 2007, American Institute of Physics. 



output signal contains spectral information about the resonator dynamics. Below we will discuss in detail two basic 
configurations of this detection scheme: optical waveguides and electron tunneling devices, including electron shuttles. 

4.2.1. Optical waveguides 

As mentioned above, light can be guided on a chip using optical waveguides. These are pieces of transparent material 
(e.g. Si, Si02 or SiN) which are carved out of the underlying substrate. Light can propagate through these waveguides 
and can be directed to any place on the chip. Just outside the waveguide the electromagnetic field is not zero and decays 
exponentially with the distance from the waveguide. When two waveguides are placed close together, light can be trans- 
mitted from one to the otheiQ. The transmission between the waveguides depends on the distance between them and they 



can be employed as a position detector. The method is also compatible with operation in water |280|, which is important 



for (bio)sensor applications, and it has also been used to make optoelectronic switches [281]. Li et al. | p9| ] have used 
the transmission-modulation principle to measure the mechanical motion of cantilever structures made in a silicon-on- 
insulator platform with high sensitivity. In particular, they placed the suspended ends of two waveguides facing each other 
at a distance of 200 nm as shown in Fig. p9| a. The transmission of photons through the gap depends on the distance and 
the exact position of the cantilevers (misalignment). The thermal motion of the cantilevers modulates the transmission 
and the transmitted light is fed into a photodiode. The spectrum clearly reveals multiple resonances of the cantilevers with 
a sensitivity of 40 fm/ VHz at room temperature. An advantage of this method over optical cavities is that the detection 
scheme allows for transduction of nanomechanical motion over a wide range of optical frequencies, instead of working 
only at certain well-defined wavelengths. 

4.2.2. Electron tunneling 

When two metallic electrodes are placed close (up to a few nanometer) to each other, electrons can tunnel through 
the gap between them. With a voltage applied between both electrodes, a net current flows that depends exponentially 
on their distance. This exponential distance dependence can be used to measure displacements of mechanical resonators. 
Flowers-Jacobs et al. [S2l| have measured the displacement of the doubly-clamped gold beam shown in Fig. 2% using 



an atomic point contact (APC) made by electromigration [282|. Since the tunneling resistance between the APC and the 



beam is high, a tank circuit is used for impedance matching to 50 Q high-frequency amplifiers. The Brownian motion 
of the beam is observed at a temperature of 250 mK with a shot-noise limited imprecision of 2.3 fm/ VHz. From the 



^''A related, but different method is demonstrated in Ref. |279|, where photons are transferred from one waveguide to another, which runs parallel to 
the former. The coupling between the waveguides depends on the distance between them (i.e., on the displacement) and photons that enter the second 



waveguide are "lost" and not detected at end of the first waveguide. 
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temperature dependence of the thermal noise spectra, a backaction temperature of 0.7 K is found, which is larger than 
expected from the momentum carried by the tunneling electrons. Kan et al. used a similar technique where an STM tip 
was positioned above a MEMS resonator [283]. The current modulation due to the resonator motion was down-mixed 
(Sec. 4-. 3. 3) to overcome bandwidth limitations. Based on these methods an interesting new type of detector is proposed 
[284|. By incorporating the movable tunnel junction into a loop threaded by a magnetic flux, the APC can be used to 



measure the momentum instead of the position of the resonator. Also, other more sophisticated schemes are envisioned 
for detecting entanglement in the mechanical quantum oscillator [285]. 

Another method is based on field emission of electrons from a vibrating tip. Consider, for instance, a multi-walled 
carbon nanotube [ 278 ] mounted under vacuum in a field emission setup (see Fig. ^9|c). A large voltage is applied between 
the carbon nanotube and an observation screen which lights up at the position where electrons are impinging. Electrons 
are accelerated from the tip of the nanotube to the screen by the electric field. Motion of the nanotube cantilever results 
in a small blurring of the spot at the screen. The thermal motion of the nanotube already blurs the spot a bit, but when 
the vibration amplitude is enlarged by applying a RF driving signal on nearby electrodes the blurring of the spot becomes 
more pronounced on resonance. The electric field also pulls on the nanotube, thereby increasing the resonance frequency 
[278]. The method has been used to build a nanotube radio 174], a mass sensor with atomic resolution [[ttI, and high Q 
silicon carbide nanowire resonators at room temperature [Q. 

In more complicated circuits shuttling of electrons has been considered [286, 125| , 287]. These devices operate as 
follows: there are two electrodes with a movable metallic island in between. When the resonator is closer to one of the 
electrodes, the transmission between the island and this electrode increases, allowing for electrons to jump on the island. 
The electrostatic force drives the island with the negatively-charged electrons to the other electrode. This way, they are 
transported to the other electrode by the resonator (mechanical transport of electrons). In steady state this shuttling occurs 
at the resonance frequency of the oscillating island. Interestingly, the voltage applied between the two electrodes can 
amplify this motion and lead to large amplitude oscillations; an example of an electomechanical instability^^ This subject 
has been extensively reviewed by Gorelik et al. ] 288 1 and we refer to this work for further reading. 



4.3. Capacitive detectors 

In many nano-electromechanical devices, movable metallic parts form capacitances with nearby metallic electrodes 
or ground planes. Displacement of the mechanical structure inevitable leads to a change in these capacitances which 
can be detected electrically if they are large enough. Likewise, the motion of the mechanical resonator can be actuated 



by applying voltages between the dififerent electrodes as we have discussed in Sec. 2.4.4 This capacitive actuation and 



read-out of mechanical motion has been appUed in several ways. Here, we will discuss three popular methods. 



4.3.1. Single-electron transistors 

One of the first mesoscopic devices that was used for position detection was the single-electron transistor (SET) [p9|. 
This device consists of a metallic island that is connected by tunnel contacts to a source and drain electrode. The total 
capacitance of the island is so small that the energy required to add one electron surpasses the thermal energy [289[ and 
electrons can only enter and leave the island one by one (sequential tunneling). The resulting single-electron current 
through the island is very sensitive to the electrostatic environment. Thus, if a nearby resonator is capacitively coupled to 
the island, and it has a diff'erent electrostatic potential, the device can be used as a displacement sensor, as illustrated in Fig. 
|o[ The first realization of such a device was made by Knobel and Cleland [^9||, where the SET measured the position of 
a 1 16 MHz beam resonator with a sensitivity of 2 fm/ VHz. This value corresponds to a position resolution a factor of 100 
above the standard quantum l imit, as is indicated in Table 0. The SET was operated as a mixer, a technique that we will 
discuss in more detail in Sec. 4.3.3 . Subsequent work by LaHaye et al. [ |30| ] demonstrated, using a superconducting radio- 
frequency SET, a sensitivity of four times the standard quantum limit. The role of the RF SET is to provide impedance 
matching for the high-frequency resonator signal (19.7 MHz). 

Using an improved version of this device (Fig. ^0|b), Naik et al. [|2|] observed the backaction of the superconducting 
SET on the mechanical resonator; a change in the capacitance leads to a change in the average charge ("occupation") 
[ ]290| ], which in turn leads to a change in the electrostatic force the resonator experiences. Depending on the bias conditions 
(gate voltage) of the SET, shifts in the resonance frequency and the damping rate were observed. Thus the backaction 
force on the detector, characterized by the function Af that was introduced in Sec. 3.4.1, is bias dependent for this 



type of detector 1240, 212, 291, 292]. The effective resonator temperature has been measured for different coupling 



strengths between the resonator and the SET. In the experiment this is done by varying the voltage difference between 
the resonator and the SET. For small couplings, the resonator temperature equals that of the sample chip i.e., close to 



This shuttling instability can thus be viewed as an example of backaction induced self-oscillations as discussed in the previous Section. When 
pursuing the quantum Hmit in these devices, people refer to the ability to transport exactly one charge earner per cycle from source to drain; the 
zero-point motion does not play a crucial role in achieving this. 
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Figure 30: (a) Illustration of the dependence of the current through a SET on the electrostatic potential. The cun'ent is almost zero in the Coulomb 
blockade regime and displays so-called Coulomb peaks where the blockade is absent. The resonator motion modulates the potential (horizontal arrow) 
and this leads to a proportional change in the current, (b) The device used in the Schwab group. The superconducting single-electron transistor (SSET) 
is made from aluminum and is connected to a source and drain electrode. The potential of the SET can be varied using a gate electrode and is also 
influenced by the position of the beam resonator. From M. D. LaHaye et al. Science 304 (2004) 14-71. Reprinted with permission from AAAS. (c) 
Illustration of self-detection via rectification. Motion of the resonator leads to a decrease in the average current on top of the Coulomb peak, whereas 
deep in the blockade region it leads to a (small) increase, (d) SEM picture of a ultra-high quality factor nanotube resonator. From G. A. Steele et al. 
Science 325 (2009) 1103-1107. Reprinted with permission from AAAS. 



the base temperature of the dilution refrigerator. For higher coupUngs, the resonator temperature is set by the backaction 
temperature of the SET, and is raised to about 200 mK. 

It should also be possible to use a SET to sense the vibrations of a suspended c arbo n nanotube, but the small capacitive 
coupling between the nanotube and the SET island makes this a challenging task [ 293]. It is more advantageous to use the 
suspended carbon nanotube itself as a self-detecting SET. Using current rectificati on a nd frequency mixing, information 
about the driven motion of the suspended nanotube has been obtained [ |8l| . 



1061 



In these experiments, a strong 
coupling between mechanical motion and the charge on the nanotube has been observed. The strong coupling is due 
to the electrostatic force generated by individual electrons tunneling onto the nanotube. For the devices used in these 
experiments, the change in equilibrium position of the nanotube after adding a single electron easily surpasses the zero- 
point motion. Typically, the single-electron tunnel rate is much larger than the resonance frequency, indicating that the 
backaction is determined by the average number of electrons ("occupation") on the nanotube. This backaction leads 
to frequency shifts and changes in damping as a function of gate voltage. To be more specific, the damping increases 
dramatically with the amount of current flowing through as the electron tunneling produces a large stochastic backaction 
force. 

In Refs. [ plj , |106[ ] readout using current rectification^ is employed. While the nanotube motion is actuated by a RF 
signal on a nearby antenna, the detected signal is at DC. The key to understand this is the notion that nanotube motion 
effectively translates into an oscillating gate voltage, which smears out the sharp features of the SET current, as illustrated 
in Fig. [295]. The technique is of special interest as it constitutes a square-law detector (see Appendix Moreover, 
it allows for the motion detection with small currents, enabling the observation of ultra high Q-factors, exceeding 100,000 
at millikelvin temperatures. The low dissipation enables the observation of single-electron tuning and frequency tuning 
oscillations, analogous to the Coulomb oscillations in th e SET current. Recently, this self-detecting rectification scheme 
has also been employed for a thin aluminum beam 1 296 ] . 



4.3.2. Quantum-point contacts 

A quantum-point contact (QPC) is a na rrow constriction in a two-dimensional electron gas, whose conductance can be 
adjusted using electrostatic potentials [ |297| ]: Every time another channel for electrons becomes available, the conductance 
increases by one conductance quantum 2e^/h. In practice these sharp steps are smoothed by temperature and when the 
QPC is biased near such a step it is very sensitive to changes in the electrostatic potential. If mechanical vibrations 
modulate the electrostatic fields, the QPC can also be used as a position detector Cleland et al. have used this principle to 



The rectific ation and mixing measuring method have in common that they both rely on the nonlinearity of the cun'ent- voltage or cun'ent-gate-voltage 



characteristics 1 294 1 ; the experimental implementation is, however, different. It sh 



d that rectification is more commonly used to detect 
.21) the dc tunnel current contains information about the 



displacement in mechanical resonators. For example, in the experiments of Ref. 
vibrating carbon nanotubes since a time-dependent displacement changes the source-drain distance and thereby the averaged current, 
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measure the vibrations of a beam made out of a single-crystal GaAs hetero structure with a sensitivity of 3 pm/ VHz [298|. 
In this case the QPC is part of the beam and frequency mixing (Sec. 4-.3.3) is used to detect the displacement through the 
voltage across the QPC. The signal is amplified through the piezoelectric effect in GaAs that will be discussed in more 
detail in Sec. ^ 

A similar sensitivity has been achieved in a completely different setup by Poggio and coworkers [^. Now a 5 kHz 
micromechanical resonator is hanging above a QPC that is located on a different substrate. Thermal noise specta have 
been recorded and the authors have mapped out the transduction factor as a function of cantilever position relative to the 
QPC. Furthermore, they observe that the cantilever Q-factor is not affected by the QPC source-drain current, indicating 
weak coupling and therefore negligible backaction. 

A QPC can also be used to probe vibrational modes of the host crystal itself [ |299| ]. This substrate is a truly three- 
dimensional resonator that consists of on the order of 10^*^ atoms. Strictly speaking, the latter example does not constitute 
of a capacitive detector as the transduction is through the piezoelectric effect; the motion directly influences the source- 
drain voltage and not the conductance of the QPC. 



4.3.3. Frequency mixing 

Frequency mixing has been adop ted in to\ 
motion into a low-frequency signal [ 298 , 29^ 



-dow n solid-state devices as a versatile technique to convert high-frequency 
300 1 . It is most often used in combination with capacitive detection tech- 



niques, although it is also used in combination with electron tunneling (Sec. |4.2[ ) and piezoelectric resonators (Sec. |4.4| ). 
The basic principle is as follows. Consider the generic relation between the input (i.e., the displacement u) and the detec- 
tor output V introduced in Sec. v(cj) = Ay(cj)u{a)). By modulating the transduction /I,, at a frequency fio (often called 
the local oscillator frequency), the displacement at frequency is converted into a signal with frequencies -i- fio and 
/r - fio- The latter component at the difference frequency can be chosen to be at a frequency far below the resonator 
frequency: for example in the kHz range. The signal at this frequency is not affected by the RC time of the measurement 
setup and can therefore be measured straightforwardly. 

The technique has become of particular interest for detecting vibrational motion of bottom-up devices. Sazonova et 
al. were the first to apply frequency mixing to suspended carbon nanotube resonators. They observed multiple gate- 
tunable resonances with Q-factors on the order of 100 at room temperature. Subsequently Witkamp et al. [ pO| ] identified 
the bending mode vibrations of a carbon nanotube. Nowadays the technique has been employed by many groups, not 
only restricted to carbon nanotubes [ [76[ [T^, 126 8^|8^, |5||, but also applied to suspended graphene sheets [ 103[ 104], 
free-hanging semiconducting nanowires [p4[ [96^ and charge-density- wave sheets [105]. Furthermore, several variations 
to the original mixing scheme have been implemented, including FM [pSh and AM modulation [126]. 



A more detailed understanding of this self-detecting method can be obtained by considering, for example, the sus- 
pended carbon nanotube of Fig. The conductance of a semiconducting nanotube depends on the induced charge, 
which is the product of the gate voltage and the gate capacitance, i.e. G = G{CgVg). As we have shown before, the 
gate capacitance is position dependent and therefore the conductance varies in time with the resonator frequency. In 
nanotube mixing experiments, the current is measured, which is the product of the conductance and the bias voltage: 
I = Vt dO/du ■ u, from which we identify /I,, - Vb dG/du. The modulation of the transduction can thus be done by 
applying an ac bias voltage: Vb = Vac cos(a)Lot), as explained in Fig. In these experiments backaction comes again 
from the electrons flowing through the nanotube | 



4.4. Piezoelectric resonators 

A different mechanism to transduce a displacement into an electrical signal uses piezoelectricity. In a piezoelectric 
material a stress induces an electric polarization, or equivalently, an electric displacement D [137[. This results in a 
voltage difference across the resonator that can be used to infer the displacement. Piezoelectric resonators are different 
from most other methods discussed in this Section as no separate detectors are needed; they detect their own motion 
instead. Moveover, they are also different from the other self-detecting resonators (e.g. the suspended carbon nanotubes 
and graphene resonators that were described above) since no external signal needs to be applied in order to measure 
the displacement (for example, in the nanotube experiments of Sections 4.3.1 and 4.3.3 one has to apply a source-drain 



voltage to transduce the displacement). Still, from the discussion in Sec. 3.4 it is clear that backaction forces should also 



act on a piezoelectric resonators when measuring the generated voltage. This happens through the converse piezoelectric 
effect, where an electric field generates a strain in the material. Referring back to the definition of stress and strain in Sec. 
^ the direct and converse piezoelectric effect are [137|: 



Di - dijko-jk. 



Jij - dijk&k- 



(71) 



Here d is the third-rank tensor with the piezoelectric coefficients. Similar to the properties of the elasticity tensor, many 
of the elements dijk are zero, or related to other elements, depending on the symmetries of the crystal structure of the 
material. A necessary requirement to have at least one nonzero element is that the material's unit cell does not have an 
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Figure ^1 : Schematic overview of tlie measurement electronics used by Witkamp et al. for position detection of a suspended nanotube using frequency 
mixing fccr. A radio frequency (RF) generator applies an ac voltage to the back-gate electrode to drive the suspended nanotube. A dc gate voltage is 
added via a bias-T (indicated by the "+"). The same generator is used to generate the ac bias voltage by mixing its output with the reference output 
of the lock-in amplifier. At the source electrode the voltage has spectral components at / 4- A/ and / - A/, whereas the gate voltage is oscillating at 
frequency /. The nanotube mixes both signals, which results in an output current at the drain electrode with spectral components at A/, / + A/, / - A/, 
2/ + A/ and 2/ - A/. The A/ part of the current flowing through the nanotube is converted into a voltage and is measured with a lock-in amplifier. The 
printed circuit board (PCB) with the sample, bias-T and 5012 terminator is located inside the vacuum chamber of a probe station. 



inversion center. For most piezoelectric materials, including GaAs and AIN, the piezoelectric coefficients are of the order 
of pC/N and they can have both positive and negative values. 

Piezoelectric displacement detection of small structures has first been employed with AFM cantilevers [ 301 , 302[ , 303 1, 
before extending the technique to the nanomechanical domain. Tang et al. used this detection scheme[^ to detect the 
motion of GaAs/AlGaAs beams [304]. The beams are made asymmetric to ensure that the piezoelectric signal is not 



nulled by the opposite stresses at both sides of the neutral plane (Sec. 2.3). Cleland et al. used piezoelectricity to detect 



the motion of GaAs beams with an integrated QPC in it. There, a deflection of the beam induces in-plane stress and the 
resulting out-of -plain field D acts as an effective gate voltage. This then changes the current through the QPC [298]. 



An important consideration is that piezoelectric materials cannot generate large currents. The resonator should therefore 
be connected to a high-impedance load, such as the gate electrode of the QPC [ [298| ], to the gate of a SET p05[ ], or a 
high-impedance amplifier [ 304 ] . 

Mahboob and Yamaguchi used the two-dimensional electron gas in a flexural beam resonator made of a GaAs/AlGaAs 
heterostructure to measure the displacement and they demonstrated parametric amplification by modulating the beam's 
resonance frequency using the converse piezo-electric effect [306|. The converse piezo-electric effect was also used to 
tune the frequency of a doubly-clamped zinc-oxide nanowire by Zhu and coworkers Okam oto e t al. found an 

enhancement of backaction effects due to excitation of carriers in an piezoelectric GaAs resonator | 307 |. Many of the 
abovementioned effects were combined by Masmanidis and coworkers who used a three-layer structure, consisting of 
p-doped, intrinsically doped, and n-doped GaAs to demonstrate piezoelectric actuation, frequency tuning, and nanome- 
chanical bit operations [308]. By electrically shifting the carriers with respect to the neutral plane, the strength of the 
actuation could be adjusted. Finally, by capacitively coupling a bulk AIN resonator to a qubit, O'Connell et al. obtained 
the high coupling rates required to perform quantum operations on a mechanical resonator [22]. This experiment was 



already mentioned in the Introduction and will be discussed in more detail in Sec. 5.2 



4.5. Flux-based position detectors 
4.5.1. SQUIDs 

Superconducting quantum interference devices (SQUIDs) are well known for detecting small magnetic signals, such 
as those generated by our brains. These devices consist of a superconducting ring in which one or more Josephson 
junctions are incorporated ]309]. The voltage across the loop depends on the amount of magnetic flux that threads the 



loop, 'I>,„ag- Thisflux-dependencej^j has been employed to sense the motion of a cantilever with a small magnetic particle 
attached to it [|T^. 

It is, however, a l so possib l e to detect dis p lacements by incorporating a mechanical resonator in the SQUID loop 
(^TTj p2[ 13131 |ll |l5l |3l|, ITtI km |3T|, |T|, |2§; In the presence of a constant magnetic field, a change in the 



In this experiment the measured signal contained both piezoelectric and piezoresistive components. 
'^SQUIDs are also used indirectly for position detection, where the SQUID is used as low-noise voltage amplifier that amplifies the signal generated 
by capacitive detection [pJ]. In this section, we only focus on direct position detection with SQUIDs. 
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Figure 32: Schematic overview of a dc SQUID position detector, including tlie suspended beam resonator and measurement setup. Tlie SQUID loop 
is indicated in brown. The output voltage of the SQUID depends on the bias current /g that is sent through the SQUID and also on the amount of flux 
<D„„g through the loop. A magnetic field B transduces a beam displacement u into a change in magnetic flux and subsequently in a change in the output 
voltage V. The flux ^i„ag is fine-tuned with a stripline current If . Backaction results from the circulating current in the SQUID loop J. Reprinted figure 
with permission from M. Foot et ciL, Phys. Rev. Lett. 105 (2010) 207203. Copyright 2010 by the American Physical Society. 

resonator displacement changes the loop area and thereby the flux through the loop. This is illustrated in Fig. ^ Recently, 
Etaki et al. have used a dc SQUID as a sensitive detector of the position of an integrated mechanical resonator. They 
detected the driven and thermal motion of a 2-MHz buckled-beam resonator with femtometre resolution at millikelvin 
temperatures [ p^ and employed active feedback cooling to cool the resonator to 20 mK [[70|]. In the present experiments, 
the sensitivity is limited by the cryogenic amplifier and not yet by the, in principle quantum limited SQUID itself. 

Backaction has also been observed, leading to tunable shifts in the resonance frequency and damping of the resonator 
Different from the backaction in the capacitive readout schemes, the backaction of the SQUID has an inductive 



[242] 



character as it is caused by the Lorentz force generated by the current circulating in the loop of the SQUID. This current 
also runs through the resonator (Fig. and the magnetic field that couples displacement and flux, generates a Lorentz 
force on the resonator 



4.5.2. Magnetomotive 

A fast and relatively ea sy way to actuate and read-out nanomechanical motion in clamped-clamped resonators is the 
magnetomotive technique [321, 322, |8, 194]. Nowadays, this method is mainly used to determine the vibrational fre- 
quencies and to characterize dynamic properties such as the Q-factor and nonlinear behavior. It is, however, of limited 
use for QEMS experiments and we will therefore only briefly describe the mechanism: An ac current lac is sent through 
the conducting (part of the) beam which is placed in a strong static magnetic field. The ac current causes an ac Lorentz 
force that drives the beam fi„B(Iac. At the same time, motion of the beam induces a time-varying voltage (Faraday's law) 
at the driving frequency Vemf - MnBtii. By sweeping the driving frequency, the resonance can be found by measuring the 
frequency response of the voltage over the beam. The required large magnetic field is typically generated using supercon- 
ducting magnets in a cryogenic environment, but recently room temperature magnetomotive actuation and detection has 
been demonstrated with a 2 T permanent magnet [ 323 1 . An important feature of this technique is that observed mechanical 
modes can be distinguished from electronic resonances by varying the magnetic field. Since both the driving force and 
the detector signal are proportional to B, the mechanical signal scales as [324]. 



4.6. Level spectroscopy in suspended quantum dots 

At low temperatures, suspended quantum dots may reveal information about their vibrational states: in the parameter 
regime dominated by Coulomb blockade physics, the quantum mechanical level spectrum of the confined electronic sys- 
tem can be characterized by transport measurements. The levels show up as steps in the current-voltage characteristics. 
A more accurate measurement involves the recording of a so-called stability diagram in which the differential conduc- 
tance dl/dV plotted as function of gate voltage Vq and source-drain voltage V. Now the excitations (due to electronic 
or vibrational degrees of freedom) show up as lines, whose energy can directly be read out [ 325 1 as illustrated in Fig. 
p3} In the case of vibrational states these are also called vibrational side bands and when the harmonics of a particular 
vibration are excited, the lines form a spectrum with equidistant spacing. Thus, electron tunneling through suspended 
quantum dots can excite vibrational modes and these modes can then be detected as steps in the current-voltage character- 
istics. Their observation involves, however, two important considerations: the energy resolution and the electron-vibron 
coupling. Concerning the energy resolution of transport spectroscopy measurements one must have hcjR > ksT so that 
measurements are always performed in the quantum limit of the mechanical mode. 

Secondly, the electron-vibron coupling must be high enough; it can be characterized by a dimensionless parameter 
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Figure 33: Stability diagrams for zero (a), intermediate (b), and large (c) electron-vibron coupling g„. A stability diagram is a grayscale plot of the 
differential conductance plotted against the gate and source-drain voltage. Coulomb peaks (Fig. pOl) appear here as dark lines, whereas regions with low 
differential conductance dl/dV are white, (a) shows a regular stability diagram without signatures of the mechanical resonator. There are regions where 
the cun'ent is blocked and regions where cun'ent flows. These are separated by the high-conductance lines (dark). In (b) a finite electron-vibron coupling 
is present and mechanical excitations appear as lines parallel to the original lines of high conductance (vibrational sidebands). The source-drain voltage 
where these lines cross the Coulomb diamond edges, equals hfn /e as indicated. For large the current is blocked at low source-drain voltage as shown 
in (c) and vibrational excitations can only been seen at high bias voltage. 




gev - {(^u/uof', where Am is the shift of the resonator position^ induced by adding one elementary charge, and uo is 
the zero point motion of the mechanical oscillator |32(:]. The parameter ge^, determines the step height in the current- 
voltage characteristic. We can consider three regimes in describing the influence of vibrational modes on transport: the 
weak electron-vibron coupling regime with g^y <s: 1, the intermediate regime with 0.1 < g^,, < 1 and the strong coupling 
(gev » 1) limit. The boundaries of the intermediate regime are, however, somewhat arbitrary. Figure ^ shows calculated 
stability diagrams for three different values of the electron-vibron coupling. In the weak coupling regime (see panel 
(a)) only the regular Coulomb step is present (no side bands are visible) and consequently vibrational modes cannot be 
detected in a transport experiment. Only for sufliciently large electron-vibron coupling, one or multiple so-called Franck- 
Condon steps can be observed in the current, which show up as lines in the stability diagram (the intermediate regime, 
panel (b)). This was first demonstrated in molecular junctions e xcitin g Ceo radial breathing modes [327|, an d lat er in 
suspended carbon nanotubes probing the radial breathing mode [328], or the longitudinal, stretching modes [329|. In 
the strong electron-vibron coupling limit, the vibrational induced excitations are only seen for larger bias voltages as the 
height of the first steps is exponentially suppressed. In panel (c) this is clearly visible: the grey lines representing the 
vibrational excitations are only present for high source-drain voltages. Importantly , the supp ression holds for any gate 



voltage and as a result the current at low bias is suppressed in the whole gate range [ |330i [331| |. Degeneracy points in the 
Coulomb diamonds are no longer visible in the stability diagrams and one speaks of phonon blockade of transpor t. Pho non 
blockade has been observed in suspended quantum dots embedded in a freestanding GaAs/AlGaAs membrane [ 332 1 and 
subsequently in suspended carbon nanotube quantum dots [333|. 

Recently, is has been shown that inelastic electron tunneling spectroscopy (lETS) can also be used to gain informa- 
tion about vibrational modes in suspended quantum dots. Hiittel et al. [334| observed a harmonic excitation spectrum 
connected to the longitudinal stretching modes of a suspended carbon nanotube in the Coulomb-blockaded regime, while 
temperature only allows the observation of a single excitation. The non-equilibrium occupation of the modes is explained 
by the pumping via electronic states, revealing a subtle interplay between electronic and vibrational degrees of freedom. 



Actually, this is the shift in the location of the minimum of the potential energy of the resonator. The tunneling of the electron is assumed to take 
place on a timescale that is much faster than the resonator can react. 
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Figure 34: A mechanical resonator coupled to a Josephson qubit. The mechanical resonator (left) is cooled cryogenically to its ground state and 
manipulation of its quantum state at the single phonon level has been demonstrated. The quantum properties of the resonator are probed by the 
superconducting qubit which is positioned on the top side of the right panel. Reprinted by permission from Macmillan Publishers Ltd: Nature 464, 
697-703, copyright 2010. 



5. Prospects 

So far, we have mainly concentrated on the limits of (linear) displacement detection. As discussed before, there 
are other ways of demonstrating quantum behavior. For example, square-law detection would directly probe the energy 
eigenstates of the mechanical resonator (Sec. ^ and App. y). Another approach is to couple a mechanical quantum 



oscillator to another quantum system such as the well-studied Josephson qubit [ 335 1: the state of the mechanical resonator 
changes the state of the qubit, which can then be probed to provide information about the mechanical states. A general 
problem in such schemes is that the coupling needs to be sufficiently strong so that the exchange of quantum states occurs 
before decoherence sets in. In April 2010, the first experimental realization of a coupled quantum system involving a 
mechanical resonator (see Fig. has been reported by the groups of A. Cleland and J. Martinis [|2|. They demonstrated 
the superposition and coherent control of the quantum states of a mechanical resonator. The key aspect of their experiment 
is the use of a piezoelectric material for the resonator which boosts the coupling; it would be extremely hard to reach 
similar coupling strengths by using electrostatic forces alone (see also Table We will briefly come back to this ground 
breaking experiment in the second part of this Section when discussing coupling mechanical systems to other quantum 



systems from a more general point of view (see Sec. 5.2) 



Also concerning the linear detection schemes many challenges lie ahead. Strong coupling, ground-state preparation 
and detection imprecision below the quantum limit have all been reached now. These achievements set the stage for 
further studies on controlling and detecting non-classical states of mechanical motion; eventually one would like to have 
a mechanical resonator in its ground state which is strongly coupled to a single photon so that quantum state of single 
phonons and photons can be exchanged. Large coupling strengths and high Q-factors are necessary ingredients for reach- 
ing this goal. It is not a priori clear which detector scheme will be the best suited. The optical and superconducting cavities 
have the advantage that the underlying concepts are well known and have studied in detail; for detection schemes based 
on other (mesoscopic) devices (e.g. superconducting SQUIDs, carbon-based resonators) the understanding is clearly not 
at the same level. For example, the backaction mechanisms in these cases are not understood in all details as these cou- 
pled detector-resonator systems cannot be mapped directly onto the Hamiltonian describing cavity dynamics (Eq. [70| ). 
More theoretical and experimental research is needed to elucidate on the underlying mechanisms and physics. In the next 
subsection, we will briefly discuss some of the issues and challenges of linear detectors in more detail. 

5.1. NEMS as quantum-limited sensitive detectors 

It is clear that in the coming years more efficient read-out and cooling techniques will be employed. For optical systems 
challenges lie in increasing the coupling strength and to construct optical setups that can be incorporated in dilution 
refrigerators. In that respect, the recently developed on-chip optic experiments are promising. In the gravitational wave 



community the next generation of LIGO will use squeezed states to further improve the displacement sensitivity [336]. 



For the electronic systems, reduction of detector noise is a major issue which can be achieved by the implementation 



of quantum-limited amplifiers (SQUIDs p37\ |338[ ], point contacts [ |206|, |339[ |340|] or Josephson parametric amplifiers 
[341, 6^]). Furthermore, nonlinear (quantum) effects (see e.g. Refs. |342, 343, 344]) have so far received less attention 
and this may become an interesting research line. 

For cooling experiments high signal-to-noise ratios are important as well as high Q-factors. Since Q-factors are limited 
by material properties, mechanical oscillators of new materials (preferably in a crystalline form) will have to be fabricated 
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Figure 35: Left: a cantilever resonator containing normal metal rings to study persistent currents. In a magnetic field, the currents produce a torque on 
the cantilever which shifts the cantilever frequency. By measuring the resonance frequency, the persistent currents could be studied with unprecedented 
sensitivity. From A. C. Bleszynski-Jayich et at. Science 326 (2009) 2 72-275. R eprinted with permission from AAAS. Right: Schematic representation 



of the diamond-based scanning spin microscope proposed in Refs. [345, 34f , 347[ |. The single spin of a NV defect in a diamond tip is a sensitive 



magnetometer for the local magnetic field with nanoscale spatial resolution. Magnetic fields near the surface shift the electron spin resonance (EPR) 
frequency of the NV center, which can for example be detected by exciting the EPR transition with a microwave field and monitoring the change in 
photoluminescence of the probe spin. Reprinted with permission from C. L. Degen, Appl. Phys. Lett. 92 (24) (2008) 2431 IL Copyright 2008, 
American Institute of Physics. 



and characterized. We should add here that the understanding of the mechanisms controlling dissipation in mechanical 
resonators is also still open. For example, the role of tension in elevating the Q-factor is not completely understood. The 
temperature dependence of the Q-factor is another open-standing problem; in some cases a very strong dependence is 
observed, while in other cases the dependence is weak. What limits dissipation and what is the role of microscopic defect 
such as two-level fluctuators [|l6[ Furthermore, for low-mass resonators new effects may start to play a role such 

as the influence of adsorbents on the resonator surface and nonlinear damping terms associated with the induced tension 



introduced in Sec. 2.4.4 



The push for refining detection schemes will undoubtedly lead to the construction of better sensors, which eventually 
will be quantum-limited. These will lead to new applications aimed at both solving fundamental science questions as 
well as the development of commercial products. An example of the former direction is the recent detection of persistent 
currents in metal rings fabricated on nanoscale cantilevers [348], illustrated in Fig. |3^a. The precision is one order of 
magnitude better than previous experiments on persistent currents, which are mainly based on SQUIDs. It is also expected 
that the newly developed optics techniques (i.e., on-chip optics) will find their way in novel applications in controlling, 
stopping and storing light [ 349| , 350 1. Another promising example is the detection of single spins using magnetic resonance 
force microscopy analogous to the magnetic resonance imaging. Single electron spins have already been detected [351 1; 
the next step is the detection of nuclear spins which requires a thousand fold increase in the sensitivity. An interesting 
recent proposal (Fig. p5|b) is the use of NV centers so that operation at room temperature becomes feasible |347|. If 
successful, this would yield a revolution in imaging: Many chemical elements carry a nuclear magnetic moment, so that a 
sensitive enough detector can deter mine t heir identity and arrangement in more complex molecules. Finally, NEMS may 



354 PV, 355 1 : the state-of-the art is that single gold atoms can 



find an application in mass sensing [ [352| , [353 
be measured (see Sec. The goal is to be able to detect masses with resolution better than 1 Dalton (the mass of a 
hydrogen atom) so that each element can be identified making it a mass spectrometer [356|. 



5.2. Hybrid quantum mechanical systems 

In most of this review we have treated optical systems with movable parts as stand-alone systems. Interesting new 
quantum systems can be built when coupling optical set-ups to other quantum systems such as an ensemble of atoms in 
a Bose-Einstein condensate. The coupled system would form a hybrid quantum system in which hybrid strong coupling 
would enable the creation of atom-oscillator entanglement and quantum state transfer A theoretical proposal appeared in 
2007 1 357 1 an d in 2008 the interaction between a Bose-Einstein condensate and the optical field in cavity was studied ex- 
perimentally [358 1 demonstrating strong backaction dynamics. Subsequent theoretical work showed that strong coupling 
between a mechanical resonator and a single atom should be feasible [ 359 ] . Strong coupling with a single atom has indeed 
been observed in the blue-detuned cavity regime: Experiments with a single, trapped Mg^ ion interacting with two laser 
beams showed the stimulated emission of centre-of-mass phonons [360|. At high driving, coherent oscillating behavior is 
observed, which can be viewed as the mechanical analogue to an optical laser, i.e., a phonon laser In another approach, 
the coupling between vibrations of a micromechanical oscillator and the motion of Bose-condensed atoms on a chip are 
mediated by surface forces experienced by the atoms placed at one micron from the mechanical structure [361 j. As illus- 
trated in Fig. ^ the Bose-Einstein condensate probes the vibration of a cantilever By adjusting the magnetic trapping 
potential, the discrete eigenmodes of the condensate can be tuned, which in turn tunes the coupling to the cantilever 
Systems of many resonators coupled to each other via atoms in an optical lattice are envisioned [362]. 
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Figure 36: Bose-Einstein condensate (BEC) on an atom-chip couples to a cantilever. The cantilever motion modulates the potential of the atoms, thereby 
coupling the cantilever position to atomic motion. Left: schematic overview of the setup. The atoms of the BEC detect the cantilever vibrations which 
can independently measured via the readout laser. Right: photograph of the atom chip. Reprinted figures with permission from D. Hunger et al. , Phys. 
Rev. Lett. 104 (2010) 143002. Copyright 2010 by the American Physical Society 



A different approach to create entangled states of mechanical motion by coupling a mechanical quantum oscillator to 
a solid state device. In the experiment of the group of Cleland and Martinis [01, a Josephson qubit served as a two-level 
system probing the mechanical states of a 3 ng FBAR resonator (c.f. Fig. |34|). With a mechanical frequency of 6 GHz, 
a mechanical Q-factor of about 100 and a coupling of 62 MHz, there was just enough time to perform a few quantum 
operations on the mechanical resonator before the quantum state decohered. 

If it would be possible to increase the Q-factor while maintaining a strong coupling, quantum mechanical manipulation 
of long-lived phonon states comes into reach. For example, with a Q-factor of 10^ and a frequency of 1 GHz, the 
mechanical state would survive for 100 /is before it significantly decoheres. This time is typically much longer than the 
coherence time of a superconducting qubit, indicating that quantum information can be stored in the mechanical motion 
and transferred back when needed. The long phonon life times should also be contrasted to the short-Uved states in optical 
cavities: with a frequency of several hundreds of THz and a state-of-the-art optical Q of 10^, the photon lifetime is a 
few ns. Also in this case the mechanical states can be used as storage or delay unit for quantum information if efficient 
exchange of photons and phonons can be achieved (i.e., strong coupling). From a more general point of view and as a clear 
direction for the future, one can envision the construction of hybrid systems involving quantum mechanical oscillators 
that exploit and combine the strengths of the individual quantum systems. The mechanical quantum resonators would 
then serve as reservoirs of long-lived states with the advantage that they can be well coupled to a variety of other quantum 
systems. 

As a final note, let us return to the crossover between quantum and classical systems as discussed in the Introduction. 
The quantum experiments of Refs. [ p2[ , p3| ] were performed on mechanical oscillators with a mass of 3 ng and 48 pg 
respectively. This means that a mechanical object consisting of about 10'^ atoms can still behave quantum mechanically. 
The crossover has thus not been reached in experiments and continues therefore to be a subject of future studies. Again 
this aspect shows that the field of mechanical systems in the quantum regime is still largely unexplored: It will still rapidly 
develop in the years to come bringing many exciting new experiments and discoveries. 
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A. Complex Green's function and displacement 

One is often interested in knowing the amplitude A(f) and phase ip{t) of the displacement when it is written as u{t) - 
A{t) cos(a>iit + ifit)). Since this is only one equation for two functions, the amplitude and phase are not uniquely defined. 
Consider, for example, the solution A(f) - u{t)/ cos(w«f) and (p{t) = 0. Already for a pure sine wave u{t) = sin{ojRt) this 
results in a rapidly varying amplitude. The usual notion of the amplitude and phase are that when the signal is close to 
sinusoidal, the amplitude and phase are that of the sine wave. This can be realized using the complex displacement, which 
is defined as udt) = Ac expiicjRt) with the requirement that Re[Mc(f)] - u{t). A convenient way of implementing this, is 
using the complex extension of the resonator Green's function, cf. Eq. 

Ucit) = Kit) ® F{t)/kR, licit) = -ie''^"' ■ e-'^&it), (A.l) 

so that Re[/jc] = hno- The amplitude and phase of the resonator displacement are in this case given by the modulus and 
argument of the complex amplitude: A(t) = \Ac\ and ip(t) = ZAc respectively. 

Another notion that is often used in the literature is that of quadratures. Now the displacement is written as u{t) = 
U{t) cos(<jJsf) + P{t) smiuRt). Here U is the in-phase component and P is the out-of -phase component (sometimes called 
quadrature). The quadrature representation is related to the amplitude-and-phase representation by: U - A cos ip, P - 
Asini^, and by A'^ - + P^, tamp - P/U. The quadratures are also readily calculated using the complex Green's 
function: U = Re[Mf] and P = Im[Mc]. Finally, we note that if the displacement is in the in-phase quadrature, i.e. 
u{t) = A cos(wfif)> then the velocity du/dt « -corA siniaiRt) is in the out-of-phase quadrature. 



B. Optimal filtering of v(0 



In the presence of both position (imprecision) and force noise, one wants to reconstruct the resonator motion in the 
absence of the detector, M,(f) = /!«(?) ® F„{t), as good as possible from the measured time trace v(f) of the detector output. 
This is done by finding the estimator u - g(t) ® v(f) that minimizes the resolution squared: Au^ = E[(m,- - u) ^ ]. U sing 
the autocorrelation functions and converting these into noise PSDs using the Wiener- Khinchin theorem [200, 207 1 the 
resolution is written as: 

1 r°° 

Am2 = /;„_„,(0) - 2R,,i,(0) + RcdO) = — J [S „,„,(t^) - 2G(co)S „,„(w) + \G(co)\^S „] dco. (B.l) 



Minimizing this w.r.t. G, yields the optimal filter Gopt = S„,JSn, [200|, where: 

+ S,„,„+2A^Re[AAH'RyS^,„^,,„], 

so that the squared resolution is given by: 

Am^ = — [Su,u,iaj) - |Gopt(w)|25„,] dw. 
Depending on the properties of the detector and the coupling A, two important limits can be distinguished: 



(B.2) 
(B.3) 

(B.4) 



The detector exerts backaction force noise, but the displacement noise is negligible: - 0. In this case the 

integral in Eq. B.4 is easily solved and one finds: 



S Y' ( 

^ ^ <t'de,.„<t'de,.„ ) 



s 



(B.5) 



The resolution thus increases (gets worse) with increasing A as the influence of the detector, i.e. the force noise 
FBA,n{t), on the resonator motion grows. For small values of A this goes as Auba A. 

The detector ad ds d isplacement imprecision noise whereas the backaction noise is very small, i.e. 5(Drf„„Od<,,„ - 0. 
In this case Eq. B.4 reduces to: 



O S 

Aul = {u]) ■ J (SNR, Q) , with SNR = ^- 

"-R >-> ll„u„ 



(B.6) 
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Figure B.37: Plots of the function /(SNR, Q) for three different values of Q. The right side of the plot where SNR 2> 1 corresponds to the situation 
where the resonator peak is well above the detector noise level and the resolution is good. On the other hand, on the left side (i.e. for SNR <*c 1 ) the 
mechan ical signal is buried in the imprecision noise and the error is as large as the resonator signal giving 7 = 1. (b) An alternative definition of Ait^ 
(cf. Eq. p.Tp is the area of the gray rectangle indicated in the noise spectrum. 



under the assumption that the detector nois e PSD referred to the detector input 5 = 5,,,_v,,/A^|/l,,p is white. 
y(SNR, Q) is a function (see the plots in Fig. B.37 1) that depends on the quality factor and the signal-to-noise ratio, 
SNR. J tends to 1 when the signal-to-noise ratio is well below unity. In that case, the signal contains so much noise 
that it hardly contains information about the displacement and using the signal v(f) is not going to give much more 
information then just assuming that the resonator is at m = 0. The average error that is made in the latter case is 

When SNR — > oo the function 7(SNR, Q) goes to zero as Q/(2^^^ SNR^^"^). The resolution improves with increasing 
coupling as the resonator signal is amplified more and more with respect to the noise floor S . 



For practical purposes it is convenient to use a slightly different definition of the resolution [ [29 1[ |30| , |19|] that does 
not involve the signal-to-noise ratio dependent estimator: 



"-^""""2 2 ~ SNR' 



(B.7) 



which is readily extracted from the measured noise spectra as indicated in Fig. B.37b. This definition is based 
on the fact that one can measure the position during a time ~ QZ/r before the resonator has forgotten its initial 
amplitude and phase (see also Fig. [l^). Note again, that this is the imprecision noise of the detector which does not 
take the effect of backaction noise into account. 



C. Square-law detection 

In a M-squared detector the output of the detector depends quadratically on the displacement. In analogy with Eq. |o| 
the detector output is now v{t) = AAf\t) ® u{t)^ + v-^\t), where the superscript indicates that this is the response of the 
square-law detector. A resonator oscillating with frequency /«, results in a detector output with frequency components at 
2/r and at dc. As with every detector, the output also contains imprecision noise, v„ . Again, by increasing the coupling 



strength A the signal-to-noise ratio can be improved. In Sec. 3.4 it was shown that for a linear detector this leads to an 
increased backaction on the resonator, so it is interesting to see if the square-law detector has backaction. The backaction 
of a detector can vanish when the commutator of the quantity that is measured (m^ - u^{a^ + a)^) at different times is zero 
In the Heisenberg representation [||| we have [M^(fi), M^(f2)] = 4/Mq sin((yR{f2 - fi)) ■ (u(ti)u(t2) + u(t2)u(ti)) + for 
the square detector Quantum mechanics thus requires that the detector exerts a backaction force noise on the resonator to 
comply with the uncertainty principle. Also the quantum limit on the resolution for a square detector can in principle be 



derived using the formalism outlined in Sec. 3.4 



In many implementations, however, the detector has a narrow bandwidth and its output does not contain the signal 
at frequency 2/r, but only the dc component (see e.g. Ref. [|T|). The operator corresponding to the M^-detection can be 
written as = Mo(a^ + fl'^ + flfl' + a'a). The first two terms oscillate at twice the resonator frequency and when they 
are discarded, one finds that this detector detects aa' + a^a - 2h + I. In other words, it detects the number of quanta or 
the energy in the resonator [ [269| ] . Using the fact that d«/df - 0, one finds that [n(fi), n(f2)] = and that this detector can 
be backaction eva ding [p]|. The measurement of the resonator energy will change the phase of resonator, but leaves the 
energy unaffected | ]269| ] . 
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